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I*.  ARSTRACT 

The  purpose  of  this  project  was  to  formulate  possible  improvements 
in  the  design  codes  concerning  building  structures  which  are  subjected  to 
earthquake  loads.  In  addition,  suggestions  were  made  concerning  ways  to 
incorporate  the  concepts  and  methods  of  discrete  mechanics,  statistical 
analysis,  as  well  as  earthquake  behavior  of  concrete  and  metal  structures 
into  the  seismic  design  code. 

Available  literature  on  existing  design  codes  is  summarized  herein. 
Design  philosophies  and  methodologies  in  earthquake  engineering  were  also 
studied.  A  "design-tree"  technique  was  then  adapted  to  present  the  recom¬ 
mended  improvements  in  the  context  of  existing  specifications.  Moreover, 
approaches  to  certain  expected  problems  for  the  implementation  of  these 
inprovements  were  outlined  along  with  a  set  of  recommendations.  Furthermore, 
literature  reviews  and  some  original  contributions  in  discrete  mechanics, 
statistical  methods,  and  dynamic  behavior  of  structures  are  given  in  three 
appendices  as  backgroung  information. 


M.  KEY  WMOS 

Seismic  Design 


Earthquake  resistant  structures:  concrete  construction 

steel  construction 


TECHNICAL  REPORT  S-10 


SEISMIC  EESIGN  OF  BUILDING  STRUCTURES 


by 

J.  T.  P.  Yao 
C.  Ctoid 'varan 
A.  Gilrpinar 
C.  L.  Hulsbos 


research  performed  at 

Uhiversity  of  New  Mexico 
CONTRACT  No. :  DACA  23-70-C-0062 


July  1972 


Department  of  the  Army 

CONSTRUCTION  ENGINEERING  RESEARCH  LABORATORY 
P.O.  Box  4005 
Champaign,  Illinois  61820 

C 

Approved  for  public  release;  distribution  unlimited. 


ABSTRACT 


The  pxirpose  of  this  project  was  to  formulate  possible 
improvements  in  the  design  code  concerning  building  struc¬ 
tures  which  are  subjected  to  earthquake  loads.  In  addition, 
suggestions  were  made  concerning  ways  to  incorporate  the 
concepts  and  methods  of  discrete  mechanics,  statistical  anal 
ysis,  as  well  as  earthquake  behavior  of  concrete  and  metal 
structures  into  the  seismic  design  code. 

Available  literature  on  existing  design  codes  is  sum¬ 
marized  herein.  Design  philosophies  and  methodologies  in 
earthquake  engineering  were  also  studied.  A  "design-tree" 
technique  was  then  adapted  to  present  the  recommended  im¬ 
provements  in  the  context  of  existing  specifications .  More¬ 
over,  approaches  to  certain  expected  problems  for  the  imple¬ 
mentation  of  these  improvements  were  outlined  along  with  a 
set  of  recommendations.  Furthermore,  literature  reviews  and 
scmie  original  contributions  in  discrete  mechanics,  statisti¬ 
cal  methods,  and  dynamic  behavior  of  structures  are  given 
in  three  appendices  as  background  information.  A  siirple 
illustration  of  the  direct  ^jproacft  to  seismic  design  is  contained  in 
the  fourth  e^ipendix. 
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This  investigation  was.  performed  in  the  Department  of  the  Civil 
Eiigineering  at  the  University  of  New  Mexico  for  the  U.S.  Anty  Construc¬ 
tion  Engineering  Research  Laboratory  (CERL)  under  Contract  No. :  DACA 
23-70-C-0062.  This  research  was  performed  under  Project  4EM78012ADK1, 
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SEISMIC  DESIGN  OF  BUILDING  STRUCTURES 


PART  I :  INTRODUCTION 


General 

j  It  is  a  well  known  fact  that,  throughout  man's  history, 

human  life  and  property  have  been  lost  during  strong-motion 
earthquakes.  As  examples,  (a)  property  damage  from  the  1964 
Alaska  Earthquake  was  on  the  order  of  $300  million,  (b)  al¬ 
most  11,000  persons  lost  their  lives  during  the  1968  Iran 
Earthqueike  (1,  2)*,  and  (c)  over  40,000  persons  died  during 
the  recent  1970  Peru  Earthquake.  Detailed  descriptions  of 
the  damage  from  the  Alaska  Earthquake  as  well  as  other  re¬ 
cent  earthquakes  can  be  found  elsewhere  (3-7)  .  While  it  is 
not  possible  at  present  to  prevent  the  occurrence  of  strong-motion 
earthquakes,  ccntinous  efforts  have  been  made  to  inprove 
the  design  of  structures  in  order  to  minimize  the  earthquake 
dcunage  of  civilian  as  well  as  military  structures. 

Existing  specifications  concern?  j  earthquake  loads 
such  as  the  Uniform  Building  Code  (8),  Recommended  Lateral 
Force  Requirements  (9) ,  and  Seismic  Design  for  Buildings 
(10),  provide  a  set  of  equivalent  static  lateral  loads.  Al¬ 
though  these  codes  are  based  on  the  dynamic  analysis  of 
earl.iquakes  (11),  Blume  (12)  has  shown  that  structural  de¬ 
formations  induced  by  loads  based  on  code  requirements  may 
be  considerably  less  than  those  computed  for  the  same  struc¬ 
ture  subjected  to  the  excitation  of  a  recorded  strong-motion 
earthquake.  As  Blume  pointed  out,  the  safety  of  the  struc¬ 
ture  against  collapse  is  greater  than  that  indicated  by  the 
linear  cuialysis  because  a  considerable  amount  of  damping  is 
introduced  by  the  failure  of  "architectural  clothing"  as 
well  as  the  elasto-plastic  behavior  of  the  frame. 

Veletsos  and  Newmark  (13)  studied  an  elasto-plastic 
single-degree-of-freedom  system  subjected  to  earthquake  mo¬ 
tions  corresponding  to  the  1940  El  Centro  and  the  1933 
Vernon  earthquakes.  Penzien  (14)  investigated  the  response 
of  an  elasto-plastic  system  to  the  ground  motion  of  the 
1940  El  Centro  earthquake,  which  confirmed  the  thought  that 
the  maximum  displacements  are  less  in  elasto-plastic  sys¬ 
tems  than  those  in  corresponding  elastic  systems.  Since 
then,  nonlinear  and  yielding  seismic  structures  have  been 
studied  by  Jennings  (15) ,  Goel  and  Berg  (16) ,  and  Penzien 


*Parenthetic  numerals  indicate  references . 
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and  Liu  (17) .  The  results  of  these  studies  seem  to  confirm 
a  generally  accepted  design  philosophy  which  aims  at  an  elas 
tic  response  in  the  case  of  small  (and  frequent)  earthquakes 
and  to  permit  inelastic  response  in  the  case  of  strong- 
motion  (and  infrequent)  earthquakes.  However,  it  is  well 
known  that  plastic  deformations  can  cause  cumulative  damage, 
which  could  lead  to  low-cycle  fatigue  failure  of  structures 
(18). 


Due  to  the  random  characteristics  of  earthquakes,  sev¬ 
eral  probabilistic  models  have  been  proposed  to  represent 
earthquake  ground  motions.  Housner  (19)  represented  the 
earthquake  as  a  series  of  impulses  that  were  random  in  time. 
Thomson  (20)  showed  that  it  is  reasonable  to  use  white  noise 
to  describe  earthquake  motions.  Bycroft  (21)  then  proposed 
a  specific  white  noise  with  a  constant  spectral  density  and 
a  fixed  duration  for  the  representation  of  a  standard  strong 
motion  earthquake.  Bogdanoff,  Goldberg,  and  Bernard  (22) 
suggested  a  nonstationary  random  process  consisting  of  a  sxjm 
of  damped  sinusoids  to  represent  earthquake  accelerations. 

In  a  subsequent  paper,  this  model  was  modified  by  Goldberg, 
Bogdanoff,  and  Sharpe  (23)  .  Lin  (24)  proposed  a  nonstation¬ 
ary  random  process  resulting  from  a  filtered  shot  noise, 
which  has  become  a  popular  model  in  the  study  of  seismic 
structures  to  date. 

Stationary  random  processes  for  earthquake  simulation 
have  been  developed  by  Housner  and  Jennings  (25)  using  a 
digital  computer,  and  by  Ward  (26)  using  an  einalog  computer. 
Earthquake  ground  motion  is  generally  regarded  as  a  non-sta- 
tionary  random  process.  Digital  computer  simulation  of  non¬ 
stationary  earthquake  motion  has  been  performed  by  Ar  n  and 
Ang  (27) ,  Jennings,  Housner,  and  Tsai  (28) ,  Shinozuka  and 
Sato  (29) ,  Rascon  and  Cornell  (30) ,  Levy  (31) ,  Hou  (32)  ,  and 
Iyengar  and  Iyengar  (33) .  The  simulation  proposed  by  Rascon 
and  Cornell  (30)  attempts  to  account  for  most  of  the  pertin¬ 
ent  physical  characteristics  of  earthquakes .  The  results  of 
Jennings,  Housner,  and  Tsai  (28)  as  well  as  Iyengar  and 
Iyengar  (33)  contain  parameters  which  account  for  different 
earthquake  types.  Recently,  Wirsching  and  Yao  (34)  com¬ 
pleted  a  Monte  Carlo  study  of  bilinear  seismic  structures 
using  analog  computation. 

The  distribution  of  earthquake  occurrences  in  time  and 
space  has  been  considered  by  Rosenblueth  (35) ,  Benjamin  (36) 
Cornell  (37) ,  and  Shinozuka  (38)  .  Benjamin  (36)  used  Bayes 
theorem  and  computed  the  posterior  probabilities  of  the 
earthquake  occurrences.  Cornell  (37)  showed  that  the  earth¬ 
quake  intensity  is  a  function  of  distance  from  a  fault.  Re¬ 
cently  Wirsching  and  Yao  discussed  the  probability  distribu- 


2 


tion  of  structural  response  to  simulated  earthquakes  (39) . 

The  traditional  techniques  for  dynamic  analysis  of  highly 
redundant  framed  structures  can  be  referred  to  as  "piece  by 
piece"  methods  of  analysis,  which  call  for  tedious  and  inde¬ 
pendent  dynamic  analysis  of  each  structure  in  the  sense  that 
any  variation  in  the  framework  properties  and  applicable 
parameters  requires  the  development  of  a  new  solution.  Re¬ 
cent  developments  in  the  field  of  mechanics  makes  it  feasible 
to  develop  explicit  f\inctional  solutions  and  formulas  for  dy¬ 
namic  characteristics  of  frames,  appliccible  to  prescribed  pat¬ 
terns  and  formations  with  arbitrary  number  of  bays  and  levels . 
The  classical  literature  concerning  functional  field  solution 
to  structures  is  due  to  Bleich  and  Melan  (40)  .  The  develop¬ 
ment  and  application  of  the  technique  to  dynamic  analysis  of 
framed  structures  is  due  to  Wah  (41) ,  Leimback  and  McDonald 
(42) ,  Ellington  and  McCallion  (43) ,  and  Dean  and  his  col¬ 
leagues  (44) .  The  technique  employs  the  calculus  of  finite 
differences  cuid  discrete  mechanic  concepts  to  derive  the  dif¬ 
ference-differential  equations  describing  the  dynamic  be¬ 
havior  of  framed  structures.  It  results  in  functional  solu¬ 
tions  which  are  exact  within  the  scope  of  linear  beam  theory 
euid  yet  describes  the  dynamic  behavior  of  various  types  of 
framed  structures. 

Although  the  existing  seismic  design  codes  are  based 
somewhat  on  a  dynamic  analysis,  the  use  of  a  set  of  "equiv¬ 
alent"  static  loads  does  not  result  in  a  solution  which  fully 
describes  the  dynamic  behavior  of  structures  under  earthquake 
loading  conditions.  The  practice  of  employing  an  "equivalent 
static  loading"  is  known  as  a  compromise  solution  with  trace¬ 
able  deficiencies.  It  was  recommended  in  lieu  of  the  com¬ 
plexity  of  a  dynamic  analysis  and  limitations  of  the  state  of 
the  art.  Moreover,  the  statistical  characteristics  of  earth- 
c[uakes  were  totally  disregarded.  Results  of  many  investiga¬ 
tions  concerning  the  dynamic  behavior  of  framed  structures 
are  now  available.  These  and  field  solution  techniques  are 
promising  as  a  supplemental  to  statistical  methods  in  serv¬ 
ing  to  provide  a  rational  basis  for  development  of  new  and 
comprehensive  codes  and  specifications. 


Objective  and  Scope 

The  objective  of  this  research  program  was  to  formulate 
possible  improvements  in  the  code  specifications  for  the 
seismic  design  of  buildings.  Available  literature  on  exist¬ 
ing  design  codes  was  reviewed  and  summarized  in  PART  II.  In 
PART  III,  design  philosophies  and  methodologies  in  earth- 
quedce  engineering  were  studied,  and  the  recommended  imp?  ve- 
ments  are  presented  with  the  use  of  a  "design-tree"  te-rnique. 
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Moreover,  ^jproadies  to  sene  problems  resulting  from  mplemantation  of  these 
inproveroents  are  outlined.  In  PART  IV,  ocnclusions  and  reoernnendatiens  are 
made.  An  introduction,  cis  well  as  seme  originad  solutions,  to  discrete 
mechanics  is  inclxided  in  i^spendix  A  ;  a  literature  review  of  statistical 
methods  in  earthquake  engineering  is  given  in  Appendix  B ;  tbe  seianic 
bdiavior  of  ooncrete  and  metal  structures  is  sunmarized  in  appendix  C ;  and 
an  illustration  for  the  direct  appa^oach  portion  of  the  design  ^ree  is  given 
in  i^jpendix  D.  . 


PART  II:  EXISTING  SEISMIC  DESIGN  CODES 


General 

Until  the  1940 's,  the  forces  induced  by  earthquake  ac¬ 
celerations  were  treated  as  static  laterial  forces  in  build¬ 
ing  codes  all  over  the  world.  The  equivalent  static  force 
was  simply  calculated  by  assuming  a  rigid-body  structure. 
During  these  past  three  decades,  engineers  have  adopted  more 
rational  and  refined  methods  to  calculate  these  forces. 

In  the  United  States,  the  first  attempt  to  better  under¬ 
stand  the  dynamic  properties  of  building  structures  was  made 
after  the  Long  Beach,  California,  earthquake  of  March  10, 
1933,  w^en  the  United  States  Coast  and  Geodetic  Survey  mea¬ 
sured  the  natural  periods  of  vibration  of  212  buildings  in 
two  major  directions  (45)  .  In  fact,  these  data  have  re¬ 
mained  the  empirical  basis  for  computing  natural  periods  of 
building  structures  to  date  (46) .  With  the  developnent  of 
high  speed  computers,  a  new  era  began  for  the  design  and 
construction  of  buildings.  Closely  connected  with  this  de¬ 
velopment  was  the  study  of  dynamic  response  of  building 
structures  to  earthquake  loads.  Any  techniques  which  could 
be  adapted  to  computer  methods  have  been  emphasized  greatly. 
Of  these,  matrix  methods  are  the  most  popular.  Either  a 
stiffness  or  a  flexibility  matrix  is  constructed  for  the 
total  structure  and  the  dynamic  responses  of  the  structure 
are  thus  calculated.  One  of  the  more  serious  limitations 
of  the  application  of  these  methods  in  earthquake  engineer¬ 
ing  was  the  fact  that  they  are  applicable  only  to  linear 
elastic  systems.  Therefore,  the  results  were  astonishing 
at  first.  The  responses  calculated  analytically  by  using 
these  methods  were  much  larger  than  those  recommended  by  the 
building  codes  (12,  4).  However,  later  it  was  understood 
that  the  difference  was  due  to  the  nonlinear  and  inelastic 
behavior  of  structures  during  strong-motion  earthquakes. 

Earthqucikes  Regulations  are  given  in  Sec.  2314  of  the 
Uniform  Building  Code  (8) .  These  regulations  are  taken  from 
the  recommendations  of  the  Seismology  Committee  of  the  Struc¬ 
tural  Engineering  Association  of  California  which  dates  back 
to  1957  (9).  Considering  the  amount  of  research,  and  the 
volume  of  literature  which  has  been  published  on  the  subject 
of  earthquake  engineering  since  that  time,  it  is  necessary 
to  look  at  these  recommendations  more  closely  and  compare 
them  with  some  of  the  other  building  codes  around  the  world. 
In  the  following,  recommendations  of  the  Uniform  Building 
Code  concerning  different  forces  and  requirements  will  be 
reviewed,  critiqued,  and  compared  with  recommendations  of 
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other  existing  building  codes . 


Recommendations  for  oomputing  the  equivalent  lateral  force, 
V,  are  a  very  importcint  part  of  the  earthquake  design  regula¬ 
tion,  because  the  seismic  coefficient  relates  the  weight  of 
the  building  to  the  equivalent  lateral  force.  In  all  the 
building  codes  which  are  reviewed  herein,  the  relationship 
between  the  weight  of  the  building  and  the  total  lateral 
force  can  be  characterized  by  a  simple  formula  as  follows : 

V  =  SW  (1) 

where 

V  =  Total  lateral  force 

W  =  Weight  of  the  building 

S  =  Seismic  coefficient 

However,  the  seismic  coefficient  S  contains  different 
factors  in  various  building  codes. 


Uniform  Building  Code  (8) 


In  the  Uniform  Building  Code,  the  relationship  is  given 

by/ 

V  =  ZKCW  (2) 

hence  we  have 

S  =  ZKC  (3) 

where 

C  =  the  dynamic  factor  which  depends  on  the  period  of 
the  building. 

K  =  a  factor  which  varies  for  different  types  of  build¬ 
ings  and  indicates  their  ductility. 

Z  =  a  factor  which  depends  on  the  earthquake  zoning  of 
the  locations  of  the  building. 

With  these  factors  in  mind,  it  is  appropriate  to  divide 
this  section  on  "Equivalent  Lateral  Force"  into  four  parts. 
The  first  three  will  discuss  the  factors  C,  K,  and  Z  sepa¬ 
rately;  the  fourth  will  be  devoted  to  the  distribution  of 
this  lateral  force  along  the  height  of  the  building. 

1.1  The  factor  C  depends  on  the  natural  period. of  vibration 
of  the  building,  T,  with  the  following  relationship: 
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r 

1 


r  =  0-05 

with  the  restriction  that. 


(4) 


C  <  0.10  (5) 

Therefore  the  computation  of  the  natural  period  of  vibration 
of  the  building  requires  attention;  hence  the  discussion 
will  start  with  an  analysis  of  T. 

In  the  cUssence  of  technical  data  and  rigorous  calcula¬ 
tion  the  Uniform  Building  Code  recommends  the  following  for¬ 
mulas  for  T: 


T  = 

O.IN, 

for  ductile  frames 

(6) 

T  = 

O.OSh 

n 

i 

otherwise 

(7) 

D 


where 

h^  =  height  of  the  building  (in  feet) 

D  =  dimension  of  the  building  in  a  direction  parallel 
to  applied  forces  (in  feet) 

N  =  number  of  stories 

T  =  period  of  vibration  (in  sec.) 

Following  the  Long  Beach,  California,  earthquake  of 
March  10,  1933,  the  United  States  Coast  and  Geodetic  Survey 
measured  the  period  of  212  (pre-1940)  buildings  (45)  ,  and 
obtained  Equations  6  and  7  on  this  basis.  However,  the  form 
of  the  formulas  has  some  theoretical  justification  as  well 
(46) .  Equation  6  is  valid  only  for  a  building  in  which  the 
lateral  resisting  system  consists  of  a  moment  resisting 
space  frame,  which  resists  100%  of  the  required  lateral 
forces  and  the  frame  is  not  inclosed  nor  adjoined  by  more 
rigid  elements,  which  tend  to  prevent  the  freune  from  re¬ 
sisting  late:.'ai  forces.  Equation  6  does  not  involve  the 
depth  D  nor  the  breadth  B,  because  (a)  the  stiffness  of 
each  floor  is  a  function  of  the  number  of  columns  on  that 
floor  which  in  turn  is  proportional  to  the  floor  area  BD, 

(b)  the  mass  of  each  story  is  also  proportional  to  the 
floor  area,  and  (c)  the  two  effects  are  said  to  cancel  each 
other  in  the  calculation  of  the  natural  period  T  (46) . 

For  a  shear-wall  type  structure  with  equal  story 
heights  £uid  stiffnesses,  the  Rayleigh  method  of  computing 
the  natural  period  yields: 
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T  =  27r  (0.63N) 


(8) 


where 

M  =  mass  of  the  floor 
k  =  story  stiffness 

It  is  to  be  noted  that  the  story  stiffness  is  a  func¬ 
tion  of  the  product  of  wall  stiffness  per  unit  length  cuid 
the  wall  length. 

Therefore,  Equation  8  can  also  be  written  in  the  form: 

T  =  27r  /|r  llsJML  (9 

w  /d~ 

where 


=  wall  stiffness  per  unit  length. 

Consequently,  Equation  7  can  be  a  reasonable  formula  if 
the  ratio  M/k^  is  a  constamt  for  all  shear  wall  buildings . 
Some  other  suggestions  for  calculating  the  natural  period 
are  summarized  below  (46) . 


Assume  that  the  first  mode  of  vibration  controls  the 
motion  and  that  the  displacement  of  the  building  in  this 
mode  is  given  by: 

y  =  Cf(x)  sin  ^  (10) 


where , 

C  -  a  constcint 
f (x)  -  mode  shape 

Using  Rayleigh's  method  and  taking  the  static  deflection 
curve  as  the  mode  shape,  we  obtain. 


T 


where 


(11) 


y^  =  deflection  produced  by  a  Ig  laiteral  load. 

For  shear-wall  type  buildings,  the  stiffness  and  mass  are 
assumed  to  be: 

k  =  2D  k^  (12) 

M  =  BDm^  +  2hm^  (B  +  D)  (13) 
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vihexe, 

s  average  mass  of  wall  per  unit  area 
ntf  »  mass  per  unit  floor  area 
h  *  story  height 

niis  leads  to  an  equation  of  the  type: 

T  =  CjN  S 

vhere,  =  a  constant 

C|  =  a  constant 
H  »  total  height. 

For  shear-beam  buildings,  the  mass  is  assuned  to  be: 

M  =  ir^  BD 

For  the  stiffiiess  of  space  frame  buildings,  three  assvnpticns  can  be 
made.  The  story  stiffness  can  be  assuned  to  be  proportional  to  the 
wei^t  above  the  story.  This  leads  to: 


where 


=  effective  stiffiiess  of  the  n  story  in  a  building  of  N  stories 

k  =  stiffiiess  per  unit  area  of  the  n^  story 
n 

kj^  =  stiffness  per  uiit  aunea  of  the  story 

Using  the  relationship  for  the  displacement  of  n^  floor  subjected  to  a  unit 
load  and  Ecjuations  15,  16,  and  17,  we  have. 


T  =  21^^  (0.8  ^)  =  Cj  ^  H  {!£ 

If,  however,  the  assutption  is  that  the  stiffiiess,  k,  of  the  story  is  the 
same  and  is  proportimad  to  the  total  nuiber  of  stories,  i.e. , 

k  =  k  BCN  (19 

n 

Again  we  have. 


T  =  2Tf^j^  (0.63  vfj)  =  C3  viT  =  (^  ,/h. 


The  third  assunption  is  that  the  story  stiffiiesses  are  the  same  and 
indf^iendent  of  the  total  nuiber  of  stories. 
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Then, 

k  =  BD  (21) 

and 

T  =  C^N  =  C^H  (22) 

Considering  the  Uniform  Building  Code,  the  assumptions 
underlying  Equations  6  and  7  can  now  be  better  understood. 
From  Equation  9 ,  it  can  be  seen  that  Equation  7  would  give 
reasonable  results  only  if  the  mass  to  wall  stiffness  per 
unit  length  ratio  of  all  shear  wall  buildings  were  constant. 
Equation  6  assumes  the  relationship  expressed  by  Equation 
21.  Both  of  these  assumptions  are  very  restrictive.  In 
fact,  for  shear-wall  type  buildings.  Equation  7  as  given  in 
the  Uniform  Building  Code  gives  the  poorest  fit  to  actual 
data  comparing  with  Equations  18,  20,  and  22  (46).  It  is 
probably  too  much  to  expect  from  a  simple  relationship  to 
give  reasonable  results  for  a  variety  of  buildings.  There¬ 
fore,  recommending  different  empirical  formulas  for  differ¬ 
ent  assumptions  about  the  mass  and  stiffness  distribution 
of  the  building  should  lead  to  more  reasoneOsle  results. 

A  plausible  approach  to  this  problem  has  been  presented 
by  Salvador!  and  Heer  (48) ,  who  proposed  to  calculate  the 
natural  period  by  combining  the  effects  of  shearing  deforma¬ 
tion,  bending  deformation,  rocking  motion  as  well  as  founda¬ 
tion  translational  motion.  For  the  Alexander  Building  in 
San  Francisco,  their  calculations  for  the  natural  period  of 
vibration  of  the  building  had  errors  of  only  between  4.5%  - 
8.7%.  By  using  extreme  values  of  all  parameters  for  normal 
buildings,  the  authors  plotted  two  curves  which  enveloped 
all  the  possible  periods  for  these  structures .  This  method 
could  be  adapted  to  give  a  range  of  value  for  T  rather  than 
a  single  value,  and  thus  improve  the  validity  of  empirical 
formulas . 


Since  the  response  is  given  as  an  explicit  function  of 
the  period  of  the  structure.  Equation  4  suggests  the  use  of 
a  response  spectrum  technique  in  its  derivation.  Therefore 
the  following  background  information  on  response  spectrum 
technique  is  presented  (49) . 


For  a  multi-degrec-of- freedom  system  and  using  the 
elastic  modal  superposition,  we  have  the  following  govern¬ 
ing  equation: 


2  ^  * 

+  2„X„a) 

n  nnnn  nn  M* 

n 
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in  which. 


*  amplitude  of  mode  "n" 

=  damping  ratio  of  mode  "n" 

0)^  =  natural  frequency  of  mode  "n* 


also. 


2 

M  *  *  y  (bf  M. 
n  .  <Pin  r 

(24) 

N 

p  *  =  y  A.  p. 

n  'in  1 

(25) 

where , 


.th 


(|i.  =  n  mode  shape. 

=  i^  mass 
=  i^  force 


In  the  earthquake  response  problem,  let 

where , 


V  = 


V  (t)  =  ground  acceleration, 
9 


and 


N 


rn 


(26) 


(27) 


The  solution  for  each  mode  can  be  written  in  the  form  of  the 
Ouhamel  integral,  or  for  earthqueike  response. 


n 


“  Mr*“ 


n  n 


V^(t) 


where , 


V  (t)  =  /  V  (T)e"^‘*’^^”^^sinu)(t-T)dT 

“  0  ’ 

We  have  the  following  relationships  and  definitions: 
"displacement"  v(t)  «  ^  V(t) 


(28) 


(29) 


(30) 
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"base  shear,"  Q{t)  =  Kv(t)  =  Mw^vCt)  =  Ma>V(t) 

"spectral  density,"  S  =  V 

max 

c 

"spectral  displacement,"  ~  ^ 

"spectral  acceleration,"  S  =  wS 

>  Si 

The  effective  acceleration  is  given  by: 


"  “n  “n 

n 


Thus  the  force  at  level  "i"  is 


(31) 

(32) 

(33) 

(34) 

(35) 


Pi„(t)  .  M^V.„(t)  =  (t)  =  jjj  a.„V„(t)  (36) 

n 

Then,  the  base  shear  becomes, 

1—1  1—1  n  n 

The  maximum  base  shear  is  obtained  using  response  spectrum 
superposition : 


2  2 

.  _  _n  ,,  c  -_£!<? 

"  M*  “n*vn  ~  M*  “an 

max.  n  n 


(38) 


Since  the  maxima  of  each  mode  do  not  occur  simultaneously, 
the  following  approximation  is  made  to  find  the  total  maxi¬ 
mum. 

Omax  =  “i  *  «2  ^  •••  °N  <”> 

The  above  procedure  is  based  on  modal  superposition  and 
the  effects  of  all  modes  of  vibration  are  taken  into  account. 
However,  the  Uniform  Building  Code  formula  for  the  coeffi¬ 
cient  C  is  based  on  a  more  approximate  method  which  only  con¬ 
siders  the  first  mode  of  vibration.  It  also  assumes  that  the 
first  vibration  mode  shape  is  a  straight  line  from  the  base 
to  the  top.  For  this  method: 

N 

L  =  I  M.4>.  (40) 

i=l  ^  ^ 

N  2 

M*  =  I  (41) 

i=l  ^  ^ 
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which  results  in: 

<W  =  <«> 

There  is  little  apparent  resemblance  between  Equation  42 
euid  the  Uniform  Building  Code  formula  except  for  the  depen¬ 
dence  on  the  natural  period  (or  frequency)  of  the  building. 
This  is  a  result  of  empirical  considerations  in  the  deriva¬ 
tion  of  Equation  4.  Another  difference  is  that  Equation  42 
is  a  result  of  purely  elastic  analysis;  therefore  it  gives 
much  higher  values  than  that  specified  in  the  Uniform  Build¬ 
ing  Code.  A  good  discussion  on  how  to  incorporate  ductility 
and  therefore  the  energy  absorption  capacity  of  buildings  in¬ 
to  response  spectrum  techniques  can  be  found  elsewhere  (50) . 

It  is  to  be  noted  that  the  inelastic  action  which  is  taken 
into  account  in  the  coefficient  C  is  further  modified  by 
another  factor  K  depending  on  the  type  of  structure  accord¬ 
ing  to  its  ductility. 

1.2  The  factor  K  represents  the  effect  of  the  type  or  ar¬ 
rangement  of  the  resisting  elements  of  the  structure,  and  is 
an  indication  of  the  structure's  overall  ductility.  In  the 
Uniform  Building  Code,  K  varies  from  0.67  for  moment  resist¬ 
ing  ductile  space  frames  to  1.33  for  box  systems  and  3.00  for 
elevated  tanks.  In  a  dual  bracing  system  where  both  shear 
walls  and  a  moment  resisting  space  freune  (to  take  at  least 
25%  of  the  total  lateral  force)  are  designed  to  carry  the 
lateral  loads,  the  factor  K  is  taken  to  be  0.80.  For  all 
other  buildings  with  normal  ductility,  the  factor  K  is  con¬ 
sidered  as  being  unity. 

The  difference  between  the  elastic  analysis  and  the  Uni¬ 
form  Building  Code  recommendations  is  due  to  the  effects  of 
inelastic  action,  which  the  code  has  taken  into  account. 
Equations  4  and  5  already  include  a  "ductility  factor"  which 
reduces  the  maximum  response  of  the  structure.  This  "ductil¬ 
ity  factor"  is  chosen  for  a  normal  building,  where  the  lateral 
deflection  is  partly  shear  deflection  and  partly  flexural. 

It  is  seen  that  for  this  type  of  building  the  "ductility 
factor"  which  is  already  incorporated  in  C  is  untouched  (i.e. 

K  =  1.00) .  However,  when  other  types  of  resisting  arrange¬ 
ments  are  used,  then  this  "ductility  factor"  is  modified  by 
increasing  or  decreasing  K.  As  an  example,  for  a  ductile 
moment  resisting  space  frame,  the  ductility  factor  increases; 
therefore  the  K  factor  reduces  (to  reduce  the  base  shear)  to 
0.67.  For  specially  designed  dual  bracing  systems  in  which 
the  space  frame  can  take  25%  of  the  lateral  load,  K  is  equal 
to  0.80.  For  more  rigid  structures  such  as  box  systems,  the 
ductility  factor  decreases  and  K  increases  to  a  value  of  1*33. 
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Whei  the  structure  in  i^uesticn  is  esqaected  to  act  as 
statically  determinate  (e.g.,  elevated  tanks),  inelastic  ac¬ 
tion  results  in  collapse;  therefore  for  these  types  of  struc¬ 
tures,  K  assvunes  its  highest  value  of  3.00. 

1.3  The  numerical  coefficients  for  the  factor  Z  depend  upon 
the  seismic  zone  map  of  the  United  States  as  given  in  the 
Uniform  Building  Code.  For  locations  in  Zone  No.  1  'Z'  shall 
be  equal  to  one-fourth.  For  locations  in  Zone  No.  2  'Z' 
shall  be  equal  to  one-half .  For  locations  in  Zone  No .  3  ' Z ' 
shall  be  equal  to  one. 

Seismic  zoning  or  regionalization  of  the  United  States 
is  a  very  difficult  task  and  is  usually  based  on  two  kinds  of 
observations.  The  first  one  is  a  statistical  approach  where 
the  history  of  earthquakes  in  a  given  location  acts  as  a 
guideline  to  estimate  the  intensity  and  the  frequency  of 
future  earthquakes.  The  best  source  of  data  which  is  avail¬ 
able  concerning  seismic  activity  in  the  United  States  is  the 
observation  of  the  frequency  of  small  shocks  for  a  given  lo¬ 
cation.  Unfortunately,  a  speculative  correlation  between  the 
number  of  small  shocks  observed  in  a  location  and  the  maximum 
expected  intensity  of  a  future  earthquake  can  be  misleading 
(51) .  Because  the  interest  of  structural  engineers  lies  with 
ijie  maximum  expected  intensity,  the  use  of  small  shock  data 
in  earthquake  engineering  is  rather  questionable.  The  second 
source  of  observation,  on  which  the  seismic  regionalization 
Ccui  be  based,  is  geological.  Usually,  both  of  these  obser¬ 
vations  have  been  evaluated  for  seismic  regionalization. 

In  the  United  States  seismic  regionalization  was  at¬ 
tempted  in  1950  and  1951,  when  the  "Seismic  Probability  Map 
of  the  United  States"  was  prepared  by  Roberts  and  Ulrich 
(52,  53).  Richter  indicated,  however,  that  this  "was  not 
strictly  a  regionalization  map,  since  it  was  directed  to 
estimate  risk  rather  than  maximum  intensity"  (51) .  Again 
quoting  from  Richter,  "The  Seismic  Probability  Map  was  of¬ 
ficially  retired  in  1952,  as  'subject  to  misinterpretation 
and  too  general  to  satisfy  the  requirements  of  many  users.' 
This  action  was  not  taken  in  consequence  of  scientific 
criticism,  but  as  a  result  of  pressure  from  a  business  group 
interested  in  lower  reinking  in  their  community . "  It  is  in¬ 
teresting  to  note  that  this  map  is  still  being  published, 
with  a  few  revisions,  as  a  part  of  the  Uniform  Building  Code. 

1.4  The  Uniform  Building  Code  recommendations  for  the  dis¬ 
tribution  of  the  total  lateral  force  V  are  as  follows : 

F^  =  0.004  V  (43) 
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with  the  following  restrictions, 

^  0.15V  (44) 

and 


F  =  0,  for  ^  £  3 

(45) 

“s 

(V-F.  )W  h 

„  _  '  t'  X  X 

(46) 

^x  n 

where , 

=  force  at  the  top  level 

D  =  the  plan  dimension  of  the  vertical,  lateral 
force  resisting  system  in  feet 

F  =  force  at  level  "x" 

X 

=  weight  at  levels  "i"  and  "x",  respectively 

h^,  hj^  =  height  at  levels  "i"  and  "x",  respectively 

Equations  43  through  46  are  not  applicable  for  one-  or 
two-story  buildings,  for  which  the  distribution  is  assumed 
to  be  uniform. 


It  can  be  seen  from  these  above  formulas  that  with  the 
exception  of  the  top  story  where  the  force  is  assumed  to  be 
a  little  greater,  the  distribution  is  linearly  increasing 
from  bottom  to  top  when  the  weight  and  story  heights  are  the 
same . 


As  it  was  discussed  in  item  1.1,  only  the  first  mode  of 
vibration  was  considered  in  the  determination  of  the  factor 
C.  The  inverted  triangular  distribution  is  a  simple  exten¬ 
sion  of  the  first  mode  assumption.  However,  assuming  the 
shape  of  the  first  mode  to  be  a  straight  line  is  overly  re¬ 
strictive.  Even  though  the  foinnulas  are  to  be  used  for  any 
type  of  structural  system,  the  straight  line  assumption  can 
be  considered  to  be  good  only  for  a  combination  of  shear 
wall  and  space  frame  building.  For  all  other  types  of 
structural  systems,  the  first  mode  shape  would  definitely  be 
represented  by  a  curve.  Furthermore,  the  first  mode  ap¬ 
proximation  holds  only  for  relatively  short  buildings.  For 
tall  buildings  the  effects  of  the  second  and  even  the  third 
mode  of  vibration  should  be  included.  The  question  of  ac¬ 
curacy  of  the  first  mode  approximation  was  considered  by 
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Clough  in  1962  (54) ,  when  the  errors  generated  by  the  first 
mode  approximation  were  computed  for  four-  to  twenty-story 
buildings  subjected  to  1940  El  Centro,.  1949  Olympia,  and 
1934  El  Centro  Earthquakes.  The  following  conclusion  were 
reached: 

(a)  The  errors  given  by  all  approximations  increase 
with  the  number  of  stories  because  the  higher 
mode  contributions  become  more  significant  as 
the  period  of  vibration  increases . 

(b)  First  mode  approximation  shows  the  greatest  dis¬ 
crepancies  in  total  shear  envelopes;  because  the 
higher  modes  are  most  effective  in  developing 
shears.  First  mode  approximation  does  not  even 
represent  the  shape  of  the  envelope. 


2.  French  Aseismic  Code  of  1964  (55) 


The  seismic  coefficients  (for  masses  along  the  height  of 
the  building)  are  given  in  the  form  of  a  product  of  four  fac¬ 
tors  which  are  discussed  as  follows: 

2.1  The  intensity  coefficient,  a,  is  somewhat  analogous  to 
the  factor  2  of  the  Uniform  Building  Code  as  discussed  in 
item  1.3.  However,  the  intensity  coefficient,  a,  is  more 
closely  related  to  the  maximum  response  because  it  is  a  func¬ 
tion  of  the  Modified  Mercalli  Intensity  scale .  The  reference 
intensity  has  been  chosen  as  VIII,  so  that  for  the  design  in¬ 
tensity  VIII,  a  is  tcJcen  to  be  unity.  It  is  assumed  that  a 
doubles  whenever  the  intensity  becomes  one  degree  higher. 

This  assumption  is  generally  considered  to  be  on  the  conser¬ 
vative  side. 

2.2  The  response  coefficient  6  depends  on  the  natural  period 
of  vibration  and  the  damping  of  the  building,  and  is  expres¬ 
sed  by  the  following  three  empirical  formulas  depending  on 
the  assumed  amount  of  damping : 

(a)  For  common  buildings  and  houses  f "normal"  damping) 

3(T)  =  1^211  (47) 

Vt 

and 

0.05  £  6(T)  £  0.10 

(b)  Buildings  with  large  spans  and  lower  density  of 
partitions,  or  with  partitions  tied  only  in  places 
("medium"  damping) 
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and 


B(T) 


0.Q85 

Vt 


(48) 


0.065  <  e(T)  <  0.13 

(c)  Structures  reduced  to  a  skeleton  with  no  external 
friction  (low  damping) 

cm  = 

V? 

and 

0.06  £  e(T)  £  0.20  (49) 

The  curves  for  Equations  4?  and  48  are  not  exactly  re¬ 
sponse  spectra.  Alterations  have  been  made  in  order  to  take 
into  consideration  higher  modes  of  vibration  as  well  as  pos¬ 
sible  plastic  actions.  For  Equation  49,  the  curve  is  much 
closer  to  the  true  spectrum,  so  that  it  can  be  applied  to 
flexible  and  slender  structures,  which  require  more  explicit 
consideration  of  higher  modes  of  vibration. 

Obviously,  B (T)  is  a  counterpart  of  C  of  the  Uniform 
Building  Code  discussed  in  item  1.1.  It  can  be  seen  that  the 
French  Aseismic  Code  gives  comparatively  more  conservative 
values  to  B(T)  tham  the  Uniform  Building  Code  does  to  C.  In 
fact,  even  Equation  48,  which  is  for  the  highest  damping,  is 
still  more  conservative  (by  30%)  than  the  single  formula 
given  by  the  Uniform  Building  Code.  An  extreme  difference  of 
100%  is  obtained  in  the  case  of  a  slender  structure  with  low 
darning,  since  the  upper  limit  of  B (T)  is  0.20  for  Equation  49 
as  opposed  to  0.1  for  C. 

The  use  of  B (T)  seems  to  be  more  versatile  and  adaptable 
than  the  coefficient  C  in  the  Uniform  Building  Code.  The 
qualitative  consideration  of  damping  is  also  another  param¬ 
eter  recognized  by  the  French  Aseismic  Code.  The  effect  of 
damping  is  recognized  by  the  Uniform  Building  Code,  but  it  is 
not  explicitly  included. 

2.3  The  foundation  coefficient  reflects  the  effect  of  (a) 
the  soil  conditions  of  the  immediate  location  and  (b)  the 
characteristics  of  the  building  foundation.  It  varies  from 
0.8  for  deep  fovmdations  on  firm  rock  to  1.3  for  piling  in 
moist  ground.  Response  curves  corresponding  to  periods  lower 
than  0 . 5  sec .  are  reduced  when  the  soil  of  the  foundation  is 
a  part  of  a  very  large  geological  formation  of  soft  ground. 

It  is  to  be  noted  that  this  coefficient  is  completely  ignored 
by  the  present  Uniform  Building  Code. 
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2.4  The  distribution  coefficient  y  is  the  coefficient  of  the 
term  r(T)/g  in  Equation  50.  It  depends  only  upon  the  char¬ 
acteristics  of  the  structure  and  describes  the  variation  of 
the  response  with  respect  to  the  level  considered. 


For  a  system  with  concentrated  masses  M(Z)  and  distri¬ 
buted  masses  m(z) ,  the  seismic  coefficient  a (h)  similar  to 
the  factor  S  defined  by  'Equation  1  applicable  to  the  masses 
situated  at  a  height  h  above  the  base  of  the  structure  is 
given  as : 

.  ZM(Z)X(2)  +  /m(z)X(z)dz  . 

a(h)  =  Y  - : - 5 -  X(h) 


ZM(Z)X  (Z)  +  /m{z)X  (z)dz 


(50) 


in  which. 


Z,  z  =  ordinates  of  concentrated  masses  M  and  distri¬ 
buted  masses  m,  respectively 

X(h)  =  a  function  giving  the  shape  of  the  fundamental 
mode 


g  =  acceleration  of  gravity 

r (T)  =  the  maximum  response  (acceleration)  of  a  single 
mass  oscillator  of  the  same  period  T  and  damping, 
as  the  fundamental  mode. 

This  coefficient,  y,  shows  a  basic  difference  between, 
the  French  Aseismic  Co  and  the  Uniform  Building  Code.  In 
the  Unifoimi  Building  Code  there  is  no  "distribution  coeffi¬ 
cient"  and  the  distribution  of  the  lateral  force  is  prede¬ 
termined  by  an  inverted  triangular  shape  regardless  of  the 
height  or  the  type  of  structure.  The  distribution  coeffi¬ 
cient  y  in  the  French  Code,  however,  varies  from  floor  to 
floor  not  depending  on  a  predetermined  function  but  on  t)ie 
mass  at  each  story  and  the  shape  of  the  mode  (not  neces¬ 
sarily  linear)  under  consideration.  In  this  respect  the 
French  Aseismic  Code  gives  more  freedom  to  the  design  engi¬ 
neer  to  use  his  judgment  about  the  shape  of  the  mode.  In 
addition,  the  French  Aseismic  Code  also,  requires  the  con¬ 
sideration  of  the  first  three  modes  of  vibration  for  slender 
structures  with  low  damping.  This  also  seem  to  be  a  step  in 
the  right  direction  if  the  discussions  of  item  1.4  are  con¬ 
sidered.  I 


3 .  Stcmdards  and  Regulations  for  Buildings  in  Seismic 
Regions  of  the  USSR  1957  (STT 

The  "Standards  and  Regulations  for  Buildings  in  Seismic 
Regions"  is  a  long  cind  involved  document  and  the  entire  dooitent 
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vail  not  ne  discussed  nere.  unly  exerpts  trom  sections  I  and  III  will  be 
discussed.  However,  the  titles  of  all  sections  are  cited  as  follows  for 
additional  information: 

X .  Earthquake  Zone  of  an  Area  or  a  Building  Site  and 
Design  Rating  of  Buildings  and  Structures 

II.  Planning  Cities  and  Towns 

III.  Seismic  Forces  for  Residential,  Civic,  Industrial, 
and  Farm  Buildings  and  Structures 

IV.  Industrial  and  Civic  Buildings  and  Structures 

V.  Water  Works  and  Sewerage 

VI.  Highways  and  Railways 

VII .  .  Hydraulic  Structures 

VIII.  Rural  Structures 

IX.  Field  Work  and  Control  of  Seismic  Requirements 

Seismic  regionalization  on  the  basis  of  maximum  expec¬ 
ted  intensity  began  in  the  USSR  as  early  as  1933  (51) .  In 
1947,  Gorshkov  came  up  with  the  first  seismic  map  of  the 
USSR  (57) .  Both  statistical  amd  geological  data  were  used 
in  the  preparation  of  the  map,  but  the  map  was  criticized  by 
Gubin  6n  geological  grounds  (58)  .  However,  the  map  was 
adopted  as  it  was  with  only  minor  alterations  (59) . 

The  seismic  regions  are  established  according  to  the 
GEOFIAN  scale  which  is  similar  to  the  Modified  Mercalli 
scale.  The  different  degrees  of  intensity  are  defined  ac¬ 
cording  to  the  conplitude  of  a  standard  pendulum.  Therefore, 
the  divisions  between  different  intensities  in  the  GEOFIAN 
scale  are  more  definite  than  those  in  the  Modified  Mercalli 
scale. 

Both  the  Modified  Mercalli  scale  and  the  GEOFIAN  scale 
contain  twelve  divisions  and  they  both  descend  from  the 
Mercalli-Concani  scale.  The  GEOFIAN  scale  can  be  summarized 
in  Table  1. 

In  the  USSR  regulations,  buildings  are  classified  ac¬ 
cording  to  their  use  and  importance  as  well  as  their  loca¬ 
tion  with  respect  to  the  seismic  regionalization.  Buildings 
are  classified  according  to  their  use  into  four  groups .  As 
excunples.  Group  I  includes  very  large  and  important  build¬ 
ings,  large  radio  stations,  etc.  For  this  group  the  design 
intensity  is  taken  to  be  one  degree  higher  than  that  indi¬ 
cated  on  the  seismic  regionalization  map.  Group  IV  consists 
of  temporary  light  buildings.  For  this  group  the  design  in¬ 
tensity  is  given  to  be  six  regardless  of  the  location. 
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TABLE  1 


GEOFIAN  Scale  for  Earthquakes  (56) 


Division 

Qualitative 

Description 

of 

Earthquake 

Instrumental  Amplitude 

Range  Mean 

6 

Strong 

1.1 

to  2 . 0  mm 

1.6  mm 

7 

Very  strong 

2.1 

to  4 . 0  mm 

3.0  mm 

8 

Destructive 

4.1 

to  8 . 0  mm 

6.0  mm 

9 

Devastating 

8.1 

to  16.0  mm 

12.0  mm 

10 

Annihilating 

16.1 

to  32.0  mm 

24.0  mm 

11 

Catastrophic 

>  32.0  mm 

— 

12 

Greatly 

Catastrcphic 

—  —  — 

— 

The  design  seismic  force  at  any  point  with  a  mass 
given  by  using  the  following  formula: 


S  = 
X 


in  which. 


= 


Q  k  3C 

X  C  X 


IS 


(51) 


dead  load  plus  0.8  times  the  live  load  at  point  x 
(except  for  warehouses,  etc.,  where  the  live  load 
is  not  reduced) . 


k_  = 


B  = 


seismic  coefficient, 
k  =  1/40;  for  8,  k 

dynamic  coefficient. 

o  _  0.9 

P  .p 

and  0.6  <  B  <  3.0 


For  design  rating  7, 

=  1/20;  and  for  9  =  1/10. 


(52) 


a  coefficient  which  is  a  function  of  the  deforma¬ 
tion  curve  resulting  from  free  vibration,  and  of 
the  position  of  the  load  Q  within  the  structure. 


5 


X 


q(\)  1  Q.q(x.) 

--  3=1. 

I  Q.q^(x.) 

j=l  ^  ^ 


(53) 


where  q(X]^)  and  q  (xj)  represent  deflections  at  levels  k  euid 
j ,  re  spect ively . 


20 


For  normal  buildings,  only  the  first  mode  of  vibration 
need  be  considered.  For  slender  and  flexible  structures,  the 
first  three  modes  must  be  considered,  and  6  is  multiplied  by 


1.6. 


The  similarities  between  the  above  regulations  and  the 
French  Aseismic  Code  are  quite  obvious.  Therefore,  compari¬ 
sons  and  discussions  which  were  carried  out  in  section  2  will 
not  be  repeated  here. 

One  of  the  strongest  points  of  the  USSR  Regulations  is 
the  attention  given  to  seismic  regionalization.  The  use  of 
the  more  exact  GEOFIAN  intensity  scale  and  the  combined  sta¬ 
tistical  and  geological  data  have  resulted  in  a  more  scien¬ 
tifically  sound  seismic  regionalization. 

Another  new  item  in  the  USSR  Regulations  is  the  incor¬ 
porating  of  the  use  or  the  importance  of  the  building  into 
the  seismic  coefficient  kc*  It  seems  to  be  desirable  to  re¬ 
duce  the  risk  where  the  building  is  massive  and  important  and 
its  collapse  may  involve  a  great  consequence  in  the  loss  of 
life  and  money. 


4 .  Aseismic  Provisions  for  the  Federal  District,  Mexico  (60) 


According  to  the  subsoil  characteristics,  the  Federal 
District  is  divided  in  two  zones  in  accordance  with  the  chap¬ 
ter  on  foundations.  In  addition,  buildings  are  classified 
into  three  groups.  A,  B,  and  C,  depending  on  their  use  and 
importance.  Group  A  includes  government  municipal  and  public 
buildings,  hospitals,  museums,  schools,  stadiums,  etc.  Group 
B  consists  of  private  habitation  and  public  places  where  con¬ 
glomerations  of  people  are  uncommon  (office  and  industrial 
buildings,  restaurants,  gas  stations,  etc.).  Isolated  and 
unoccupied  buildings  are  in  Group  C.  Another  classification 
is  made  according  to  the  structural  characteristics .  Three 
types  are  thus  classified.  Type  1  calls  for  framed  struc¬ 
tures  in  which  the  freime  is  designed  to  take  50%  of  the  shear. 
Type  2  consists  of  structures  where  the  frame  is  designed  to 
take  25%  of  the  shear.  Type  3  includes  water  tanks,  chimney 
stacks,  one  column  structures,  etc. 

The  base  shear  coefficients  are  tabulated  for  group  B 
buildings.  According  to  the  zone  and  the  type  of  structure, 
this  coefficient  varies  between  0.15  and  0.04.  These  tabu¬ 
lated  values  are  to  be  multiplied  by  1.3  for  group  A  build¬ 
ings.  On  the  other  hand.  Group  C  buildings  (mostly  unoc¬ 
cupied,  isolated  buildings)  do  not  require  aseismic  design. 
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Three  different  types  of  cuialyses  are  suggested:  sim¬ 
plified,  static  and  dynamic.  The  selection  of  the  method 
depends  on  the  problem  at  hand  and  is  described  in  detail  in 
the  code.  The  simplified  method  includes  buildings  with 
height  to  length  ratio  of  not  more  than  1.5  plus  other  re¬ 
quirements.  The  static  analysis  also  assumes  the  inverted 
triangular  shear  distribution  along  the  height  of  the  build¬ 
ing. 


For  dynamic  analysis  a  coefficient  "a"  is  defined  as 
follows : 


(a) 

In  the 

zone  of  low  compressibility, 

for  T  ^ 

0.5,  a  =  1.0 

(54) 

for  T  > 

0.5,  a  =  ^ 

(55) 

(b) 

In  the 

zone  of  high  compressibility. 

for  T  £ 

1.0,  a  =  0.5(1+T) 

(56) 

for  2.5 

^  T  ^  1.0,  a  =  1.0 

(57) 

for  T  ^ 

2.5,  a  =  ^ 

(58) 

For  the  purpose  of  computation,  it  is  assumed  that  the  ground 
suddenly  undergoes  a  constant  acceleration,  which  is  equal 
to  (ag)  times  the  base  shear  coefficient.  It  is  to  be  noted 
that  the  dynamic  analysis  requires  the  use  of  the  root-mean- 
square  of  the  responses  corresponding  to  each  mode. 

One  of  the  new  items  in  this  code  which  has  not  been 
covered  so  far  in  the  codes  already  mentioned  is  the  sug¬ 
gested  three  methods  of  analysis  for  different  conditions. 
This  permits  the  mode  to  be  simple  for  simple  cases  and  elab¬ 
orate  and  not  too  simplistic  for  buildings  which  require 
special  attention.  Another  new  point  is  the  superposition 
of  the  modes  of  vibration  of  the  structure  by  the  root-mean- 
square  method.  This  criterion  is  based  on  probability  anal¬ 
yses  (61)  . 


5 .  Earthquaike  Resistant  Design  Provisions  of  Other  Countries 

In  the  preceding  sections  2  through  4 ,  earthquake  regu¬ 
lations  of  France,  the  USSR  and  Mexico  were  discussed  and 
compared  with  the  recommendations  of  the  Uniform  Building 
Code.  In  this  section  the  provisions  of  the  building  codes 
of  four  more  countries  (Romania,  Japan,  Turkey,  Israel)  will 
be  mentioned.  However,  a  lengthy  discussion  or  a  comparison 
is  not  undertaken  because  the  recommendations  of  these  codes 
are  qualitatively  covered  in  the  three  codes  which  have  been 
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nenticned  already.  The  main  reason  to  choose  the  codes  dis¬ 
cussed  in  sections  2  through  4  as  a  basis  for  comparison  with 
the  Uniform  Building  Code  was  due  to  their  uniqueness  as  well 
as  their  accessibility. 

The  following  will  then  be  a  summary  of  the  earthquake 
regulations  of  the  aforementioned  countries  concerning  the 
total  lateral  force  and  its  distribution. 


5.1  Romanian  Regulations  for  Aseismic  Design*  (62) 


The  shear  force  S  is  the  product  of  the  total  weight  Q 
and  the  seismicity  coefficient  c.  This  seismic  coefficient 
is  a  product  of  five  parameters. 

(59) 


=  kg-n; 


with  the  restriction, 
c  <  0.02 
where , 

k  =  indicates  the  seismic  intensity  based  on  the 
®  Modified  Mercalli  scale 

n  =  depends  on  the  soil  resistance 

6  =  a  dynamic  coefficient  indicating  the  response 
spectrum  (acceleration) 


(60) 


6  = 


0.9 

T 


(61) 


6  = 


a  correction  coefficient  compensating  for  the 
single  degree  of  freedom  assumption  and  given  by, 

h 


k^k 


6  = 


n 

I  Qj 
1  * 


h 

I  Qj 
1  * 


(62) 


where , 
n 
h 

'ik 


=  number  of  stories 
=  story  under  consideration 

=  the  deflection  of  the  multi-degree  of  freedom  sys¬ 
tem  for  the  considered  mode  of  vibration 


Q,  =  the  mass  of  the  considered  point  cf  deflection 


*Note : These  regulations  were  established  taking  into  con¬ 
sideration  the  different  codes  used  in  other  countries, 
especially  the  USSR  and  USA. 


23 


Ip  =  a.  measure  of  damping 

The  distribution  of  this  shear  force  is  given  by; 

Sk  =  s  (63) 

^  ®k^k 

where  is  the  deflection  of  the  point  k  under  the  Qj^  load¬ 
ings.  For  simplification,  qj^  can  be  taken  as  a  parabola. 

5.2  Some  Recommendations  from  the  Japanese  Building  Code  (63) 

The  seismic  force  is  calculated  by  multiplying  the  sum 
of  dead  and  live  loads  with  the  coefficient  of  horizontal 
force,  which  is  not  an  explicit  function  of  period  T.  It  is 
taken  to  be  0.2  for  buildings  up  to  16  meters,  then  0.01  is 
added  to  the  coefficient  for  each  additional  4  meters  above 
16  meters.  For  wooden  buildings  on  soft  soil,  water  tanks 
and  projected  chimneys,  this  coefficient  is  increased  to  0.3. 

A  reduction  coefficient  which  can  be  used  to  reduce  the 
coefficient  of  horizontal  force  is  designated  for  combination 
of  type  of  construction  and  the  soil  conditions. 


Table  2 

Reduction  Coefficients  in  Japanese  Code  (63) 


Type  of 

Const. 

Kind  of 

Ground 

Wood 

steel 

R/C  or 
Composite 

Kind  I* 

0.6 

0.6 

CD 

• 

o 

Kind  II** 

o 

• 

00 

o 

• 

00 

0.9 

*  Kind  I:  Ground  consisting  of  rock,  hard  sandy  gravel, 

etc.  classified  as  Terhary  or  older  strata  over 
a  consideredjle  area  around  the  structure. 

**Kind  II;  Ground  consisting  of  sandy  gravel,  sandy  hard 
clay,  loam,  etc.  classified  as  diluvial,  or 
gravelly  alluvium,  about  5  meters  or  more  in 
thickness,  over  a  considerable  area  around  the 
structure . 
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Another  teible  is  provided  to  reduce  the  coefficient 
further  either  by  0.8  or  by  0.9  depending  on  the  location. 
This  reduction  is  a  result  of  seismic  regionalization. 

5.3  Regulations  for  Construction  in  Regions  of  Hazard  in 
Turkey  (64) 

The  total  seismic  coefficient  C  is  made  up  of  four  dif¬ 
ferent  coefficients. 

C  =  C^a6\|j  (64) 

C  =  coefficient  of  seismic  regionalization  and  can 
take  on  three  different  values 

a  =  a  coefficient  of  foundation  soil  conditions 

3  *  coefficient  of  building  importance 

1))  =  dynamic  coefficient 


for  T  <  0.5, 


=  1.0 


for  T  >  0.5, 


ula: 


ip  <  0.3  (67) 

The  story  weights  are  calculated  according  to  the  form- 


W.  =  G.  +  nP. 

W.  =  total  weight  of  story  "i" 

=  dead  load  of  story  "i" 

=  live  load  of  story  "i" 

n  =  is  taken  to  be  1  for  public  buildings  and  0.5  for 
others . 


The  total  shear  is  to  be  distributed  according  to  the 
following  formula: 


F.  =  F  - 
1  n 


w^hi 


I  W.h. 

F^  =  force  at  story  "i" 
h^  =  height  of  story  "i” 


(69) 


5.4  .  Israel  Standard  Code  Seismic  Loads  on  Buildings 

(Proposed)  (65) 

The  total  horizontal  force  H  is  calculated  by: 

H  =  C(G  +  aP) 
where , 

C  =  depends  on  the  seismic  region  (there  are  two  re¬ 
gions)  ,  type  of  building  and  fovindation  soil  con¬ 
dition  (six  combinations  are  considered) 

G  =  dead  load 

P  =  live  load 

a  =  varies  from  0  to  1  depending  on  the  importauice  of 
the  structure 


Overturning  Moment 

The  following  are  the  reocmnendaticns  of  the  1967  tftiiform  Buildino 
Code  (without  revision)  oonoeming  the  occputation  of  the  overturning 
mcments.  For  the  base  of  the  structure  the  overturning  moment  M  is 
given  by: 


M  =  J(F^h^  + 


n 


T  0.5 
■ 

and 

0.3  <  J  <  1.00 


(70) 

(71) 


(72) 


The  overturning  moment  at  level  "x" : 


=  J  + 


(73) 

(74) 


where  all  the  symbols  have  been  defined  in  section  1  of  this 
Part. 


It  C2U1  be  seen  that  the  overturning  moments  are  calcu¬ 
lated  statically  from  the  story  shear  forces  reduced  by  the 
factor  Jjj.  The  reduction  factors  and  J  for  the  overturn¬ 
ing  moment  at  the  base  decrease  as  the  period  of  vibration 
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of  the  building  increases  (or  as  the  building  gets  taller 
and  more  slender) .  Assuming  that  T  =  O.IN,  where  N  is  the 
number  of  stories  of  the  building,  we  get  the  following  ap¬ 
proximate  relationships:  J  =  1.0  for  a  three-story  building, 
J  =  0.5  for  a  ten-story  building  and  J  =  0.33  for  an  eighteen 
story  building.  Justifications  for  this  reduction  for  more 
slender  buildings  are  qualitative  in  the  following  manner 
(47,  50).  For  low  buildings  (i.e.,  short  T)  the  assumption 
that  the  structure  acts  as  a  cantilever  is  a  good  one.  But 
the  application  of  this  concept  to  tall  buildings  is  over¬ 
conservative,  because  maximum  total  moments  on  a  section 
through  the  building  and  maximum  shears  on  such  sections  do 
not  occur  for  the  same  combinations  of  modal  components.  The 
formula  for  the  distribution  of  overturning  moment  provides 
simply  for  a  straight  line  variation  of  moment  from  M  at  the 
base  to  zero  at  the  top,  even  though  this  is  not  consistent 
with  the  static  moments  computed  from  the  lateral  force  dis¬ 
tribution.  However,  it  is  reasonable  for  slender  structures 
(50)  . 


Another  point  made  is  that  no  tall  building  has  been  ob¬ 
served  to  simply  "topple  over"  during  an  earthquake,  whereas 
for  shorter  buildings  this  is  not  an  uncommon  occurrence. 
However,  this  argument  cannot  be  considered  valid  because  the 
number  of  observations  on  tall  building  behavior  during 
earthquakes  is  still  limited.  In  addition,  the  fact  that  no 
tall  building  was  observed  to  have  "toppled  over"  does  not 
mean  that  the  overturning  moments  do  not  excite  a  different 
form  of  failure. 

The  following  comments  against  the  code  formulas  were 
made  by  Derrick  in  his  commentary  to  Reference  47.  "In  the 
SEAOC  proposal  the  reduced  values  Fj^  are  used  to  compute  M, 
the  overturning  moment.  The  true  value  of  M  is  the  sum  of 
the  moments  of  the  equivalent  inertia  forces  capable  of  pro¬ 
ducing  the  maximum  dynamic  distortion  at  maximum  response. 

The  results  of  both  F^  and  Fj^j^  act  at  approximately  2/3  of 
the  height  of  the  structure  and  from  the  ratio  Vn/V  the  F^h^ 
is  (Vn/V)  (Fjjj^hjj)  ,  where  Fjjn  denotes  the  distributed  shear 
of  nth  mode  on  xth  story.  Hence  to  evaluate  the  maximum 
overturning  moment  from  F^,  there  should  be  a  correction 
factor  increasing  the  F^h^  by  (Vn/V) .  The  SEAOC  proposal 
provides  a  correction  factor,  but  it  decreases  the  siammation 
by  J  s=  0.5/^Vt^  with  limits  1.0  to  0.33."  Derrick  further 
comments  that  the  factor  J  moves  the  position  of  IFx  to  a 
point  lower  than  what  the  triangular  distribution  assumption 
would  predict.  As  an  example,  the  triangular  distribution 
assumption  would  place  the  resultant  at  2H/3 .  However,  with 
the  use  of  J-f actor,  the  resultant  is  placed  at  2H/9 .  Der¬ 
rick  claimed  that  this  difference  was  not  adequately  ex- 
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plained . 


In  the  above  discussion  the  subscript  "n"  refers  to  the 
dynamic  analysis  response  and  therefore  Vn/V  d.s  the  ratio  of 
the  shear  value  obtained  from  a  dynamic  analysis  to  the  rec¬ 
ommended  SEAOC  value. 

In  an  analytical  study,  Bustamante  calculated  the  ratios 
of  dynamic  to  static  overturning  moments  (66)  .  Both  direct 
addition  and  the  root  mean  square  values  of  the  spectral 
superposition  were  computed.  These  ratios  were  then  compared 
with  the  J  factor  of  the  Uniform  Building  Code.  It  was  con¬ 
cluded  that  there  could  be  large  differences  between  the  com¬ 
puted  ratios  and  the  reduction  factor  of  the  SEAOC  code .  The 
errors  involved  could  be  as  large  as  100%  even  for  moderate 
periods  of  less  than  0.35  seconds.  However,  all  errors  were 
fovind  to  be  on  the  safe  side.  The  considerations  were  formu¬ 
lated  by  assuming  that  there  were  no  errors  in  the  estimation 
of  T.  As  far  as  the  distribution  was  concerned,  the  author 
concluded  that  "the  ratio  of  static  to  dynamic  moments  in¬ 
creases  pronouncedly  towards  the  top.  The  tendency  is  op¬ 
posite  the  one  concerning  shears ,  for  which  the  larger  fac¬ 
tors  of  safety  were  at  the  bottom."  This  study  by  Bustamante 
seems  to  confirm  the  qualitative  arguments  of  Clough  (47)  and 
Blume,  Newmark  and  Corning  (50)  . 

The  foreign  building  codes  which  were  cited  previously  do 
not  provide  any  regulations  with  regard  to  the  overturning 
moment.  Therefore  it  is  assumed  that  the  overturning  moments 
are  simply  calculated  from  the  story  shears  for  which  the 
provisions  of  these  codes  are  extensive.  One  exception  is 
the  "Aseismic  Provisions  for  the  Federal  District  of  Mexico" 
(60) .  But  even  in  this  code  the  provisions  for  overturning 
moment  are  restricted  to  the  "Static  Analysis,"  and  are  as 
follows:  "For  design  purposes,  overturning  moments  may  be 

reduced,  but  for  each  frame  or  group  of  resisting  elements 
it  shall  not  be  taken  smaller  than  the  design  shear  at  the 
elevation  considered  times  its  distance  to  the  center  of 
gravity  of  all  the  corresponding  masses  above." 


Torsion 


In  the  Uniform  Building  Code,  the  static  method  is  used 
to  calculate  torsional  moments.  Torsional  analysis  is  re¬ 
quired  for  all  buildings  (symmetrical  buildings  included) 
and  the  minimum  design  eccentricity  is  to  be  taken  as  5%  of 
the  maximum  dimension  of  the  building. 

Considering  the  importance  of  the  problem,  the  Uniform 
Building  Code  treats  it  very  lightly.  Four  references  will 


be  cited  below  about  the  problem  of  torsion  and  some  of  their 
important  observations  and  conclusions  will  be  quoted.  One 
area  of  unanimous  agreement  is  the  importance  and  the  poten¬ 
tial  danger  of  the  problem. 

The  following  observations  were  made  by  Housner  and 
Outinen  (67) : 

"It  has  been  shown  that  an  unsymmetrical  building 
undergoing  torsional  accelerations  can  be  expected 
to  sustain  stresses  in  the  more  flexible  wall  that 
are  higher  than  those  predicted  by  the  customary 
static  method  of  analysis.  If  the  objective  of  the 
design  is  to  keep  the  maximvim  stress  within  the  usual 
allowable  limits ,  a  correction  should  be  made  to  the 
results  given  by  the  usual  static  method  of  analysis . 

On  the  other  hand,  if  the  objective  of  the  design  j  3 
to  provide  a  certain  ultimate  strength,  the  relative 
rigidity  of  the  wall  is  not  so  important  a  factor  as 
is  the  ability  of  the  wall  to  absorb  energy." 

Similar  and  additional  conclusions  are  also  reached  by 
Bustamante  and  Rosenblueth  (68) . 


"(a)  Dynamic  eccentricity  may  exceed  statically 
computed  values.  The  excess  is  particularly  im¬ 
portant  when  the  polar  moments  of  inertia  are  close 
to  their  critical  values;  these  are  attained  within 
the  range  of  usual  characteristics  of  buildings. 

Dynamic  eccentricity  is  defined  as  follows: 


! 

i’ 


(75) 


where , 

=  torsional  moment  in  the  i^  mode 

=  shear  force  in  the  i^  mode  of  a  statically 
eccentric  structure , 


" (b)  A  rough  estimate  of  torsional  dynamic  effects 
in  multistory  buildings  can  be  obtained  from  the 
response  of  a  single  story  structure  with  similar 
characteristics . 


" (c)  Excessive  dynamic  eccentricity  due  to  close¬ 
ness  to  the  critical  polar  moments  of  inertia  can 
be  greatly  reduced  by  changing  the  critical  values . 
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This  may  be  accomplished  by  increasing  the 

stiffness  of  perimetral  frames." 

From  these  observations  we  can  conclude  that  the  "core" 
type  buildings  where  the  resisting  elements  are  not  at  the 
perimeter  of  the  structure  are  more  susceptible  to  failure 
due  to  torsion.  This  effect  should  be  incorporated  in  the 
recommendations  of  building  codes. 

Shiga  (69)  concurs  with  Bustamante  and  Rosenblueth  in 
the  equivalent  single  story  concept  and  further  adds  that 
"It  should  be  noted  that  when  the  torsional  rigidity  is  less 
than  the  lateral  rigidity  the  fundamental  vibration  is  rota¬ 
tional  cuid  torsion  is  serious;  when  the  building  is  not  rigid 
against  the  ground  motion,  the  torsion  becomes  more  serious 
and  vibration  pattern  becomes  more  complicated."  In  addition 
Skinner,  Skilton  and  Laws  add  that  close  torsional  and  trans¬ 
lational  periods  are  dangerous  (70) . 

In  nearly  all  the  foreign  codes  cited  before,  torsion  is 
passed  rather  lightly  and  in  none  is  a  dynamic  torsion  anal¬ 
ysis  suggested.  In  the  Aseismic  Provisions  for  the  Federal 
District  of  Mexico  (60) ,  the  torsional  eccentricity  is  to  be 
calculated  for  all  stories  and  the  design  eccentricity  is  to 
be  1.5  times  the  computed  value  plus  or  minus  0.05  times  the 
maximum  plan  dimension.  The  1.5  factor  is  based  on  Housner’s 
suggestions  for  single  story  buildings  (67) . 


Allowable  Stresses 


The  ACI  and  AISC  codes  provide  for  a  33-1/3%  increase  in 
the  allowable  stress  of  concrete  for  dynamic  loads  (wind  and 
earthquake)  in  Article  1004  (71) .  This  increase  in  the  al¬ 
lowable  stresses  (both  concrete  and  steel)  has  also  been 
adopted  by  the  Uniform  Building  Code.  Also  for  the  ultimate 
stren^  design  of  ductile  reinforced  ccncrete  nonent  resisting  space  frames 
in  seismic  zones  2  and  3  the  ultimate  capacity  U  is  given  by; 

U  =  1.40  (D  +  L  +  E)  (76) 

or 

U  =  0.90  D  +  1.25  E  (77) 

whichever  is  grei ier.  Where, 

D  =  dead  load 
L  =  live  load 
E  =  earthquake  load 
U  =  ultimate  capacity 
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The  building  codes  all  around  the  world  have  recognized 
the  increase  in  the  allowable  stresses  of  different  materials 
under  temporary  loads  such  as  wind  and  earthquakes.  However, 
the  magnitude  of  the  increase  varies  greatly. 

In  the  USSR  a  service  factor  of  1.4  is  given  for  steel 
and  wooden  structures, 1.2  for  reinforced  concrete  and  1.0  for 
prestressed  concretes.  A  service  factor  of  1.4  is  equivalent 
to  increasing  the  design  strength  of  the  member  by  40%  (from 
translator's  note  in  Reference  56). 

In  Turkey  the  increase  in  the  allowable  stresses  depends 
on  the  foundation  soil  conditions.  The  increase  is  taken  to 
be  50%  for  good  conditions  and  30%  for  fair  foundation  soil 
conditions.  If  the  conditions  are  poor,  no  increase  is  per¬ 
mitted  (64) .  The  Japanese  Building  Code  allows  an  increase 
in  the  cillowabl'3  stresses  of  50%  for  reinforcing  steel  and  100%  for 
concrete  in  reinforced  concrete  (63) .  Also,  increases  of  50%  for  structural 
and  reinforcing  steel  and  timber,  and  33%  for  concrete  are 
recognized  by  the  Mexican  provisions  (60) . 

An  interesting  treatment  of  this  problem  is  in  the 
French  Aseismic  Code.  There  is  an  option  of  choosing  one  of 
the  two  methods  suggested.  The  so-called  "ordinary  method" 
increases  the  allowable  stresses  for  steel,  reinforced  con¬ 
crete  and  prestressed  concrete  by  50%.  The  second  method 
is  the  "state  limit"  method.  State  limit  is  defined  as  the 
state  beyond  which  the  construction  can  no  longer  remain  ap¬ 
propriate  for  normal  use.  This  corresponds  to  the  elastic 
limit  for  steel  structures  and  to  the  yielding  of  the  steel. 

The  seismic  forces  corresponding  to  the  nominal  design  in¬ 
tensity  Ifj  are  considered  as  being  normal  loads,  and  those 
corresponding  to  1^  +  1/2  as  being  exceptional  ones .  In 
addition,  for  various  types  of  loads,  different  "ponderation 
coefficients"  are  assigned  to  the  capacity  of  the  structure 
for  inelastic  action.  This  way,  a  more  uniform  design  can 
be  obtained  (55) . 

In  all  cases,  the  increases  in  the  allowable  stresses 
apply  only  to  the  combination  of  gravity  and  lateral  loads . 
Also,  they  are  given  for  working-stress  designs.  In  fact, 
this  is  the  basic  reason  for  the  inconsistencies  in  the  in¬ 
crease  permitted  for  the  allowable  stresses  by  the  above 
cited  codes.  Blume,  Newmark  and  Corning  explain  these  dif¬ 
ferences  as  follows  (50) ; 

"It  can  be  shown  that  the  reserve  capacity  of 
structural  members  to  resist  severe  earthquake 
loading  will  depend  upon  the  additional  stress 
requirement  to  cause  yielding  and  failure.  Thus 
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a  bracing  member  proportioned  only  for  earthquake 
stress  (no  dead  or  live  load  participation)  will 
tend  to  have  less  reserve  strength  and  energy 
capacity  than  a  member  the  .  size  of  which  is  de-. 
termined  by  dead  and  live  load  requirements  as 
well  as  earthquake.  The  result  is  an  inconsis¬ 
tency  in  the  strength  provided  for  different  mem¬ 
bers  in  the  same  structure,  when  a  working  stress 
or  even  a  modified  working  stress  design  proced¬ 
ure  is  used.  A  design  procedure  based  on  ultimate 
strength  or  energy  capacity  would  avoid  these  in¬ 
consistencies  . " 


Other  Considerations 


1.  Lateral  Force  on  Parts  of  Buildings  and  Appendages 

The  Uniform  Building  Code  formula  for  the  "Lateral  force 
on  parts  or  portions  of  buildings  and  other  structures"  which 
is  designated  by  Fp  is  as  follows : 

=  ZC  W  (78) 

P  P  P 

where 


the  weight  of  the  portion  under  consideration 


C 


P 


a  coefficient  given  by  Table  23-D  of  the  Uniform 
Building  Code  and  which  varies  from  1.00  for  "canti¬ 
lever  parapet  and  other  cantilever  walls  (except  re¬ 
taining  walls) "  and  "exterior  and  interior  ornamen¬ 
tations  and  appendages"  to  0.10  for  "tank  plus  ef¬ 
fective  mass  of  its  contents,  when  resting  on  the 
ground . " 


This  formula  is  analogous  to  the  base  shear  formula,  but 
the  coefficient  Cp  is  always  equal  to  or  greater  than  C. 

Since  there  is  no  dependence  of  Cp  on  the  period  of  vibration 
Equation  78  is  a  result  of  an  equivalent  static  analysis. 

The  problem  of  appendages  on  buildings  is  often  a 
serious  one.  The  attached  portions  are  subject  to  the  ac¬ 
celeration  of  the  building  which  in  general  is  much  larger 
than  ground  acceleration.  Normally  all  such  appendages  have 
different  periods  of  vibration;  therefore  a  rigorous  analysis 
is  generally  very  difficult.  Another  factor  is  that  usually 
these  appendages  are  statically  determinate  and  therefore 
cannot  develop  much  ductility,  which  is  desirable  in  buildings 
designed  to  resist  earthquakes.  This  is  vhy  a 
factor  which  is  equivalent  to  K  in  Equation  2  has  been  omit¬ 
ted  from  Equation  78. 
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The  following  conclusions  are  drawn  about  appendage  be¬ 
havior  during  earthquakes  by  Penzien  and  Chopra  (72) . 

(a)  The  two-degree-of-freedom  method  of  analysis 
quite  accurately  predicts  the  maximum  dynamic  re¬ 
sponse  of  an  appendage  attached  to  the  top  of  a 
multi-story  building  even  when  its  period  of  vibra¬ 
tion  coincides  with  the  fundamental  period  of  the 
building;  thus,  the  availed>ility  of  two-degree-of- 
freedom  response  spectra  makes  this  method  of  anal¬ 
ysis  practical. 

(b)  The  single-degree-of-freedom  method  of  analysis 
becomes  considerably  in  error  when  the  period  of  the 
appendage  is  near  the  period  of  one  of  the  lower 
building  modes;  therefore,  this  method  of  analysis 
should  not  be  used  in  such  cases. 

(c)  To  greatly  reduce  the  seismic  forces  in  an  ap¬ 
pendage,  it  should  be  designed  so  that  its  period 
of  vibration  differs  considerably  from  the  first 
mode  of  vibration  of  the  building  and  also  does  not 
coincide  with  other  lower  building  modes. 

(d)  The  seismic  forces  developed  in  an  appendage, 
even  when  designed  in  accordance  with  the  recom- 

.  mendations  of  (c)  above,  are  larger  than  code  values; 
therefore,  the  desirable  effects  of  inelastic  defor¬ 
mations  must  be  considered  as  a  standard  practice  in 
the  design  of  buildings. 

With  regard  tw  this  point  there  is  not  much  that  the 
other  international  codes  would  contribute  since  their  pro¬ 
visions  are  generally  in  the  same  form  as  that  of  the  Uni¬ 
form  Building  Code.  The  Turkish  Code  gives  a  coefficient  as 
a  fijnction  of,  in  fact  three  times,  the  coefficient  C  for 
the  whole  building  (64) .  The  code  of  Israel  specifies  a  co¬ 
efficient  Cp  which  can  be  as  much  as  twenty  times  the  coef¬ 
ficient  C  or  the  whole  building,  although  the  relationship 
is  not  explicit  (65).  The  USSR  provisions  specify  the  prod¬ 
uct  (see  section  3  on  page  18  for  definitions)  to  be  at 
least  five  for  portions  of  buildings  and  attached  structures. 

2.  Effect  of  Vertical  Seismic  Force 


In  most  building  codes,  including  the  Uniform  Building 
Code,  the  effects  of  the  vertical  seismic  force  are  totally 
neglected.  Generally  the  reasons  given  for  ignoring  the 
vertical  acceleration  are:  (a)  the  large  safety  factor 
(1.65)  which  the  building  is  designed  for  in  the  vertical 
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direction,  (b)  the  general  acceptance  of  the  fact  that  the 
vertical  acceleration  is  considerably  less  than  the  hori¬ 
zontal  accelerations,  and  (c)  rareness  of  failure  due  to 
vertical  acceleration  from  observations  at  the  field.  All 
these  reasons,  however,  are  delusive  and  at  best  simplistic. 
It  is  true  that  most  building  structures  are  very  stiff  in 
the  vertical  direction  and  for  this  reason  the  dyijamic  ampli¬ 
fication  in  the  vertical  direction  is  almost  negligible. 
Therefore,  the  axial  stresses  in  the  vertical  members  usually 
are  not  critical  owing  to  the  factor  of  safety.  However,  the 
main  problem  is  not  one  of  stresses  but  one  of  dynamic  in¬ 
stability.  The  general  notion  that  the  vertical  acceleration 
is  considerably  less  than  the  horizontal  acceleration  is  also 
false.  Obviously  the  average  intensity  of  the  vertical  and 
horizontal  components  of  an  earthquake  depends  upon  the  re¬ 
lationship  between  the  epicentral  distance  to  the  seismograph 
and  the  focal  depth.  If  these  two  distances  were  equal,  the 
average  intensities  of  the  vertical  and  horizontal  components 
would  also  be  equal.  Therefore  the  problem  is  magnified 
especially  for  deep  focus  earthquakes.  However,  this  rela¬ 
tionship  can  only  be  assumed  for  the  average  intensities . 

For  the  extreme  values  of  acceleration  such  simple  relations 
do  not  hold,  since  the  vertical  and  horizontal  accelerations 
are  actually  components  of  a  random  process.  From  the  data 
that  he  takes  from  Housner  (73) ,  Despeyroux  points  to  this 
fact  in  the  following  way  (55) ; 

"At  a  distance  of  48  km.  from  the  epicentre  of  the 
El  Centro  earthquake  of  1940,  the  focal  depth  being 
24  km.,  the  maximum  vertical  acceleration  was  60% 
of  the  maximum  horizontal  one;  and  for  the  Taft 
earthquake  of  1952,  it  was  also  60%  at  a  distance 
of  64  km.  for  the  same  focal  depth  of  24  km." 

Therefore,  if  the  failure  criterion  is  taken  to  be  the  maxi¬ 
mum  value  and  not  the  cumulative  effect,  then  even  for  shal¬ 
low  earthquakes  and  considerable  epicentral  distances  the 
vertical  acceleration  cannot  be  neglected  in  comparison  to 
the  horizontal  component.  In  fact,  for  a  stability  analysis 
the  maximum  value  would  be  the  failure  criterion.  As  for 
the  lack  of  evidence  that  failure  does  not  occur  due  to 
vertical  acceleration,  the  following  argument  can  be  pre¬ 
sented.  First,  the  eimount  of  data  that  we  have  adsout  the 
failures  during  earthquakes  is  not  enough  to  exclude  any 
mode  of  failure.  Second,  the  evidence  of  the  Good  Friday 
earthquake  of  Alaska  as  well  as  that  of  Venezuela  are  suf¬ 
ficient  to  conclude  exactly  the  opposite  (74)  .  However, 
even  if  there  were  no  evidence  whatsoever  about  stability 
failure  due  to  vertical  earthquake  acceleration,  this  prob¬ 
lem  still  cannot  be  discarded.  The  problem  is  a  coupled  one 
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between  flexural  lateral  deflection  and  axial  vibration. 
Therefore  the  transfer  of  energy  from  one  mode  to  the  other 
is  quite  probeible.  In  fact,  it  has  been  shown  that  the 
energy  of  one  mode  can  almost  be  totally  transferred  into 
the  other  one  (75) .  Therefore,  a  seemingly  lateral  failure 
can  actually  be  induced  by  vertical  forces .  This  is  nothing 
new,  and  it  is  a  natural  outcome  of  any  coupled  system. 

As  was  mentioned  earlier,  most  codes  do  not  include  any 
provisions  for  vertical  forces.  Two  exceptions  are  the 
Romanian  and  the  French  Aseismic  Codes.  The  recommendations 
of  the  Romanian  regulations  for  aseismic  design  read  as  fol¬ 
lows  (62):  "Besides  horizontal  forces,  supplementary  verti¬ 
cal  forces  are  also  taken  into  account  by  increasing  the 
vertical  forces  from  25  to  100  per  cent  according  to  the 
degree  of  intensity."  Also  the  recommendations  of  the 
French  Aseismic  Code  are  given  as  (55) : 

"Vertical  and  horizontal  seismic  forces  are  con¬ 
sidered  to  act  simultaneously.  As  the  maximum 
horizontal  and  vertical  responses  are  not  likely 
to  occur  at  tlie  same  time,  the  vertical  component 
tciken  into  account  is  only  a  part  of  the  maximum 
response.  However,  for  structural  elements  in 
which  lateral  forces  are  of  no  concern,  like  vault- 
ties  and  cantilevers ,  higher  values  might  be  im¬ 
posed,  with  the  idea  of  penalizing  projected 
facades  and  so  on." 
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PART  III:  POSSIBLE  IMPROVEMENTS  IN  CODE  FORMAT 


General 

The  Uniform  Building  Code  (8)  specifies  equivalent 
static  loads  for  design.  Although  these  static  loads  are 
based  on  dynamical  considerations,  discrepancy  in  their  ac¬ 
curacy  is  noted  in  many  areas.  For  example,  Kanai  and 
Yoshizawa  (76)  have  shown  that  the  equation  given  in  the 
building  code  for  the  natural  period  of  vibration  of  build¬ 
ings  is  not  in  agreement  with  experimental  results .  Simi¬ 
lar  conclusions  have  been  drawn  by  Arias  and  Husid  (77) . 
Blume  (12)  has  shown  that  the  cumulative  displacement  due  to 
dynamic  loads  could  critically  exceed  static  load  deflec¬ 
tions  specified  by  codes.  A  review  of  existing  codes  as 
given  in  Part  II  indicates  that  the  Uniform  Building  Code 
takes  into  account  the  dynamic  effects  in  a  different  maui- 
ner  from  French  and  other  codes.  At  the  present  time,  there 
is  a  considercible  amount  of  research  work  being  done  in 
earthquake  engineering  (78) .  Nevertheless,  the  current  re¬ 
search  programs  are  primarily  concerned  with  specific  prob¬ 
lems.  It  seems  to  be  desirable  to  review  the  philosophies 
behind  these  design  codes. 

In  Part  III,  some  design  philosophies  including  the 
concept  of  structural  control  are  examined  in  the  context 
of  earthquake  engineering.  Design  methodologies  are  also 
reviewed  with  emphasis  on  the  design  tree  and  the  fully 
systematic  method,  which  are  later  applied  to  the  formula¬ 
tion  of  a  new  code  format.  Moreover,  the  expected  problems 
in  the  implementation  of  the  proposed  code  format  are  listed 
along  with  an  outline  of  necessary  procedures  for  their  so¬ 
lutions  . 


Some  Design  Philosophies  and  Methodologies 

According  to  Asimow  (79) ,  there  are  the  following  three 
parts  in  a  philosophy  of  engineering  design: 

1.  a  set  of  consistent  principles  and  their  logical 
derivations , 

2.  an  operational  discipline  leading  to  action,  and 

3 .  a  feedback  mechanism  measuring  the  advantages , 
detecting  weakness,  and  pointing  out  the  di¬ 
rection  for  improvement. 
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The  principles  and  their  logical  derivations  being  used 
in  support  of  existing  design  codes  are  primarily  determin¬ 
istic  in  nature.  Whenever  it  is  possible,  dynamic  effects 
are  represented  by  "equivalent  static  forces."  Moreover, 
the  seismic  structures  are  designed  to  remain  elastic  or 
nearly  so  under  the  influence  of  moderate  and  frequent 
earthquakes,  but  to  allow  local  yielding  without  collapsing 
under  severe  and  infrequent  earthquakes  (50,  80) .  However, 
the  effect  of  repeated  plastic  deformation  is  known  to  be 
cumulative  (81,  82) .  Recently,  it  was  found  to  be  possible 
to  have  low-cycle  fatigue  failure  of  structures  under  earth¬ 
quake  loads  (18)  . 

The  relations  between  pure  research,  applied  research, 
and  practical  applications  are  given  by  the  Committee  on 
Earthquake  Engineering  Research  as  shown  in  Figure  1  (2) . 
Because  of  differences  in  objectives,  the  pure  research  in 
seismology  does  not  contribute  directly  to  engineering  de¬ 
sign  under  earthquake  conditions.  As  an  example,  the 
structural  engineer  is  interested  to  know  the  structural  de¬ 
formations  caused  by  strong-motion  earthquakes,  which  can¬ 
not  be  measured  with  overly  sensitive  seismographs  utilized 
by  seismologists.  In  general,  however,  research  results  in 
various  related  specialties  such  as  seismology,  geology,  and 
geophysics  are  helpful  to  earthquake  engineer!:  g  ii  terms  of 
basic  information  such  as  (a)  the  locations,  ma-:;nitudes ,  and 
frequencies  of  occurrence  of  earthquakes,  (b)  the  distinc¬ 
tion  between  active  fault  zones  and  potentially  active 
faults,  (c)  the  locations  and  amounts  of  surface  fault  move¬ 
ments,  and  (d)  the  nature  of  crustal  deformation  (2) .  The 
applied  research  as  shown  in  Figure  1  may  be  described  as  a 
feedback  mechanism,  which  serves  to  point  out  the  direction 
for  improvement  of  the  practical  application.  A  design  code 
may  be  considered  as  a  component  of  practical  applications 
as  shown  in  Figure  1. 

Several  attempts  have  been  made  to  "control"  the  struc¬ 
tural  response  recently.  These  include  (a)  the  application 
of  "static  control"  by  Zuk  (83) ,  (b)  "active  systems"  by 
Wright  (84)  and  Nordell  (85) ,  and  (c)  "adaptive  systems"  by 
Yao  (86)  .  As  an  exeimple  of  the  adaptive  system,  the  bi¬ 
linear  system  was  studied  is  some  detail  (34,  87,  88). 

Along  the  same  line  of  thought,  a  statistical  study  of  pas¬ 
sive  motion-reducing  devices  in  earthquake  engineering  was 
completed  very  recently  (89)  .  While  these  studies  indicate 
possible  improvements  in  the  seismic  structural  design,  the 
available  results  are  not  yet  directly  applicable. 

Generally  speaking,  specifications  refer  to  (a)  func¬ 
tional  requirements  and  (b)  physical  limitations  in  the  en¬ 
gineering  design  (90)  .  This  information  can  be  provided  by 
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clients,  authorities  or  agencies  concerned  with  public 
safety,  and  engineering  organizations.  In  earthquake  en¬ 
gineering,  design  specifications  have  been  proposed  by  pro¬ 
fessional  engineering  organizations  such  as  the  Structural 
Engineers  Association  of  California  (47)  and  adopted  as  a 
part  of  the  Uniform  Building  Code  (8) .  The  guidance  for 
the  design  of  seismic  structures  is  also  given  in  the  form 
of  manuals  as  the  one  prescribed  by  the  Departments  of  the 
Arity,  the  Navy  and  the  Air  Force  in  1966  (10)  . 

A  comprehensive  summary  of  design  methodologies  is 
given  by  Eder  (91).  These  include  (a)  expenses,  (b)  modi¬ 
fication  and  running  redesign,  (c)  checklists,  (d)  design 
trees,  (e)  the  fully  systematic  method,  and  (f)  the  system 
search  methods.  In  this  section,  the  design  trees  and  the 
fully  systematic  methods  are  examined  in  some  detail,  be¬ 
cause  they  are  used  in  the  development  and  presentation  of 
the  new  code  format. 

In  engineering  design  (92) ,  the  cyclic  design  sequence 
of  (a)  analysis  of  the  problem,  (b)  proposing  solutions,  (c) 
delineating  these  solutions,  and  (d)  modifying  them  can  be 
plotted  as  a  design  tree  as  shown  in  Figure  2,  where  a 
problem  is  denoted  by  a  vertical  line  and  a  solution  is  de¬ 
noted  by  a  slanting  line.  The  tree  is  started  with  a  state¬ 
ment  of  the  primary  problem,  which  is  referred  to  as  the 
zeroth  order,  Qo«  The  alternative  solutions,  first  order, 
are  represented  by  slanting  lines  Ai,  A2/  etc.  For  each 
alternative  solution,  such  as  A2f  there  may  be  several  prob¬ 
lems  to  be  solved,  i.e.,  Q21»  022/  etc.  Any  one  or  more 
of  the  alternative  solutions  are  sufficient  to  solve  the 
problem.  On  the  other  hand,  it  is  necessary  to  solve  all 
problems  before  a  lower  order  alternative  solution  is  ob¬ 
tained.  As  an  example,  the  bold  lines  in  Figure  2  indicate 
one  possible  set  of  solutions. 

In  the  fully  systematic  method,  a  list  of  factors  rele¬ 
vant  to  the  problem  or  its  solution  factors  must  be  decided 
in  the  beginning.  It  is  important  to  include  all  ideas 
without  criticism  and  to  allow  duplications.  Each  factor, 
described  by  a  single  statement,  should  be  numbered  in  the 
order  of  appearance.  Furthermore,  any  factor  that  contains 
a  possible  solution  of  a  sub-problem  should  begin  a  new 
category.  Eventually,  all  factors  must  be  classified  into 
a  set  of  categories.  The  influence  of  one  category  upon 
another  can  be  studied  with  the  use  of  an  interaction  matrix, 
which  can  be  transformed  into  interaction  nets.  The  result¬ 
ing  interactions  can  now  be  transformed  into  performance 
specifications,  which  are  used  to  define  performance  (e.g., 

"a  column  must  support  the  floor  load",  not  "a  column  must  be 
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of  a  certain  size").  The  objective  of  the  performance  speci¬ 
fications  is  to  completely  separate  the  problem.  Therefore, 
these  specifications  must  not  make  reference  to  shape,  mate¬ 
rial,  or  design.  Each  performance  specification  can  now  be 
considered  to  conceive  the  idea  of  solutions  and  to  list 
them  in  a  morphological  chart  (93)  .  The  method  of  the  de¬ 
sign  tree  is  used  in  the  process  of  obtaining  an  improved 
seismic  design  code  in  this  report. 


Direct  Approach 

The  "primary  question,"  Qq,  in  the  present  research  pro¬ 
gram  is  "what  constitutes  a  seismic  design  code?"  There  are 
two  first-order  alternative  answers,  namely  A^:  the  "equiva¬ 
lent  static  forces"  approach  which  summarizes  all  of  the 
major  existing  codes  as  shown  in  Figure  3,  and  A2:  a  direct 
approach  which  is  shown  in  Figure  4.  All  detailed  "answers," 
A's,  in  both  Figures  3  and  4,  are  tabulated  in  Table  3  im¬ 
mediately  following  these  two  figures.  The  equivalent  static 
forces  approach  was  discussed  in  detail  in  Part  II.  There¬ 
fore,  only  the  "direct  approach"  is  elaborated  and  discus¬ 
sed  herein. 

To  date,  the  existing  seismic  design  codes  specify  a 
set  of  equivalent  static  forces  for  use  by  structural  en¬ 
gineers  in  their  designs.  This  approach  has  the  advantage 
of  being  easy  to  use,  because  the  engineer  is  spared  the 
need  to  understand  structural  dynamics.  On  the  other  hand, 
the  simple  formulae  given  in  the  code  cannot  possibly  cover 
all  possible  situations  accurately  and  precisely.  With  (a) 
more  sophisticated  education  of  young  engineers,  and  (b) 
widespread  use  of  tools  such  as  the  electronic  digital  com¬ 
puter  in  recent  years,  it  is  believed  that  the  present  gene¬ 
ration  of  engineers  can  be  trusted  with  more  analysis  on 
their  own  than  ever  before.  The  "direct  approach"  as  sug¬ 
gested  is  also  a  more  flexible  one,  which  is  expected  to 
yield  safer  as  well  as  more  economical  structural  designs. 

As  shown  in  Figure  4,  there  can  be  three  second-order 
questions  in  the  direct  approach  A2,  namely,  Q21:  excita¬ 
tion,  Q2Z‘  structural  model  and/or  analysis,  and  Q23:  re¬ 
sponse  limitations  for  functioning  and  safety.  There  can 
be  two  alternative  third-order  answers  to  each  of  these 
second-order  questions.  For  Q2i»  the  excitation  can  be 
specified  by  A211:  a  single  sample  function  (deterministic) , 
or  by  A212*  ^  family  of  sample  functions  (probabilistic) . 

For  Q22»  the  code  can  give  A221J  the  methods  of  computing 
parameters  for  the  mathematical  model  of  the  structure,  or 
A222*  the  model  as  well  as  computed  responses  using  dis¬ 
crete  mechanics.  The  alternate  approach  A222  can  be 
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Figure  4.  "Design  Tree"  for  a  Direct  Approach 
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Table  3  (continued) 


^11231* 


^1X241 ‘ 

^11311* 

^11321* 

^11331* 

^11341’ 

^11411' 


.065  <  B 


.13;  Large  Span  &  Little 
Partitions  (medium 
damping) 


.095 

.06  6  =  -  <  .20;  Skeleton,  no  external 

-  friction  (low  damping) 


q(h)  =  shape  £cn.  of  the  fundamental  mode, 

r (T)  =  max.  response  (acceleration  of  a  SDF 

system  with  the  same  period  T  and  damp¬ 
ing  as  the  fundamental  mode* 


♦The  first  three  modes  are  considered  for 
slender  structures  with  low  deunping. 

f  =  0.8  for  deep  foundation  of  firm  rock 
=  1.3  for  piling  in  moist  ground 

z  =  0.8  or  0.9  depending  on  seismic  zone 


s  =  0.6  for  steel  freune  or  wooden  structure 
on  firm  ground, 

=  0.9  for  semi-firm  ground  and  reinforced 
concrete  structure, 

=  1.0  for  soft  ground 

c  -  0.2  for  structures  up  to  16  m.,  and  add 
0.01  for  each  4  m.  in  addition  to  16  m. 


W  =  total  weight 

Type  1  =  structures  having,  at  right  angles 

with  the  direction  being  considered, 
two  or  more  elements  capable  of  re¬ 
sisting  shear,  and  whose  deformation 
under  lateral  loads  is  essentially 
due  to  flexure  in  members. 
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^11421* 

^11431* 

^*11511’ 

^11521' 


^11531' 


TaUble  3  (continued) 


Type  2  s  Structures  whose  deformations  under 
lateral  loads  are  due  to  shearing 
stresses  or  axial  forces  in  members . 

Type  3  =  Elevated  tanks,  chimney  stacks,^ etc. 


k  *  1.3  for  Group  A,  government,  municipal 
and  public  buildings 

=  1.0  for  Group  B,  private  housing 

I 

0  for  Group  C,  isolated  and  .unimportant 
buildings 

W  *  total  weight 

Z  =  D.L.  +  0.8  L^L.  (for  warehouses,  use  full 
^  L.L.) 

GEOFIAN  Scale _ 7 _ ^ 

a  1/40  1/20  1/iO 

4  Types  of  Bldgs:  ' 

I;  Massive  bldgs,  of  great  importance,  in¬ 
crease  rating  by' 1 

II:  Public  Bldgs,  industrial  bldgs,  of 
primary  importance,  same  rating 

III:  Industrial  bldgs,  of  secondary  importance, 
one-story  houses,  same  ratings  for  "6", 
"7",  reduce  "8"  and  "9"  by  1. 

IV:  Use  rating  "6"  for  all  farm  bldgs., 
barns ,  sheds ,  etc . 

0.6  <  6  =  <  3.0,  emd  use  1.6  6  for  slender 

and  flexible  structures 
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Table  3  (continued) 


^11541* 


A 

^121* 

A 

^122* 

A 

^131* 


N 


N 


I  Q.q  (X.) 


i=l 


use  1  mode  for  normal  buildings,  and  use 
the  first  3  modes  for  slender  and  flexible 
bldgs . 


=  S^d.S  e^  +  0.5  B^) 

shear  eccentricity  max.  dimension 
e_  =  1.5  e^  ±  0.05  D 

D  C 

May  be  increased  by  33% 


Aj^22‘  increased  by  50% 

May  be  increased  by  50%  for  reinforcing  steel,  and 
100%  for  ocncrete  in  reinforoed  oczicrete 

^134*  increased  by  50%  for  steel,  and  33% 

for  concrete 

^135*  increased  by  40%  for  steel,  20%  for 

^  reinforced  concrete,  and  0%  for  prestressed 

concrete 

^136’  increased  by  50%  for  good  soil  con¬ 

ditions,  30%  for  fair  soil  conditions,  and 
0%  for  bad  soil  conditions 


A,^,:  French  and  Romanian  codes  mentioned  it.  How¬ 

ever,  it  is  not  clear  as  to  how  to  do  it  in 
designs 

^71111*  Take  a  typical  set  of  earthquake  records, 

e.g.,  1940  El  Centro 

A-...,:  Find  the  least  favorable  excitation,  e.g., 

Drenick 


^21117111*  Zoning  (deterministic) 
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Table  3  (continued) 


^21112121* 

Soil  Conditions  (deterministic) 

^21211* 

Take  a  set  of  artificial  functions,  e.g., 
Shinozuka  and  Sato,  Amin  eu:id  Ang,  etc. 

^21212- 

Take  several  sets  of  artificially  generated 
functions,  e.g.,  Jennings,  Housner  and  Tsai 

^2122111* 
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^22111' 
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Analytical 

^22122* 

Experimental 
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♦discrete  mechanics,  based  on  the  calculus  of 
finite  differences,  can  be  used  to  obtain  field 
or  functional  solutions  for  systems,  which  are 
more  accurately  represented  by  a  lattice  or  a 
network  of  elements. 

^2212'  ^22222'  ^22232'  ^22242*  Nonlinear  analysis 
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Table  3  (continued) 


^2311231 ‘ 
^2312111 • 
^312121' 
^2321111* 
*“2321121  • 
*2321211* 
*2321221* 
*2321222* 
*2321231’ 
*2322111’ 
*2322121’ 


Stable 

Haunmering  Prevented 
Comfortable 

i  V  -  Pi 

P (Unstable)  ^  P2 

P<lwl  iyilmit*  iP3 

P(I  >  1)  < 

i  ®cr.'  i  Ps 
P (Unstable)  ^  Pg 

P  (Hammering)  p^ 

P (Discomfort)  ^  Pg 

Pi's  should  be  functions  of  type  of  structure, 
function  (importance)  of  structure,  consequence 
of  failure,  material  properties,  etc. 


considered  as  an  intermediate  code  between  the  existing  ones 
cuid  that  given  by  A221«  For  Q23r  “tho  response  limitations 
for  functioning  and  safety  can  be  specified  by  A231:  a  sin¬ 
gle  value  (deterministic)  or  by  a  probability  (probabilis¬ 
tic)  .  The  fourth-order  questions  related  to  these  third- 
order  approaches  are  discussed  in  the  next  section,  which 
also  includes  fifth-order  cuiswers. 
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Implementation  Problems  and  Possible  Approaches 


1.  Excitation  (Q21) 


To  specify  the  excitation  for  the  dynamic  analysis  of 
structures,  time  record(s)  of  the  gound  motion  is.^are) 
needed.  One  possible  answer,  which  has  been  used  by  many 
consulting  engineers,  is  use  a  typical  set  of 

strong-motion  earthquake  records  such  as  those  of  the  famous 
1940  El  Centro  earthquake.  An  alternate  approach,  A2iij2» 
is  to  make  use  of  a  "least  favorable  excitation."  Drenick 
(94)  showed  that,  for  a  class  of  excitation  x(t)  with  a 
finite  total  energy  (a)  x(t)  =  (M/N)h(-t)  causes  the 

least  favorable  response  of  a  linear  structure  with  the  im¬ 
pulse  response  function  h(t),  and  (b)  the  maximum  energy 
response  is  MN,  where 

=  /  "h^(t)dt 
*00 

More  recently,  Shinozuka  (94)  showed  how  to  estimate  the 
least  favorable  structural  response  to  probabilistic  and 
earthquake-like  excitation.  With  digital  or  analog  com¬ 
puters,  it  is  also  possible  to  apply  Monte  Carlo  methods 
using  artificially  generated  earthquake- like  motions  (27- 
34) .  All  these  time  functions  should  be  modified  acording 
to  the  location  of  the  structure,  which  in  turn  can  be  ex¬ 
pressed  as  a  function  of  the  zoning  and  local  soil  condi¬ 
tions  as  shown  in  Figure  4 . 


2.  Structural  Model  and/or  Analysis  (022^ 


To  analyse  a  structural  system  subjected  to  dynamic 
loads,  it  is  necessary  to  have  a  mathematical  model.  With 
the  recent  advances  in  digital  computation,  the  discrete 
model  is  a  possible  and  useful  one  (A221  Figure  4)  .  Fre¬ 
quently,  the  structure  can  be  represented  by  a  set  of  ordi¬ 
nary  differential  equations  with  parameters  representing 
the  inertia  (translational  and/or  rotational) ,  the  damping 
(linear  or  nonlinear) ,  and  stiffness  (linear  or  nonlinear) . 
Although  the  study  of  structural  dynamics  has  reached  a 
mature  stage  with  certain  given  models  such  as  the  linear 
shear-type  building  frames,  theoretical  and  experimental 
work  seems  to  be  needed  in  relating  these  parauneters  to  real 
structures.  For  example,  though  many  studies  have  been  made 
on  material  and  structural  damping  (96,  97),  research  and 
development  work  is  needed  to  specify  damping  coefficients 
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for  various  types  of  structures  and  partitions.  Moreover, 
the  effect  of  cracking  or  cumulative  damage  due  to  previous 
earthgucdces  on  these  parameters  is  not  clearly  understood 
at  the  present  time. 

An  intermediate  seismic  design  code  might  be  the  speci¬ 
fication  of  response  spectra  k^222  Figure  4) . 

The  dynamic  analysis  of  multi-story  building  structures 
subjected  to  earthquake  excitations  is  a  complicated  problem 
ir  the  field  of  mechanics.  Redundancy  of  such  structures, 
random  nature  of  loads,  difficulties  in  determining  contri¬ 
bution  to  damping  by  architectural  clothing,  behavior  of 
foundation  under  shakedown,  questions  related  to  dynamic 
stability  of  structure,  significance  of  nonlinear  analysis 
accounting  for  secondary  stresses,  and  complexity  of  elastic- 
plastic  behavior  of  frames  are  among  problems  which  contrib¬ 
ute  to  the  intricacy  of  the  subject  area. 

A  reasonable  objective  for  the  codes  and  specifica¬ 
tions  in  line  with  economic  factors  appears  to  be  that 
seismic  design  codes  should  provide  guidelines  for  design  of 
structures  : 

(a)  to  withstand  frequent  earthquakes  without  signif¬ 
icant  damage  to  structure  or  facings, 

(b)  to  avoid  total  collapse  i’'  major  earthquedces  , 

(c)  to  have  sufficient  stiffness  to  exclude  sensible 
vibration  when  subjected  to  wind  load  in  normal 
conditions . 

Subject  to  a  few  restrictions,  it  is  feasible  to  pre¬ 
dict  the  dyneunic  behavior  of  almost  any  structure  by  the  use 
of  high  speed  digital  computers,  as  well  as  mathematical 
manipulation.  The  alternative  approach  A222  <^®~ 

velopment  of  working  formulas  which  are  simple  enough  to  be 
included  in  the  codes  for  use  by  practicing  engineers  and 
yet  describe  the  dynamic  behavior  of  a  broader  category  of 
structures.  A  more  detailed  outline  of  the  problem  areas 
will  include  the  following: 

(a)  Formulas  or  charts  for  response  spectrum  of  rigid 
emd  braced  multi  degree-of-freedom  freunes  subject 
to  decided  critical  earthquake  excitation  to  in¬ 
clude  (i)  displacement  response  spectrum  in  shear, 
torsi jn,  or  bending,  (ii)  shear  response  spectrum, 
(iii)  bending  response  spectrum,  (iv)  torsion  re¬ 
sponse  spectrum,  (v)  determination  of  influence 
on  response  spectrum  (shearing  forces,  bending, 
torsion,  and  displacements)  accounting  for  the 
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interaction  of  vertical  earthquake  excitation  with 
the  horizontal  accelerations,  (vi)  determination  of 
influence  on  response  spectra  (displacements, 
shearing  forces,  bending,  and  torsion)  in  a  non¬ 
linear  cuialysis  accounting  for  secondary  stresses 
and  the  interaction  between  lateral  and  vertical 
forces  known  as  the  P-A  effect,  (vii)  dynamic  re¬ 
sponses  of  elastic-plastic  frames,  (viii)  deter- 
minaticn  of  time -dependent  excitation  most  cri¬ 
tical  for  any  particular  structure,  and  (ix)  formu¬ 
lation  of  dynamic  stability  of  structures  under 
shakedown . 

(b)  Energy  absorption  capacity  of  structures: 

(i)  formulas  for  natural  frequencies  of  vibration 
of  rigid  and  braced  structures .  Structures  with 
orderly  pattern  and  boundary  conditions  as  well  as 
stepped-up  structures  and  structures  with  append¬ 
ages  need  to  be  considered.  (ii)  Influence 
charts  for  the  effects  of  variations  in  mass  dis¬ 
tribution,  ratio  of  stiffnesses,  shapes,  damping 
factors,  sizes,  etc.,  on  energy  absorption  capa¬ 
city  of  structures.  (iii)  Optimum  distribution 
of  ductility  to  aid  maximum  energy  absorption 
capacity . 

(c)  Problems  related  to  foundations  and  structural 
materials  and  their  properties' 

(i)  necessary  modifications  in  response  spectra 
due  to  nature  of  subfoundation  soil  layers,  (ii) 
liquefaction  of  sand  and  silt  under  intense  verti¬ 
cal  acceleration,  (iii)  the  uneven  foundation  pro¬ 
perties  and  their  effects,  (iv)  damping  factors 
most  closely  descriptive  of  different  construction, 
(v)  dynamic  behavior  of  reinforced  concrete  mem¬ 
bers,  joints,  precast  or  prestressed  members,  (vi) 
dynamic  characteristics  of  multistory  reinforced 
concrete  buildings. 

Besides  those  approximate  methods  which  fall  under  the 
classification  of  static  equivalent  methods,  two  methods  of 
solution  are  available.  Clough  (98)  suggested  matrix  formu¬ 
lation  of  each  independent  structure  and  use  of  generally 
applicable  but  somewhat  lengthy  computer  programs .  Newmark 
(80)  formulated  the  response  spectra  for  single  degree-of- 
freedom  system.  Then,  a  multi  degree-of -freedom  system  is 
treated  by  trcuisfoimiation  to  an  equivalent  single  degree-of- 
freedom  system.  At  the  present  time,  this  technique  is 
being  pursued  at  the  University  of  Illinois  (80) . 
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Research  results  under  the  present  project  have  indi¬ 
cated  the  suitcdsility  of  employing  an  alternate  approach  to 
formulate  the  response  spectra  of  multi  degree-of-freedom 
systems.  The  proposed  method  of  solution  bypasses  much  of 
the  difficulty  inherent  in  the  earlier  lines  of  thinking  by 
assuming  a  discrete  model  representation  of  highly  redundant 
structures  and  using  the  calculus  of  finite  differences  and 
the  concepts  of  discrete  mech2aiics  to  derive  a  solution. 

These  concepts  are  illustrated  in  Appendix  A,  where  a  single 
bay  multi-story  frame  is  analysed.  Closed  functional  solu¬ 
tions  are  found  for  natural  frequencies  of  vibration  and  for 
lateral  force  distribution.  Furthermore,  the  formulation  of 
the  dynamic  response  in  combined  torsion  and  biaxial  shear 
is  given. 

Earthqucikes  produce  torsional  oscillation  and  stress  in 
structures  due  to  nonsymmetric  mass  and  stiffness  prop¬ 
erties  of  buildings.  There  is,  however,  a  lack  of  rational 
solutions  for  dyneunic  behavior  of  structures  in  torsion  due 
to  complexity  of  the  problem.  The  problem  of  dynamic  be¬ 
havior  in  torsion  and  biaxial  shear  has  received  little  at¬ 
tention  even  though  the  interaction  of  bieixial  shear  dis¬ 
placement  with  torsion  is  apparent.  A  general  formulation 
of  the  dyneunic  response  of  structures  in  torsion  and  biaxial  shear 
is  given  in  Appendix  A.  This  formulation  as  well  as  other 
zk  search  work  based  on  discrete  mechanics  performed  under 
this  research  program  serve  to  illustrate  the  potential  of 
the  techniques  of  discrete  mechanics. 

3.  Response  Limitations  (Q23) 


For  the  consideration  of  functioning  and  safety  of 
structures,  limitations  can  be  specified  in  terms  of  either 
maximum  response  (A231)  reliability  or  probability  of 
failure  (A232) •  These  specifications  can  be  made  with  the 
use  of  statistical  methods  and  on  the  basis  of  seismic  be¬ 
havior  of  structures .  A  review  of  statistical  methods  in 
earthquake  engineering  is  given  in  Appendix  B.  Moreover, 
the  seismic  behavior  of  concrete  and  metal  structures  is 
sximmarized  and  discussed  in  Appendix  C. 

When  the  structure  behaves  in  a  linear  manner,  the  de¬ 
formation  is  proportional  to  either  stress  or  strain.  There¬ 
fore,  strength  limits  can  be  specified.  On  the  other  hand, 
deformation  limits  must  be  given  in  the  case  of  inelastic 
structural  behavior.  Because  the  damage  due  to  plastic  de¬ 
formation  is  cumulative,  it  is  also  desirable  to  restrict 
the  allowable  damage  in  building  structures. 
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The  st.ruct.ural  stability  under  the  vertical  as  well  as 
horizontal  ground  motions  is  a  complicated  problem.  There 
appears  to  be  little  work  done  in  this  regard  at  the  present 
time.  The  hammering  effect  of  closely-spaced  buildings 
against  one  another  should  be  prevented.  Moreover,  it  is 
desireible  to  minimize  the  discomfort  of  human  inheibitants 
during  earthquakes . 

In  both  the  deterministic  specification  A231  and  the 
probedjilistic  approach  ^222/  desirable  to  medce  sta¬ 

tistical  analysis  of  available  data.  In  addition,  the  de¬ 
cision  analysis  can  be  a  useful  tool  in  the  choice  of  these 
limitations  (99) . 

In  Appendix  B,  following  an  introduction  to  basic  defi¬ 
nitions,  topics  covered  inc.^ude  statistics  of  strong-motion 
earthquakes,  random  response  of  seismic  structures,  struc¬ 
tural  safety  in  earthquake  engineering,  and  decision  anal¬ 
ysis.  In  Appendix  C,  a  review  is  made  of  the  seismic  be¬ 
havior  of  various  types  of  structures  during  recent  earth¬ 
quakes.  In  Appendix  D,  a  simple  design  problem  is  givai  as  an 
illustration  of  the  "direct  ^aproach"  portion  of  the  design  tree. 
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PART  IV:  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


A  survey  of  exlst.ing  and  availcd>le  seismic  design  codes 
was  made.  The  "design  tree"  technique  was  used  to  summarize 
these  existing  specifications  as  well  as  to  present  possible 
improvements .  Moreover ,  some  original  work  has  been  per¬ 
formed  applying  the  discrete  mechanics  to  structural  dynamics 
as  given  in  Appendix  A.  Literature  reviews  of  statistical 
methods  in  earthquake  engineering  and  seismic  behavior  of 
concrete  and  metal  structures  are  given  in  Appendices  B  and 
C,  respectively. 

A  "direct  approach"  to  the  seismic  design  code  is  pro¬ 
posed  herein.  With  this  direct  approach,  the  code  would 
specify  (a)  the  excitation,  (b)  the  structural  model,  and 
(c)  the  response  limitations.  It  is  believed  that  engineers 
using  this  type  of  code  will  be  able  to  design  building 
structures  in  a  more  precise,  realistic,  and  economic  manner. 
However,  it  might  be  desirable  to  perform  this  type  of  anal¬ 
ysis  and  design  enough  times  to  determine  if  siitpler 
analysis  would  provide  the  required  information  at  cost  and  time  savings. 


Recommendations 


The  following  is  a  list  of  problems,  which  are  delin¬ 
eated  for  possible  further  studies  by  CERL  personnel  or 
interested  persons  elsewhere. 

1.  Specification  of  Ground  Motions  -  the  concept  of  the 
least  favorad>le  excitation  originally  developed  by 
Drenick  (94)  and  recently  modified  by  Shinozuka  (95) 
seems  to  be  very  promising.  Nevertheless,  further  work 
is  required  before  the  present  highly  mathematical  forms 
can  be  incorporated  into  any  design  code.  In  case  it 

is  desirable  to  specify  several  artificially  generated 
records  for  statistical  analysis,  development  work  is 
also  required  for  code  adoption  of  existing  techniques. 
Moreover,  studies  should  be  made  to  develop  methods 
with  which  these  design  ground-motions  can  be  modified 
with  respect  to  the  location  and  soil  conditions  at 
the  site  of  the  structure. 

2.  Computation  of  Pareimeters  for  Given  Structures  -  Mathe¬ 
matical  models  should  be  as  representative  of  the 
actual  structure  as  possible.  In  the  iterative  design 


process,  the  engineer  should  be  given  a  set  of  fomulae 
for  the  computation  of  structural  parameters  such  as 
damping  coefficients.  Personnel  of  a  research  facility 
such  as  CERL  are  in  a  iinique  position  of  conducting 
analytical  as  well  as  experimental  investigations,  which 
could  result  in  successfully  relating  abstract  and 
mathematical  models  to  real  structures. 

Response  Spectra  -  As  an  intermediate  step,  formulae 
and/or  charts  can  be  developed  for  response  spectra  of 
multistory  frame  structures  subjected  to  various  earth- 
gucJce  loads .  These  solutions  can  be  obtained  as  func¬ 
tions  of  number  of  floors,  ratio  of  stiffnesses,  dif¬ 
ferences  in  natural  frequencies  of  vibration,  mass,  etc. 
Criteria  for  the  intensity  and  distribution  of  loads  to 
be  suggested  in  the  codes  must  be  provided  based  on  the 
formulation  of  the  response  spectra.  Displacement  re¬ 
sponse  spectrum  in  shear,  torsion,  bending  and  their 
combination;  as  well  as  shearing  force  response  spec¬ 
trum,  bending  response  spectrum,  torsion  response  spec¬ 
trum  need  to  be  formulated  and  the  recommended  load 
distribution  for  each  action  needs  to  be  defined.  The  ex- 
eunples  as  given  in  Appendix  A  (some  are  original  con¬ 
tributions  resulting  from  this  research  program)  il¬ 
lustrate  the  potential  usefulness  of  discrete  mechanics 
in  structural  dynamics. 

Dynamic  Stability  of  Seismic  Structures  -  Although 
methods  for  three-dimensional  dynamic  analysis  of 
seismic  structures  are  available,  the  stability  of 
seismic  structures  has  not  been  studied  extensively. 

To  minimize  the  chance  of  structural  failure  due  to  in- 
stcdjility  under  earthquake  loads,  it  is  necessary  to 
develop  new  knowledge  in  this  regard. 

Inelastic  Behavior  of  Structures  -  Because  the  princi¬ 
ple  of  superposition  is  not  applicable  to  inelastic 
structural  analysis,  most  available  solutions  are  re¬ 
sults  of  numerical  analyses .  Perhaps ,  improvements  car. 
be  made  in  these  numerical  methods  as  well  as  problem- 
oriented  computer  languages . 

Low-Cycle  Failure  of  Seismic  Structures  -  Experimental 
results  are  availcd>le  for  concrete  members  (100)  as 
well  as  steel  members  (101)  subjected  to  repeated  ap¬ 
plications  of  prescribed  large  deformations .  The  fati¬ 
gue  behavior  of  structures  under  random  loads  needs 
further  study  (102) .  Moreover,  there  seems  tr  be  a 
need  for  experimental  studies  on  the  fatigue  behavior 
of  structures  under  earthquake-like  random  loads. 


though  some  data  are  available  in  aerospace  engineering 
fields  for  coupon  specimens  (103) . 

7.  Reliability  Considerations  -  Structural  reliability  has 
become  a  major  field  of  study  in  structural  engineering 
since  Freudenthal ' s  pioneering  works  (104,  105)  some 
twenty-five  years  ago.  In  the  determination  of  accept¬ 
able  probabilities  of  failure,  the  extended  reliability 
concept  of  Ang  and  Amin  (106,  107)  seems  to  provide  a 
reasonable  basis  for  further  development,  which  incor¬ 
porates  subjective  decision-making  into  objective  prob¬ 
abilistic  analyses.  In  this  regard,  the  decision  anal¬ 
ysis  such  as  the  applications  of  Benjcimin  (108,  109) 
can  be  useful  in  considering  the  economic  consequences 
of  failure. 

8.  Importance  of  Structures  -  It  seems  to  be  desirable  for 
the  design  engineer  to  make  use  of  a  quantitative  as¬ 
sessment  of  the  importance  of  structures.  Obviously, 
the  safety  of  a  school  is  more  important  than  that  of 

a  warehouse.  Perhaps,  the  acceptable  probabilities  of 
failure  can  be  made  as  functions  of  these  "measures  of 
importance"  for  various  structures. 

9 .  Soil-Structure  Interaction  -  The  structural  response  is 
certainly  dependent  on  the  foundation  and  surrounding 
soils.  This  problem  should  also  be  studied  further. 

10 .  Interaction  of  Structural  and  Non-structural  Components 
-  Although  it  is  known  that  non-structural  components 
affect  the  overall  structural  response  to  earthquake 
loads,  these  do  not  appear  to  exist  in  any  concise 
solution  to  this  problem.  It  seems  to  be  desirable  to 
study  interaction  of  structural  and  non-structural  com¬ 
ponents  both  cinalytically  and  experimentally  at  CERL. 

These  recommendations  for  further  study  are  not  meant 
to  be  exhaustive.  Nevertheless,  they  seem  to  be  the  signif¬ 
icant  ones  in  relation  to  the  development  of  a  new  seismic 
code,  with  which  safer  arid  more  economical  building  struc¬ 
tures  can  be  designed. 
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APPENDIX  A 


REVIEW  OF  DISCRETE  MECHANICS 
AND  SOME  ORIGINAL  SOLUTIONS 


General 

The  "discrete  mechanics”  or  "discrete  field  mechanics" 
can  be  used  to  obtain  field  or  functional  solutions  for 
structural  systems,  which  can  be  represented  with  a  lattice 
or  a  network  of  elements.  In  discrete  mechanics,  the  mathe¬ 
matical  model  is  a  set  of  difference  equations  or  recurrence 
relations  among  lattice  nodes .  These  equations  are  then 
solved  in  closed-form  with  the  use  of  the  calculus  of  finite 
differences  (110) ,  which  is  concerned  with  functions  of  in¬ 
teger  variables.  The  calculus  of  finite  differences  is  the 
suitcible  mathematical  tool  in  discrete  mechanics  in  the  seune 
manner  as  the  differential  calculus  is  effective  in  contin¬ 
uum  mechanics . 

From  the  viewpoint  of  continuum  mechanics,  the  infini¬ 
tesimal  elements  of  an  elastic  body  are  used  to  express  com¬ 
patibility  and  dynamic  conditions,  which  along  with  the 
material  properties  lead  to  a  set  of  differential  equations. 
On  the  other  hand,  in  discrete  mechanics,  a  finite  number  of 
elastic  segments  are  examined  to  satisfy  the  compatibility 
and  dynamic  requirements  at  the  joints,  which  yield  a  set  of 
difference  equations  or  difference-differential  equations 
(111-113) .  Functional  solutions  to  these  equations  are  field 
solutions,  which  describe  the  dynamic  behavior  of  the  nodal 
points  in  the  structure. 

For  further  references.  Dean  (114)  gave  a  comprehensive 
outline  of  basic  concepts  in  discrete  mechanics.  Leimback 
and  McDonald  (42)  applied  the  technique  to  investigate 
steady  state  vibration  of  certain  frames.  Wah  (41)  formu¬ 
lated  natural  frequencies  of  vibration  of  grids  vibrating 
out  of  plane  using  the  technique.  Mauivi,  Duvall  and  Lowell 
(112)  formulated  the  wave  propagation  in  lattices.  Lin 
et  al.  (115)  has  adopted  a  technique  resembling  discrete  ap¬ 
proach  to  formulate  ribbed  panel  vibration  problem. 

Deaui  and  Ugarte  (44)  applied  this  technique  to  the  ana¬ 
lysis  of  building  frames.  Omid' varan  (116)  investigated  be¬ 
havior  of  slender  high-rise  frames  using  finite  difference 
approach . 

The  formulation  of  governing  difference-differential 
equations  of  latticed  frames  based  on  energy  approach  was 
given  by  Omid 'varan  (117) .  Omid 'varan  (118)  and  Guthowski 
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(119)  and  Dean  and  Ugarte  (44)  have  all  investigated  the 
stability  of  frameworks.  Space  lattices  and  structural  nets 
were  treated  by  Deem  (120,  121)  and  Dean  and  Ugarte  (122- 
124) .  The  basis  for  investigation  of  behavior  of  composite 
discrete  cuid  continuous  elastic  media  is  given  by  Decui, 

Omid' varan  and  Gangaro  (125,  126).  Additional  application 
of  the  technique  which  more  represents  the  historic  develop¬ 
ment  of  the  technique  is  found  in  references  (127-129) . 

In  the  following,  examples  are  given  to  illustrate  the 
application  of  the  calculus  of  finite  differences  in  dis¬ 
crete  mechanics.  In  addition,  original  solutions  to  several 
problems  in  structural  dynamics  are  given. 


Example  1  -  Continuous  Beams* 

A  closed  functional  solution  was  derived  for  deforma¬ 
tions  and  stresses  at  the  supports  of  a  beam  continuous  over 
N  panels  herein.  The  moments  at  the  supports  represent  a  fi¬ 
nite  number  of  ordinates  of  a  continuous  function,  namely  the 
moment  in  the  continuous  becun.  As  stated  earlier,  discrete 
mechanics  can  be  employed  to  develop  closed  functional  solu¬ 
tions  to  a  finite  number  of  ordinates  of  a  continuous  func¬ 
tion.  Given  a  continuous  beam  loaded  in  an  arbitrary  manner 
as  shown  in  Figure  5,  an  equivalent  system  as  far  as  joint 
deformation  and  support  moment  is  concerned  carries  concen¬ 
trated  moments  m  at  the  supports. 


Figure  5 .  Continuous  Beam 


m^  =  Mj,(r)  +  M^(r-l) 


(79) 


where  m  denotes  the  sum  of  fixed  end  moments  at  support  "r” 
due  to  loading  of  the  panel  preceding  and  the  panel  succeed- 


♦driginally  presented  by  D.  L.  Dean. 
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ing  joint  "r".  Mp  (r)  and  denote  fixed  end  moments 

at  the  support  "r."  The  free  body  diagram  of  the  joint  is 
shovm  in  Figure  6.  ’  , 


Figure  6 . 


The  Free  Body  Diagram  of  Joint  Loaded 
Equivalent  System 


Considering  the  joint  free  body  diagram  at  the  support  as 
shown  in  Figure  6,  we  have  '  . 


M  + 
r 


^r-l 


m^  =  0 


(80) 


The  panel  force-deformations  are; 


-Mr- 

■(E+2)  ' 

=  2k 

M' 

r  _ 

(2E+1) 

(81) 


In  the  cJsove  presentation  k  denotes  the  panel  flexural  stiff¬ 
ness  auid  E  denotes  the  Bool's  difference  operator.  Operators 
such  as  E,  delta  A,  Nebla  V,  Devla  AV,  and  multa/!E7are  among 
the  most  frequent  difference  operators  in  the  calculus  of 
finite  differences  and  are  defined  as  follows: 


E(0j.) 

=  ®r+l'  “  =  ®r-l 

(82) 

A(-) 

=  (E-l)(*);  e.g.,  (A0^=0^^j^-0^) 

(83) 

1 

V(*) 

=  (1-e‘^)(-) 

(84) 
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A7(*) 

2  (E-2+e"^)(*) 

(85) 

ZH7(*) 

2  1  (E-e“^)  (•) 

(86) 

Introduce  Equation  81  into  Equation  80  to  obtain  the  govern¬ 
ing  difference  equation  for  6^.  The  result  for  the  case 
when  panel  stiffness  properties  are  identical  (i.e.,  for  k  = 
constant  and  for  azbitraiy  loading)  is  given  by 

(AV  +  6)0j.  =  ^  (87) 

This  second  order  difference  equation  is  the  counterpart  of 
the  following  differential  equation  in  differential  calculus 

(D^  +  6)  y(X)  =  ~  m(x)  (88) 

2  2  2 

where  D  =  d  /dx  .  The  solution  to  Equation  87  based  on  cal¬ 
culus  of  finite  differences  will  contain  a  homogeneous  and  a 
particular  solution 

0^  =  (-1)^(C2_  sinh  ar  +  C2  cosh  ar)  +  (89) 

where 


a  *  cosh“^(2)  (90) 

The  particular  solution  0p(r)  may  be  evaluated  for  any  load¬ 
ing  condition.  The  most  general  loading  case  is  that  of  a 
unit  impulse  moment  at  r=b.  All  other  solutions  can  be  ob¬ 
tained  by  superposition  of  this  unit  impulse  solution. 

'  =  iv^sinh'  a 

where  H(r-b)  denotes  discrete  step  function  such  that, 

!o,  for  r<b 

(92) 

1 ,  for  r^b 

For  fixed  end  conditions  0  =0=0,  these  constants  are 

given  b)'  ° 


C2  =  -0p(O), 


C  =  . 

1  sinh  aN 


[(-l)^*'”  0p(N) 


+  0n(r)  cosh  aN] 

tr 


(93) 


Given  the  closed  functional  solution  to  0  , 
be  used  to  express  as  follows:  ^ 


Mj.  =  2k  (0J. 


2Qrtl) 


Equation  81  can 
(94) 
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Results  of  this  example  showed  that  discrete  mechsmics  rep¬ 
resents  an  extension  of  finite  element  method.  In  the  fi¬ 
nite  element  method,  the  properties  of  individual  segments 
are  employed  to  express  joint  equilibrium  (or  dynamic  re¬ 
quirements)  and  compatibility  requirements.  The  result  is  a 
set  of  linear  equations  which  are  solved  numerically.  In 
discrete  mechanics  the  segment  properties  are  employed  to 
derive  a  difference  equation  describing  the  behavior  of 
nodal  points  at  the  joints  between  the  segments.  A  closed 
solution  is  derived  based  on  a  solution  to  difference  equa¬ 
tions  based  on  calculus  of  finite  differences . 


Example  2  -  A  Truss 

A  cross-braced  pin  connected  truss  as  shown  in  Figure  7 
is  a  simple  example  of  a  discrete  elastic  body.  Field  solu¬ 
tions  to  stresses  are  given  to  illustrate  our  point  of  view. 
We  shall  formulate  the  stresses  in  the  chord  segment  and  the 
bracing  identified  in  a  sequence  as  the  chord  by  bracing  "r" . 
These  are  denoted  by  and  respectively  as  shown  in  Fig¬ 
ure  8.  In  this  manner,  and  Gj.  are  functions  of  an  inte¬ 
ger  variable  "r". 


Figure  7 .  Cross-Braced  Truss 


Consideration  of  the  equilibrium  requirements  of  the  joint 
"r"  as  shown  in  Figure  8  of  a  symmetrically  loaded  truss 
will  give: 


(Gj.  +  Gj._j^)  cosa  +  =  0 


(95) 


(Gr  -  G^_j^)  sina  +  =  0 


(96) 
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Figure  8.  Free  Body  Diagram  of  Joint  "r" 


Equations  95  and  96  represent  a  set  of  two  first-order 
difference  equations,  which  will  be  solved  for  and  N^. 
For  constant  a  Equation  96  yields. 


AG^  = 


_vr+l 

sin 


where  A  denotes  the  first  foarward  difference  operator,  l.e., 
AG^  =  ®r+l  *" 

The  solution  will  Include  a  homogeneous  solution  which 
In  this  case  Is  a  constant  and  particular  solution  Gx>(r) 
which  depends  on  loading  term  P^: 

Gr  »  +  Gp(r)  (98) 

For  uniform  loading  Pj.  =  P  and 


=  ilHS  ‘’5 

The  constcuit  Ci  Is  determined  In  view  of  the  end  condition 
at  r  =  0. 

R  *  NP  =  G  sina,  or 
o  o 


G  = 
o  slna 


From  Equation  100,  we  obtain 

r  -  PN 
'"1  slna 


(100) 


(101) 


Thus,  a  functional  field  solution  to  the  stress  In  bracing 
"r"  is  given  by 
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G  =  P(N-r) 

sinci 


(102) 


To  obtain  a  closed  solution  to  stresses  in  the  chords ,  Equa¬ 
tions  95  and  102  will  be  considered  From  Equation  95, 

=  -cos  (Gj.  +  (103) 


In  view  of  Equation  102, 

AN^  =  -P  cota  (2N-2r-l) 


(104) 


The  solution  to  Equations  104  will  include  a  summation  con¬ 
stant  C2  and  a  particular  solution,  as  follows. 

Nr  =  C2  +  P  cota(r^-2Nr)  (105) 

Using  the  condition  N^  =  -PN  cota, 

C2  =  -PN  cota  (106) 

As  a  result,  field  solution  to  stresses  in  the  chords  is 
fo\ind  to  be: 


N^  =  P  cota(-N+r  -2Nr) 


(107) 


Since  the  truss  considered  in  this  example  is  a  static¬ 
ally  determinate  structure,  the  stiffness  properties  of  the 
individual  segments  did  not  have  to  be  considered  in  formu¬ 
lation  of  the  stresses .  The  stiffness  properties  become 
necessary  should  we  choose  to  formulate  a  field  solution  to 
truss  deformations . 


The  following  example  will  illustrate  the  application 
of  discrete  mechanics  in  the  dynamic  analysis  of  a  static¬ 
ally  indeterminate  system. 

Example  3  -  Steady  State  Vibration  of  a  Lumped  Beam* 

The  frequencies  of  vibration  of  a  system  of  N  concen¬ 
trated  equal  masses  connected  by  massless  elastic  bars  of 
equal  lengths  as  shown  in  Figure  9  will  be  formulated  to  il¬ 
lustrate  the  dynamic  application  of  discrete  mechanics.  A 
multi-story  shear-frame  in  vibration  remotely  resembles  a 
lumped  cantilever  beam. 

The  free  body  diagram  of  the  mass  M  at  the  node  point 
"r"  is  shown  in  Figure  10.  The  notation  employed  to  desig¬ 
nate  the  moment  and  shear  in  the  elastic  segment  preceding 
the  r*^  mass  is  shown  in  Figure  11. 


♦Prepared  by  C.  Omid' varan  in  collcdjroation  with  C.  Bacchus. 
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Figure  9.  Lumped  Mass  Becun 


Figure  10,  The  Free  Body  Diagrcun  of  Mass  r 


Figure  11.  The  Free  Body  Diagram  of  Panel  r 
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The  lateral  displacement  and  rotation  of  mass  r  are  denoted 
by  Y^(t)  amd  0j.(t)  respectively,  and  Y^(t)  and  ©^(t)  are 
discrete  ftinctions  of  r  and  continuous  functions  of  t,  the 
time  element. 


Neglecting  rotational  inertia,  the  dynamic  consideration 
of  mass  r  as  shown  in  Figure  10  will  require  that: 


M'  +  M  =  0 

r  r+1 

(108) 

<i^Y^(t) 

''r  - 

(109) 

The  stiffness  properties  of 
masses  are  given  by: 

elastic  bars  connecting  the 

Mr  =  k[40r.i  +20^-1 

(110) 

«r  =  ''f20r.i  ^  4ar  -  ^ 

(111) 

Vr  =  V'  =  i  (Mr  +  Mp 

(112) 

where  k  =  EI/L. 


Eliminate  the  joint  rotations  from  Equations  110  and 
111  for  a  difference  equation  relating  M^,  and  Yj..  The 
result  may  be  introduced  into  Equation  108  to  obtain: 


(AV^+6)M^  - 


AVr^r  =  0 


(113) 


The  difference  operation  AVj.  (Deb la)  ,  is  an  operator 
concerning  typical  three-node  relations. 


^^r^r  =  ^r-1  '  ^^r  ^  ^rfl 
AVrM-  .  M-_^  -  2m;  + 


(114) 

(115) 


Equation  113  represents  a  generalized  three-moment  equation 
for  non-zero  displacement  at  the  node  points. 

In  view  of  Equation  108,  Equation  112  will  give. 


V  =  V  =  f  (M'-M*  ,) 
r  r  L  r  r-1 


(116) 


Introduce  Equation  116  into  Equation  109,  to  obtain. 


d^Y. 


AV  M'  +  mL 
r  r 


=  0 


dt^ 


(117) 
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This  is  a  second  order  difference-differential  equation.  For 
a  steady  state  harmonic  motion  with  angular  frequency  o) ,  we 
have 


ri  = 


(118) 


As  a  result  of  Equation  118,  Equation  117  can  be  transformed 
into  the  following  difference  equation: 


(119) 


The  solution  to  a  set  of  two  difference  equations  as  given  by 
Equations  113  and  118  for  and  represents  a  closed  func¬ 
tional  solution  to  the  problem.  In  a  matrix  form,  these  are 


(AVj.  +  6) 


-mLo) 


(120) 


For  simple  support  conditions  at  r=0  and  r=N,  both  M*  and  Y 
must  vanish  at  the  boundaries. 

The  following  finite  series  will  satisfy  the  required 
end  conditions  naturally. 


(121) 


^  Jo 


(122) 


Introduce  Equations  121  and  122  into  Equation  120,  we  obtain 


2 (2+cos6^) 


2(cosB^  -  1) 


where 


-12  ^  (cose.  -  1) 

ii  X 


-mLu) 


=  0  (123) 


fl  _  Iff 

®i  ■  iT 


(125) 


For  a  non-trivial  solution  the  determinate  of  coefficients 
must  be  identical  with  zero,  and  we  obtain, 
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The  N  different  angular  frequencies  of  vibration  of  a  lumped 
becua  are  given  by  Equation  125.  The  cuigular  frequencies  are 
found  to  approach  those  of  continuous  beams  rapidly  as  the 
number  of  lumped  masses  replacing  a  continuous  beam  increase. 
For  10  equal  masses  w  =  9.87  in  contrast  with  for  a  con¬ 
tinuous  beam  given  that  k/mL2  =  l. 


Vibration  of  Multi-Story  Frames 


A  single-bay  multi-level  fraitfewith  uniform  properties  as 
shown  in  Figure  12  is  used  to  illustrate  the  derivation  of  a 
closed  functional  solution  to  dynamic  properties  and  response 
of  structures.  This  work  illustrates  the  effectiveness  of 
discrete  mechanics  in  resolving  several  problems  in  the  im¬ 
plementation  of  a  new  code. 


Figure  12 .  Discrete  Model  of  a  Multi-Story  Frame 

The  masses  are  assumed  to  be  identified  and  denoted  with  a 

sequence  of  numbers  ...  r-1,  r,  r+1  _ _  The  forces  acting 

on  the  mass  identified  as  mass  "r"  in  Figure  12  include  the 
inertia  force  Fj,  the  damping  force  F^,  the  elastic  force 
Fg,  and  the  forcing  function  P(t). 

t 
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Considering  the  free  body  of  mass  x,  we  have 

-Fj  +  Fjj  +  Fg  =  P(t)  (126) 

Denoting  the  displacement  of  this  mass  by  V  (t) ,  the  inertia 
force  is  given  by 

Fj  =  -mVj.(t)  (127) 

The  elastic  force  on  mass  ”r"  will  depend  on  displacement  of 
mass  (r-1)  preceding  and  mass  (r+1)  following  this  mass  shown 
in  Figure  14. 


Figure  14.  Shearing  Forces  on  Mass  r 


Let  k  denote  the  frame  shearing  stiffness  relating  the  shear¬ 
ing  force  and  relative  end  displacements  of  any  two  succes¬ 
sive  levels. 


Sr  +  s;  =  k(V^^^  -  v^) 

Then  the  elastic  force  is  given  by 

Fg  =  .  k(V^(t)  -  - 


K[Vr.i(t)  -  Vj.(tH 


Thus 


Fg  =  -kAV^Vj.(t) 


(128) 


where  AV  ,  Debla,  denotes  the  second  central  difference 
operatorf  such  that 

AVrVr(t)  H  -  2Vj.(t)  +V^_j^(t)l 


(129) 


(130) 


(131) 


The  damping  force,  Fq,  will  in  a  similar  memner  depend  on 
damping  coefficient  C  and  second  central  difference  of  velo¬ 
cities  at  three  nodes 
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Fjj  =  -CAV^Vj.(t) 


(132) 


iitroducing  Equations  127,  130,  and  132  into  Equation  126,  we 
obtain : 

mV^(t)  -  kAVj.V^(t)  -  CAVj,V^(t)  =  P(t)  (133) 

The  second-order  difference-differential  equation  as  given  by 
Equation  133  represents  the  mathematical  model  for  the  frame 
considered.  It  is  a  differential  equation  with  respect  to 
continuous  element  t,  and  a  difference  equation  with  respect 
to  integer  varic±>le  r. 


First  the  steady  state  harmonic  motion  will  be  investi¬ 
gated.  A  solution  will  be  given  for  an  undamped  structure 
(C  =  0)  subject  to  earthquake  excitation  (P(t)  =  0) .  Then 
Equation  133  will  reduce  to 

mV^(t)  -  kAV^V^(t)  =  0  (134) 


In  view  of  harmonic  nature  of  structure  in  the  absence  of 
damping  we  may  assume  that, 

Vr(t)  =  F^  sincjt  (135) 

where  u  denotes  the  angular  frequency  of  vibration.  Substi¬ 
tuting  Equation  135  into  Equation  134,  we  obtain, 


(AV^-23)F^  =0;  6  =  -  ^ 

The  solution  to  Equation  136  is  found  to  be: 

F^  =  A  sinyr  +  B  cosyr,  for  (-2  <  3  <  0) 


(136) 

(137) 


F^  =  (-1)^  [Aj^  coshnr  +  A2  sinhnr] ,  for  (3  <  -2)  (138) 


where 


y  =  cos”^(l+3),  n  =  cosh“^{-l-3)  (139) 

The  siimmation  constants  A,  B,  Ai  and  A2  will  first  be  de¬ 
termined.  Then  Equation  135  will  be  used  to  express  func¬ 
tional  solution  to  displacement  V^(t) .  For  the  purpose  of 
illustration,  the  solution  corresponding  to  (-2  <  3  <  0) 
only  is  presented  herein. 


For  a  steady-state  vibration  of  a  "N"  story  frame  fixed 
at  r  =  0,  the  end  conditions  are  assumed  to  be: 


V^(t)  =  0  (140) 

=  0  (141) 
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V^Ct)  =  0  (142) 

[V^^l(t)  -  =  0  (143) 

Equation  143  implies  that  the  shear  at  r  =  N  vanishes.  Equa¬ 
tion  140  will  given  B  =  0.  Then  Equations  141  and  142  are 
also  satisfied.  Equation  143  will  give: 

A[sinY  cosyN  +  (cosy-l)  sinyN]  =  0  (144) 

For  a  nontrivial  solution  (A  ^  0)  we  must  have 

siny  cosyN  +  (cosy-1)  sinyN  =  0  (145) 


Equation  145  is  a  function  of  three  parameters,  N,  k/M,  and 
u),  angular  frequency.  A  computer  program  can  be  employed  to 
plot  0)  versus  k/M  for  different  choices  of  N.  Such  a  graphi¬ 
cal  solution  is  rationally  based  and  reflects  considerable 
improvement  over  existing  solutions  to  frame  natural  fre¬ 
quencies  . 

A  functional  solution  to  the  shear  force  at  differential 
levels  of  a  multi-story  frame  due  to  a  unit  impulse  accelera¬ 
tion  at  the  base  can  also  be  formulated.  The  response  due  to 
any  particular  earthquake  acceleration  is  the  result  of  inte¬ 
gration  of  response  due  to  unit  impulse  load. 


The  boundary  conditions  for  a  unit  impulse  acceleration 
acting  at  the  base  of  a  multi-story  frame  are  assumed  to  be; 


o 

n 

> 

(146) 

v^(o)  =  1 

(147) 

(148) 

Equation  146  is  naturally  satisfied  in  view  of  Equation 
135.  Equations  147  and  148  yield; 


-  _  [siny  sinyN  -  3  cosyN] 

[siny  cosyN  -  3  sinyN] cu 

Given  summation  constants  A  and  B,  a  closed  functional 
solution  to  displacement  response  Vj-(t)  is  obtained,  and  we 
will  be  in  a  position  to  formulate  shear  response  spectrum. 
The  shearing  force  Qr(t)  is  equal  to  the  total  sum  of  in¬ 
ertia  forces  acting  on  the  masses  above  the  rth  level. 


(149) 

(150) 


Q_(t)  =  -  I  mw'^V.  (t) 
i=r+l  ^ 


(151) 
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In  view  of  Equation  135,  Equation  151  will  give, 

2  N 

Q  (t)  =  -mo)  sin  t  J  F.  (152) 

^  i=r+l  ^ 

Calculus  of  finite  difference  (110)  provides  means  of  evalua¬ 
tions  of  finite  sums  of  discrete  functions  in  the  manner  dif¬ 
ferential  calculus  helps  with  integrals  of  continuous  func¬ 
tions.  These  will  be  employed  to  simplify  Equation  152. 

Q^(t)  -mw^  sinoJt  [aT^  ^i^i=r+l  (153) 

which  will  give: 

2 

Q  (t)  =  [-A  COSY  (i-1/2) 

2  sinl 

+  3  siny  (1-1/2)  (154) 

or 

2 

Q  (t)  =  -  [A  COSY  (N+ 1/2)  -  A  cosy  (r+1/2) 

^  4  sinX 

+  3  siny  (N+1/2)  -  3  siny  (r+1/2) ]  (155) 

Equation  155  represents  a  solution  to  the  shear  force  at  any 
level  as  afunction  of  time  due  to  a  tinit  impulse  acceleration 
at  the  base.  In  addition,  it  can  be  used  for  design  pur¬ 
poses,  taking  into  account  factors  which  are  neglected  at  the 
present  time. 


Torsion  and  Biaxial  Shear 


Inspection  of  buildings  subjected  to  earthquakes  has 
indicated  the  presence  of  torsional  oscillations  and  tor¬ 
sional  damage  in  these  structures.  Most  structures  are  not 
symmetric  in  mass  and/or  stiffness.  Further,  shearing  stiff¬ 
ness  properties  of  structures  are  different  along  principal 
stiffness  axes  in  normal  cases.  The  nonsymmetrical  features 
will  effect  the  induced  stresses  subjected  to  a  dynamic  ex- 
citaticn  in  an  arbitrary  direction.  Biaxial  shear  displacement 
as  well  as  torsion  will  result  in  structures. 

In  the  cas«_  ^f  one  degree  of  freedom  systems  it  is  pos¬ 
sible  to  assume  an  equivalent  static  force  acts  through  the 
center  of  mass  of  structure  and  determines  the  forces  based 
on  elastic  analysis.  However,  for  multi-degree  of  freedom 
systems,  it  is  difficult  to  conceive  of  a  static  method  of 
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analysis  cind  there  is  a  dearth  of  rational  solutions  for 
torsional  oscillation  of  multi-story  frames,  and  little  is 
known  regarding  combined  torsion  and  biaxial  shear  behavior 
of  high  rise  structures.  Consciousness  of  their  importance 
has,  nevertheless,  led  to  adoption  of  provisions  in  most 
codes  and  specifications.  Housner  and  Outinen  (67)  have 
investigated  the  dynamic  behavior  of  single-story  structures 
in  torsion.  Bustamante  and  Rosenblueth  (68)  have  made  a 
general  examination  of  code  provisions  and  available  litera¬ 
ture.  Shiga  (69)  has  presented  a  manner  of  formulation  of 
dynamic  behavior  of  multi-story  frames  in  torsion  in  cases 
when  the  centers  of  mass  and  rigidity  coincide.  However, 
the  method  is  not  easily  applicable  since  it  requires  pre¬ 
determination  of  rigidity  ratios  in  shear  and  torsion  of  each 
level  in  relation  to  all  other  levels. 

The  presentation  here  illustrates  a  method  of  formula¬ 
tion  of  dynamic  properties  and  response  in  torsion  and  bi¬ 
axial  shear.  The  illustration  assumes  identical  mass  and 
geometry  for  all  levels. 

The  displacements  of  level  "r"  are  designated  by  u^,  v^, 
and  0J..  Ur  and  v^  denote  displacement  parallel  with  x  and  y, 
the  principal  elastic  axis.  0r  denotes  the  angle  of  distor¬ 
tion  and  for  convenience  will  be  transformed  into  a  parameter 

Wr: 

^  =  9^  (156) 

R  =  radius  of  gyration  of  mass  m  of  each  level. 


2 


X 


Figure  15.  Principal  Elastic  Axis 

First  influence  coefficients  will  be  defined  for  dis¬ 
placement  and  force  interaction  between  any  two  successive 
levels,  kx  will  be  used  to  designate  the  force  in  x  direc¬ 
tion  of  any  one  floor  due  to  a  unit  displacement  in  y  di¬ 
rection  for  a  rotation  corresponding  to  unit  w  of  following 
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floor,  kv  will  denote  torque  in  z-direction  for  a  rotation 
corresponding  to  unit  w  of  following  panel. 


j^x  _  reactive  force  direction 

y  ~  unit  displacement  direction 


,(157) 


Given  ,  the  ratio  of  rigidity  against  shear  in  y  direction 
to  that  in  x  direction;  we  have 


=  (Xk* 

y  X 


(158) 


In  view  of  the  choice  of  x  and  y  as  principal  elastic  axes' 
we  have 


k^  =  k^  * 
Y  X 


(159) 


Equation  159  will  serve  as  the  basis  for  choice  of  principal 
axes  X  and  y.  6»  Y»  ©2  will  be  used  to  designate 

other  influence  coefficient  ratios  as  defined  by  the  follow¬ 
ing  relations: 


k"-  = 

X 

k.^  = 


®2^z 


k^ 

z 


6k 

Yk 

nk 


(160) 

(161) 

(162) 

a63) 

(164) 


In  the  absence  of  damping,  inertia  and  elastic  forces 
in  X,  y  and  z,  directions  will  be  active.  The  inertia  forces 
are  given  by: 

(165) 


*• 

'  -"“r 

ly  ’ 


(166) 

(167) 


The  shearing  elastic  forces  acting  on  panel  r  are  a  function 
of  displacements  of  the  mass  preceding  the  mass  succeeding 
this  panel.  Vj^,  the  shearing  force  in  x  direction,  is  given 
by  Figure  16. 

-  S'.,  (168) 


In  view  of  influence  coefficients  defined  earlier,  we  find 
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3 


Figure  16 i  Shearing  Forces  in  x  Direction 


(169) 

r-1  =  <  "^r  ^  "^r 

(170) 

Where  A  and  7  denote  the  first  forward  and  backward  dif¬ 
ference  ot^erators 


u 


r+1 


u 


r 


u 


r-1 


(171) 


Considering  Equation  159,  and  introducing  Equation  169 
and  170  into  Equation  168,  we  obtain 


Vx  =  k*  [AV  u^  +  yAVw^]  (172) 

The  shearing  force  V  in  the  y-direction  is  found  in  a 
similar  manner  to  be:  ^ 

Vy  =  ak*  [AVv^  +  riAVw^]  (173) 

The  torque  T^  is  found  to  be: 

T^  =  3k*  [AVWj.  +  e2AVVj.  +  e^^AVu^]  (174) 

It  will  be  assumed  that  the  principal  mass  and  elastic  axes 
are  parallel  as  shown  in  Figure  17. 


Denoting  the  eccentricities  between  center  of  mass  c 
and  elastic  center  G  by  a^  and  a2,  consideration  of  the  faree 
body  diagram  of  mass  r  as  shown  in  Figure  17  will  require 
that : 


I 


75 


Equations  178-180  represent  a  set  of  three  difference- 
differential  equations  for  u^r  v^./  and  w^..  The  mathematical 
model  for  torsion  combined  with  biaxial  shear  of  a  shear 
£r2une  is  a  second  order  differential  equation  with  respect 
to  time  element  and  second  order  difference  equation  with 
respect  to  displacement  elements. 

In  view  of  the  fact  that  damping  is  assumed  to  be  ab¬ 
sent,  we  have. 
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Uj.  =  u^q(t),  Vj.  =  v^q(t),  Wj.  =  Vl^q(t) 


(182) 


where 

q(t)  *  A  coswt  +  B  sinujt 
Considering  Equations  182  and  183,  we  find 


(183) 


U_  =  -0)  U  ,  V  =  -0)  V  ,  w  =  -u  w 
r  r  T  x'  X  r 


(184) 


Equation  184  will  serve  to  separate  difference  and  differen¬ 
tial  mathematical  models.  In  view  of  Equation  184,  Equations 
178-180  will  give: 


(AV  +  Moj^) 


(aAV  +  Md)  ) 


anAZ 
Aj^AS/  + 


=  [0]  (185) 


A  solution  to  vector  difference  Equation  185  to  meet  physical 
boundary  conditions  of  multi-story  frame  is  feasible  by  the 
aid  of  the  stress  ^unction  approach.  For  fixed  end  condi¬ 
tions  at  r=0  and  r=N,  a  solution  is  easily  feasible  and  will 
be  given  to  conclude  this  illustration.  Fixed  end  conditions 
can  be  met  naturally  through  assumption  of  the  following  fi¬ 
nite  series  representation  for  displacements: 


%  =  •'k  =  , 


(186) 


Introduce  Equation  183  into  Equation  182,  where  we  obtain  the 
following: 


(2  cosbj^-2) +Ma)^ 


2A2 (cosb^-1) 


2a (cosb^-1) +Mu^ 
2A^ (cosb^-1) 


2y (cosb^-1) 
2an (cosb^-1) 
2Aj^  (cosbj^-1)  +Mu)' 


=  [0] 


i=l, . . . ,N) 


(187) 


The  N  different  frequencies  of  vibration  of  frame  (i=l 
to  N)  in  torsion  and  biaxial  shear  are  contained  in  the  solu¬ 
tion  to  the  linear  equation  derived  from  the  expansion  of  the 
determinate  in  Equation  187 . 


The  strength  of  the  technique  in  discrete  mechanics  il¬ 
lustrated  in  these  examples  will  be  more  fully  appreciated 
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considering  the  fact  that  torsional  vibration  of  frameworks 
is  recognized  as  one  of  the  more  complex  problems  in  earth- 
quate  engineering.  Scarcity  of  literature  regarding  the 
subject  is  an  evidence  of  this  point  of  view.  This  illus¬ 
tration  has  been  concerned  with  more  than  torsional  vibra¬ 
tion  in  formulating  the  case  of  torsion  combined  with  bi¬ 
axial  shear  displacement. 
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APPENDIX  B 


STATISTICAL  METHODS 


General 


There  exist  vincertainties  in  the  prediction  of  the  oc¬ 
currence  of  future  earthquakes  as  well  as  the  estimation  of 
the  resistance  of  structural  materials.  Therefore,  it  is 
desirable  to  apply  statistical  methods  to  problems  of  earth 
quake  engineering. 

The  objective  of  this  appendix  is  to  review  and  summa¬ 
rize  available  literature  in  this  regard.  The  terminology 
of  statistics  is  introduced  in  this  section  as  well  as 
wherever  it  is  necessary.  The  statistical  characteristics 
of  strong-motion  earthquakes  and  their  computer  simulations 
are  then-  discussed.  The  random  response  of  seismic  struc¬ 
tures  is  reviewed.  Finally,  the  general  problem  of  struc¬ 
tural  safety  in  earthquake  engineering  is  described  and  dis 
cussed,  along  with  a  brief  introduction  to  decision  anal¬ 
ysis. 


Although  the  magnitude  of  an  earthquake  can  be  repre¬ 
sented  with  the  use  of  the  Richter  scale  and  the  intensity 
can  be  described  with  the  use  of  the  Modified  Mercalli 
scale,  the  most  basic  and  useful  data  of  earthquake  engi¬ 
neering  are  the  recordings  of  ground  accelerations  during 
strong-motion  earthquakes  (130).  The  record  of  an  indi¬ 
vidual  earthquake,  say  x(ii)j^,t),  can  be  considered  as  a  sam¬ 
ple  function  of  a  random  process  X((jj^,t)  ,  the  definition  of 

which  along  with  some  other  terminologies  are  introduced 
herein  (131) . 

A  random  process  (or  stochastic  process  or  random 
function)  is  a  family  of  random  variables  indexed  by  an  ar¬ 
gument  oj  in  a  sample  space  n  and  a  parameter  t  in  an  index 
set  T,  i.e.,  X(t)  =  {X((i)  ,t)  :  weft,  teT}.  In  earthquake 

engineering  applications,  the  parameter  t  denotes  time. 

For  a  fixed  value  u  in  Q,  say  X(u)j^,*)  =  X(a)j^)  is  a 

function  of  t  representing  a  possible  observation  of  the 
random  process.  The  function  X(u)j^)  is  called  a  realization 

or  a  sample  function  of  the  process.  For  a  fixed  value  of 
t,  say  tj^,  X(t^,*)  =  X(tj^)  is  a  random  variable,  which  is 
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a  function  on  the  sample  space  il  such  that  a  specific  value 
is  assigned  to  each  and  every  outcome  w.  The  probability 
that  a  random  variable  is  less  than  or  equal  to  a  specific 
value  X  is  called  the  distribution  function,  i.e., 

Fx(t^)(x)  =  P(X(tj^)  <  X)  (188) 

Equation  188  is  called  the  first-order  distribution  of  the 
random  process  X(t)  with  its  corresponding  density  function 
as  follows: 


Px(tj^)  3x 


(189) 


Given  two  times  instances  t^^  and  t2,  there  are  two  random 
variables  X(tj^)  and  X(t2).  The  second-order  distribution 
and  density  functions  are  respectively  as  follows: 

^X(t^)X(t2)  ~  P(X(t^)  ^  1x2)  (190) 

32 

Px(tj^)X(t2)  ”  3x]Tx^  ^X(tj^)X(t2) 

The  description  of  a  random  process  can  also  be  made  by 
the  following  moment  functions: 


yx(ti)  =  E[X(tj^)]  =  /  X  p. 

«00 

<|)xx(ti,t2)  =  E[X(t^)X(t2)  ] 


X(t^) 


(x)  dx 


(192) 


-  /  fxj^x^  Px(tT  )X(t,)  dx^dx2 

—00  1  ^ 


(193) 


where  the  first  moment  second  moment  ‘l>xx^^l'^2^ 

are  called  the  mean-value  function  and  the  autocorrelation 
function,  respectively.  The  averaging  operation,  E,  de¬ 
notes  expectation.  Another  important  second-order  moment 
is  the  autocovariance  function  defined  as  follows: 


^XX^^l'^2^  =  E([X(tj^)  -  [X(t2)  -  yx(t2)]) 

If  for  a  random  process  X(t),  for  all  values  of  n, 


(194) 
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Px(t,)X(t,)  ...  X(t  )  **" 

X  z  n 

“  Px(t,+a)  X(t,+a)  ...  X(t  +a)  •**' 

1  z  n 

(195) 

the  process  is  called  completely  stationary  or  strongly  sta¬ 
tionary  or  stationary  in  the  strict  sense.  If  Equation  195 
is  appliccdjle  for  n=l  cmd  n~2,  the  process  is  said  to  be 
weaJcly  stationary  or  stationary  in  the  wide  sense.  Then  the 
expected  value  is  a  constemt  and  the  auto  correlation  is  de- 

-  t^ .  Al¬ 
though  seismic  ground  motions  are  nonstationary,  they  are 
sometimes  represented  with  stationary  random  processes  for 
the  sake  of  simplicity,  as  it  will  be  discussed  later. 

For  a  weakly  stationary  process  X(t),  a  frequency  de¬ 
composition  of  the  autocorrelation  function 

made  as  follows : 


noted  by  which  is  a  function  of  x  =  t2 


exp(iu)T)  dw  (196) 

00 

1  " 

“^XX^*^^  “  ^  ^XX^^^  exp(-ita)T)  dx  (197) 

«»ao 

If  X(t)  is  a  real-valued  process,  is  a  non-negative 

and  even  function  of  w.  Equations  196  and  197  are  usually 
called  the  Wiener- Khintchine  relations.  Note  that, 

<30 

Rxx(O)  =  I  ^XX^“^  "  E[X^(t)]  (198) 

—00 

Equation  198  cein  be  interpreted  to  mean  that 

scribes  the  density  of  the  total  mean-square  value  over  the 
frequency  domain.  Therefore,  the  function  called 

mean-square  spectral  density.  It  is  also  called  the  power 

2 

spectrum  because  E[X  (t) ]  is  often  a  measure  of  average 
energy. 


Statistics  of  Strong-Motion  Earthquakes 

Housner  (19)  represented  the  earthquake  as  a  series  of 
impulses  which  are  random  in  time.  The  ground  acceleration 
X(t)  is  given  by 
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(199) 


X(t)  =  I  Vj  6(t  -  jtg  -  tj) 

where  is  a  random  variable  representing  the  magnitude  of 
the  impulse,  t^  is  an  average  period  computed  by  divid¬ 

ing  the  total  time  by  the  total  number  of  pulses,  jtjj  +  t^ 
is  the  time  elapsed  between  the  beginning  time  and  the  ar- 

rival  of  the  pulse,  and  6 (t)  is  the  unit  impulse  func¬ 
tion.  The  random  variable  V.  is  specified  with  given  mean 

2  ^ 

E[V]  and  mean-square  E [V  ].  Later,  Thomson  (20)  showed  that 
the  power  spectrum  of  a  random  process  as  given  by  Equation 
199  and  with  E[V]  =  0  is  as  follows: 


4>(a)) 


g  (oj) 


2 


E[V^] 


(200) 


where  g((D)  is  the  Fourier  transform  of  the  pulse  shape. 

When  the  pulse  duration  is  short  in  comparison  to  the  period 
of  vibration,  !g(to)|  approaches  unity  irrespective  of  the 
pulse  shapes.  Consequently,  the  power  spectrum  was  foxind 
to  be  proportional  to  the  weighted  sum  of  the  mean-square 
values,  which  is  the  same  as  the  mean-square  value  of  all 
the  pulses.  In  other  words,  the  power  spectrum  in  this 
case  becomes  a  constant  regardless  of  the  pulse  shape  and 
time  spacing. 


Goodman,  Rosenblueth,  and  Newmark  (61)  used  a  velocity 
pulse  model,  in  which  the  pulses  are  either  uniformly  or 
randomly  occurring  in  time.  Later  Rosenblueth  (132) ,  and 
Rosenblueth  and  Bustamante  (133)  used  an  acceleration  pro¬ 
cess  which  was  generated  from  closely  spaced  velocity  steps 
of  random  sign  and  magnitude  with  a  constant  intensity  per 
unit  of  time. 


When  the  power  spectrum  of  a  random  process  is  found 
to  be  constant,  the  process  is  called  "white  noise"  because 
the  energy  content  of  the  process  is  considered  to  be  uni¬ 
formly  distributed  over  the  entire  frequency  range.  By¬ 
croft  (21)  proposed  such  a  white  noise  with  a  constant 

2  4 

spectral  density  of  0.75  ft  /sec  cps  and  a  duration  of  30 
seconds  for  the  representation  of  typical  strong-motion 
earthquakes.  The  basis  of  comparison  for  this  study  is  the 
response  spectrum  obtained  by  Housner  (134).  Ward  (26) 
adopted  a  similar  white  noise  process  tor  analog  simula¬ 
tions  with  some  consideration  to  emphasize  power  at  certain 
frequencies. 
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Because  observations  showed  that  certain  frequencies 
predominate  the  ground  motion  for  given  soil  conditions, 
Kanai  (135)  obtained  a  response  process  of  the  soil  over¬ 
burden  to  a  white  noise  excitation  at  bedrock.  Tajimi  (136) 
suggested  determining  spectral  density  from  spliced  and  re¬ 
peated  sections  of  the  actual  earthquake  records.  Caughey 
and  Stumpf  (137)  used  this  type  of  model  and  found  root- 
mean-square  response  to  the  1940  El  Centro  earthquake. 


Bogdanoff,  Goldberg,  and  Bernard  (22)  studied  the  use 
of  a  nonstationary  random  process  in  representing  earth¬ 
quake  accelerations  as  follows: 


X(t) 


n 

I  t  a . (exp [-a . tl )  cos(3.  t  +  Y^)»  t  ^  0  (201) 

j=l  3  J  33 


where  a . 


^1  < 


aj. 


and  Bj  are  given  sets  of  positive  numbers  with 
6_ ,  and  are  n  independent  random  vari- 


'2  •••  "n' 

ables  uniformly  distributed  over  an  interval  [0,2Tr].  The 
mean  value  is  zero  and  its  variance  function  has  a  shape  as 
expected  from  an  earthquake  record.  Later,  Goldberg,  Bog¬ 
danoff,  and  Sharpe  (23)  used  the  same  process  except  making 
B.  also  independent  random  variables,  each  uniformly  dis¬ 


tributed  over  the  interval  [6,46].  This  modification  causes 
little  effect  in  the  appearance  of  sample  functions,  but 
there  seems  to  be  considerable  difference  in  their  response 
spectra. 


Lin  (24)  proposed  the  use  of  a  nonstationary  random 
process ,  which  is  obtained  by  passing  a  nonstationary  shot 
noise  through  a  linear  filter.  The  resulting  process  is 
also  called  a  filtered  Poisson  process  (138,139).  A  shot 
noise  S(t)  refers  to  a  series  of  random  pulses,  the  arrivals 
of  which  are  defined  by  a  Poisson  process  N(t),  i.e., 

N(t) 

S(t)  =  I  Y  •  6(t  -  T.)  (202) 

k=l  ^  ^ 

where  Yj^  are  mutually  independent  random  variables  and  5  ( • ) 
denotes  the  unit  impulse  function.  Furthermore, 

E(S(t)]  =  0  (203) 

E[S(t^)S(t-)  1  =  I(tj)5(t2  -  (204) 

In  addition  to  generating  an  earthquake-like  process,  the 
linear  filter  can  also  be  used  in  matching  the  frequency 
content  of  the  process  to  that  of  recorded  ground  motions. 
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Because  it  is  frequently  difficult  to  obtain  analytical 
solutions  for  structural  problems  in  earthquake  engineering, 
statistical  answers  can  be  obtained  by  the  simulation  of 
these  problems  with  the  use  of  digital  or  analog  computers. 
Simulation  techniques  for  earthquake  motions  as  stationary 
processes  were  developed  by  Housner  and  Jennings  (25)  using 
a  digital  computer  and  by  Ward  (26)  using  an  analog  com¬ 
puter.  Nonstationary  earthquake  motions  have  been  simulated 
with  the  use  of  digital  computers  by  Shinozuka  and  Sato 
(29),  Amin  and  Ang  (27),  Jennings,  Housner  and  Tsai  (28), 

Hou  (32) ,  Rascon  and  Cornell  (30) ,  Iyengar  and  Iyengar  (33) , 
and  Levy  (31).  As  it  is  summarized  by  Levy  (31)  ,  the  typi¬ 
cal  operation  of  the  digital  simulation  is  (a)  to  generate 
an  equally  spaced  time  sequence  of  uncorrelated  random  num¬ 
bers  by  drawing  numbers  following  a  normal  distribution  with 
zero  mean,  (b)  to  multiply  these  numbers  by  a  time-dependent 
multiplier,  and  (c)  to  filter  these  numbers  with  a  given  im¬ 
pulse  response  function  to  obtain  the  sample  function  of  the 
earthquake-like  acceleration. 

In  earthquake  engineering,  passive  analog  computers 
have  been  used  by  various  investigators  including  Alford, 
Housner  and  Martel  (140)  ,  Naka,Kato,  and  Yuasa  (141)  ,  and 
Murphy,  Bycroft,  and  Harrison  (142).  Bycroft  (21)  used  an 
active  analog  system  with  amplifiers  in  studying  the  re¬ 
sponse  of  linear  as  well  as  elasto-plastic  single-degree- 
of-freedom  systems  to  white  noise  excitation.  Other  inves¬ 
tigators  using  the  analog  computer  in  earthquake  simulations 
include  Merchant  and  Hudson  (143),  Cherry  and  Brady  (144), 
Brady  '145) ,  Ward  (26)  ,  and  Wirschinq  and  Yao  (146) .  A 
comprehensive  treatment  of  random-process  simulation  in 
general  can  be  found  in  Korn  (147) . 


Random  Response  of  Seismic  Structures 

Until  approximately  10  years  ago,  the  principal  appli- 
:ations  of  random  vibration  studies  have  been  to  missiles, 
satellites,  and  other  space  vehicles.  Recently,  the  ran- 
lom  vibration  of  seismic  structures  has  been  emphasized, 
rrespective  of  the  type  of  structures  studied,  the  ulti- 
late  objective  of  the  study  of  random  vibration  is  to  en¬ 
ure  the  maximum  reliability  of  the  structure.  Some  as¬ 
sets  of  structural  reliability  will  be  discussed  in  the 
ext  section.  In  this  section,  the  fundamentals  of  random 
ibration  will  be  reviewed  along  with  its  application  in 
arthquake  engineering. 

When  the  excitation  can  be  represented  by  a  random 
rocess,  the  problem  of  computing  certain  statistics  of  the 
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response  belongs  to  the  study  of  random  vibration.  During 
the  past  decade,  many  text  or  reference  books  concerning 
this  subject  area  have  been  published  (131,  148-153). 


A  linear  system  can  be  characterized  by  its  response  to 
(a)  the  sinusoidal  forcing  function,  and  (b)  the  impulsive 
forcing  function.  As  an  example,  consider  a  mechanical  sys¬ 
tem  with  the  following  equation  of  motion: 


*.  2  f(t) 

X  +  2^0)  X  +  u^x  =  — — 

o  o  m 

where 


k 

m 


^  2(»)  m 

o 

m  =  mass 

X  =  displacement 
c  ==  damping  constant 
k  =  spring  constant,  and 
f(t)  =  forcing  function. 

When  the  forcing  function  is  sinusoidal , 


i.e. , 


f  (t) 


Ae 


i'jit 


r 


(205) 


let  the  particular  solution  be 

Xp(t)  =  Ke^^^^ 

Substituting  Equation  206  into  205,  we  obtain. 


K  = 


m(u)2  - 


+  2i^w^w) 


A 

Z  (o))' 


=  AH(co) 


(206) 

(207) 


where  Z(w)  is  called  the  "mechanical  impedance"  of  the  sys¬ 
tem,  and  H(a))  is  called  the  "frequency  response  function" 
or  the  "transfer  function"  of  the  system. 


When  the  forcing  function  is  impulsive,  i.e., 
f  (t)  =  C  6  (t)  , 

the  particular  solution  becomes , 
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x^(t)  =  C  h(t)  =  - ^  exp[(-c+i  Vl-?^ 

for  t  >  0  (208) 

The  function  h(t)  in  Equation  208  is  called  the  "impulse  re¬ 
sponse  function." 

;  I  :  ’ 

It  is  to  be  noted  here  that  the  response  of  a  stable 
linear  system  to  an  arbitrary  excitation  does  not  depend  a 
great  deal  on  the  initial  conditions  after  a  fairly  long 
time.  This  is  known  as  the  steady-stat^  condition,  during 
which  the  following  relations  exist  with  the  use  of  the 
principle  of  superposition: 

CO 

x^  (t)  =  /  H(a))  f(a))  exp(ib)t)  d 

—00 

t 

=  /  f(T)  h(t  -  T)  dT  (209) 

o 

where  ' 

00  ! 
f(w)  -  ^  !  f(t)  exp(i«t)dt,  and  (210) 

o 

00  .00 

f(t)  =  /  f^(w)  exp(ia)t)da)  =  /  f(T)h(t  -  T)dT  (2il) 

-00  o  , 

t 

The  integral  /  f(T)h(t-T)dT  is  called  the  "Duhamel  inte- 
o 

gral."  From  Equation  209,  the  following  convolution  theo¬ 
rem  of  Fourier  transformation  can  be  obtained: 

1  ” 

h(u)  -  JtF  ^  H(u)  exp  (iu)u)da),  '  (212) 

—00  • 

and 

00 

H(a))  =  /  h(t)  exp(-ia)t)dt  (213) 

—  00 

When  the  excitation  is  represented  by:  a  random  process 
F(t),  the  response  also  becomes  a  random  process  X(t). 
Correspondingly,  Equation  205  becomes 
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1 


(214) 


X(t)  +  2?WQX(t)  +  0)^  X(t)  =  i  F(t) 

assuming  that  (a)  the  derivatives  of  X(t)  exist  in  mean- 
square  sense  (131)  and  (b)  the  random  excitation  begins  at 
t  =  0. 

Because  the  principle  of  superposition  holds  for  a 
linear  system,  we  have 


X(t)  =  /  F(t)  h(t  -  T)dT 


(215) 


Also  because  integration  is  a  linear  operator,  we  obtain. 


E[X{t)3  =  /  EIF(t)1  h(t  -  T)dT 


(216) 


h  4 


4* ( b ,  t2 )  /  /  ( "^ ^  /  "^2 ^  ^  ^ ^ ^  2 ^  1^"^ 2 


o  o 


If  F(t)  is  stationary  in  the  wide  sense,  i.e., 
E[F(t)]  *  jip,  =  constant 


(217) 


(218) 


E[F(t)F(t  +  T)1  =  RFF(b) 

it  can  be  shown  that 

E[X(t)]  =  “  constant 


(219) 


(220) 


ElX(t)X(t  +  T)]  =  Rxx^"^^ 


(221) 


Equations  220  and  221  indicate  that  the  response  of  a  linear 
system  to  weakly  stationary  excitations  is  also  weakly  sta¬ 
tionary.  Furthermore,  the  spectral  density  of  the  response 
process  is  equal  to  the  product  of  the  spectral  density  of 
the  forcing  process  and  the  "transmittancy  function”  or  the 

2 

"system  function"  |h(w)|  ,  i.e.. 


(222) 


For  the  given  system  subjected  to  a  white  noise  excitation 
with  $pp((»))  =  Wq,  the  mean-square  response  is  found  to  be 
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as  follows: 


■9  "5  '^*^0 

ElX^(t)]  =  RxX^°^  "  I  |h(u)|  do)  =  - l-j  (223) 

—00  2^a)^in 

More  advanced  topics  in  random  vibration  are  summarized  and 
discussed  by  Lin  (131) . 

Rosenblueth  and  Bustamante  (133)  found  a  distribution 
function  for  the  structural  response  to  an  idealized  earth¬ 
quake  excitation,  which  was  represented  by  a  series  of  im¬ 
pulses  random  in  time.  Bogdanoff,  Goldberg,  and  Schiff 
(154)  studied  the  effect  of  ground  transmission  time  on  the 
response  of  long  structures.  Shinozuka,  Hakuno,  and  Ita- 
gaki  (155)  evaluated  the  response  of  a  multi-story  frame 
structure  to  nonstationary  random  excitation.  A  general  dy¬ 
namic  analysis  for  both  stationary  and  nonstationary  re¬ 
sponse  of  linear  structures  was  proposed  by  Shinozuka  and 
Yang  (156) .  A  numerical  method  for  the  analysis  of  complex 
structures  subjected  to  nonstationary  random  excitation  was 
also  made  available  by  Yang  and  Shinozuka  (157) .  Recently, 
the  evolutionary  power  of  random  processes  and  its  applica¬ 
tion  to  earthquake  engineering  was  presented  by  Shinozuka 
(158)  and  Brant  and  Shinozuka  (159) .  The  application  of 
random  vibration  theory  to  earthquake  engineering  problems 
was  also  contributed  by  Bogdanoff,  Goldberg,  et  al.  (22,23) 
and  Penzien  and  Liu  (17) .  Moreover,  the  problem  of  random 
vibration  of  nonlinear  systems  has  been  studied  by  many  in¬ 
vestigators  (160-167) . 


Structural  Safety  in  Earthquake  Engineering 

The  basic  problem  of  structural  safety  can  be  described 
in  the  following  manner.  Consider  a  structural  component. 
Usually,  the  resistance  is  considered  as  a  function  of  the 
applied  force,  and  failure  is  assumed  to  occur  whenever  the 
force,  S(t),  exceeds  the  resistance  R(t).  In  general,  both 
force  and  resistance  are  random  in  nature.  Therefore,  the 
"life"  of  this  structural  component  can  be  represented  by 
random  variable  T.  The  reliability  function  of  this  struc¬ 
tural  component,  L,p(t),  is  defined  as  the  probability  of 

the  event  that  its  useful  life  will  be  at  least  t,  i.e., 

L^(t)  =  P(T  >  t) 

=  P(R(t)  >  S(t),  0  <  t  <  t)  (224) 

Because  of  wear,  cumulative  damage,  and  increasing  chance  of 
encoxontering  larger  loads  with  increasing  time,  the  relia- 
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bility  of  a  structure  is  usually  a  monotonically  decreasing 
function. 

By  definition,  the  distribution  function  for  the  life, 

T  is  then, 

F^(t)  =  P(T  <  t)  =  1  -  L^(t)  (225) 

For  a  special  case  where  R(t)  =  r  =  constant,  and  S (t)  is  a 
known  random  function,  the  problem  of  finding  F^(t)  belongs 

to  the  "first  passage"  (or  "barrier"  or  "crossing")  problem, 
which  will  be  discussed  later.  Another  special  case  is  that 
when  R(t)  =  R,  S (t)  =  S,  both  R  and  S  are  random  variables 
with  known  distributions.  Freudenthal  first  proposed  the 
basic  concepts  of  structural  safety  more  than  twenty  years 
ago  (104).  Later,  Asplund  (168),  Brown  (169),  and  Freuden¬ 
thal  and  Shinozuka  (170)  derived  the  expression  for  the 
probeibility  of  failure.  Shinozuka  (171)  then  found  the 
probability  of  survival  under  n  applications  of  load.  A 
comprehensive  review  of  the  classical  studies  of  structural 
reliability  was  made  by  Freudenthal,  Garrelts,  and  Shinozuka 
in  1964  (105).  In  1967,  Cornell  (172)  gave  reliability 
bounds  for  cases  where  the  loads  and  the  resistance  are  not 
statistically  independent.  Amin  and  Ang  (173)  established 
a  monotonic  property  of  the  hazard  function  of  structures, 
which  is  subjected  to  a  sequence  of  random  loads.  Recent 
studies  of  structural  reliability  include  the  use  of  an  un¬ 
certainty  factor  (106,107),  redtindant  systems  (174,175), 
statically  indeterminate  structures  (176,177),  and  plastic 
collapse  of  structural  frames  (178). 

In  addition  to  its  use  as  a  common  basis  for  comparing 
the  relative  safety  among  alternative  designs,  the  proba¬ 
bility  of  failure  can  be  used  as  a  constraint  in  the  opti¬ 
mum  design  of  structures.  Hilton  and  Feign  (179)  used  a 
Lagrange  multiplier  technique  to  find  the  minimum  weight  of 
a  simple  structure  with  statistically  distributed  force  and 
resistance.  This  optimization  technique  was  later  refined 
by  Kalaba  (180).  In  1964,  the  simplex  linear  programming 
method  was  used  to  solve  the  optimization  problem  following 
the  development  of  a  linear  relationship  between  the  weight 
and  the  probability  of  failure  (181).  Then  Moses  and 
Kinser  (182)  developed  computational  techniques  for  the  de¬ 
termination  of  failure  probabilities  of  elastic  frames  and 
trusses  under  multiple  loading  conditions,  which  were  used 
in  the  optimum  structural  design.  More  recently,  the  ef¬ 
fect  of  proof-load  testing  to  reliability-based  structural 
optimization  was  studied  by  Shinozuka,  et  al.  (183,184). 

In  earthquake  engineering,  structural  failures  can  re- 
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suit  from  excessive  deformation,  dynamic  instability,  or  fa¬ 
tigue  deunage.  Collapse  of  yielding  structures  during  earth¬ 
quakes  was  studied  by  Jennings  and  Husid  (185) .  The  failure 
probability  of  seismic  structures  was  treated  as  a  barrier 
problem  by  Shinozuka,  et  al.  (155,  186-193)  and  others  (22, 
23,194,195) . 

The  deterministic  problem  of  dynamic  stability  has  been 
studied  by  many  investigators  (196-202) .  When  the  axial  ex¬ 
citation  of  the  column  is  represented  by  a  stationary  Gauss¬ 
ian  process,  various  stability  conditions  were  obtained  by 
Ariaratnam  (203) ,  Kozin  (204) ,  and  Caughey  and  Gray  (205) . 
Later,  Caughey  and  Dickerson  (206)  studied  the  stability  of 
systems  subjected  to  a  narrow-band  random  excitation.  Re¬ 
sults  from  a  recent  study  showed  that,  for  a  concentrically 
loaded  column  without  sidesway,  dynamic  instability  due  to 
the  vertical  component  of  the  earthquake  excitation  is  not 
likely  to  be  a  problem  (207) .  However,  significant  lateral 
motion  of  the  structure  occurs  and  thus  may  cause  the  col¬ 
umns  to  become  unstable  during  an  actual  earthquake. 

To  improve  the  design  of  earthquake-resistant  struc¬ 
tures,  it  is  desirable  to  know  more  about  the  failure  mech¬ 
anisms  of  these  structures.  Although  it  is  possible  to  have 
low-cycle  fatigue  failure  in  seismic  structures  (18,208), 
such  effects  have  not  been  studied  statistically.  To  date, 
the  studies  of  random  fatigue  consist  of  cases  with  station¬ 
ary  loading  processes,  which  were  summarized  recently  by 
Tang  and  Yao  (102).  Miles  (209)  was  the  first  to  compute 
the  mean  value  of  cumulative  fatigue  damage  for  a  narrow- 
band  random  process.  An  approximation  was  suggested  by 
Powell  (210)  for  the  case  where  the  stress  did  not  follow  a 
narrow-band  process.  The  variance  of  the  cumulative  damage 
resulting  from  narrow-band  random  stresses  were  computed  by 
Crandall,  Mark  and  Khabbaz  (211) ,  Bendat  (212) ,  and  Shino- 
zuka  (213).  Parzen  (214)  and  Freudenthal  and  Shinozuka 
(215)  applied  the  renewal  theory  to  the  study  of  random  fa¬ 
tigue.  Sweet  and  Kozin  (216)  incorporated  physical  proper¬ 
ties  into  the  renewal  model ,  which  was  compared  to  experi¬ 
mental  results.  Rice  (217)  developed  an  approximate  density 
function  of  the  rise  and  fall  statistics ,  which  was  used  for 
the  prediction  of  fatigue  crack  propagation  by  Rice,  Beer, 
and  Paris  (218).  A  comprehensive  review  of  experimental 
studies  in  the  field  of  random  fatigue  was  made  by  Swanson 
(219) . 


Decision  Analysis  in  Earthquake  Engineering 
In  engineering  problems,  the  engineer  must  decide  on 
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the  action  to  take  among  a  space  of  possible  acts  A  =  {a}. 
Also,  it  is  usually  assumed  that  there  exists  a  space  of 
possible  "states  of  nature  (the  way  things  really  are)" 

0  *  {6}.  In  order  to  obtain  information  about  the  state  of 
nature,  engineers  can  conduct  certain  experiments  from  a 
family  of  possible  experiments  E  =  {e}.  The  family  of  ex¬ 
periments  E  is  said  to  include  the  "dummy"  experiment  which 
refers  to  a  decision-making  without  experimentation.  The 
collection  of  all  possible  outcomes  for  various  experiments 
is  denoted  by  Q  =  {u}.  For  each  combination  (a,6,e,a))e 
AXOXEXfl,  it  is  assumed  that  the  engineer  can  assign  a  con¬ 
sequence  ceC. 

As  an  excimple,  the  action  to  be  decided  upon  might  be 
whether  to  accept  or  reject  a  lot  of  reinforcing  steel  bars. 
The  possible  states  of  nature  might  be  (a)  that  all  bars 
meet  the  specifications,  (b)  that  majority  meet  the  speci¬ 
fication,  (c)  that  only  few  bars  meet  the  specification, 
and  (d)  none  meets  the  specification.  To  obtain  some  in¬ 
formation  concerning  the  state  of  nature,  the  engineer  can 
perform  certain  experiments  on  a  selected  number  of  these 
bars.  The  result  of  these  experiments  can  be  either  satis¬ 
factory  or  unsatisfactory.  For  each  possible  combination 
of  (a,6,e,u),  the  consequence  can  range  from  excessive 
cost  to  defective  structure  as  an  end  product.  The  deci¬ 
sion  analysis  can  be  used  to  guide  a  decision  maker  such 
as  the  engineer  in  his  choosing  an  action  under  uncer¬ 
tainty. 

A  comprehensive  introduction  to  decision  analysis  with 
a  minimum  of  mathematical  demands  is  given  by  H.  Raiffa 

(220)  .  Other  introductory  references  include  G.  Hadley 

(221)  ,  H.  Raiffa  and  R.  Schlaifer  (222)  ,  Tribus  (99)  ,  and 
H.  Chernoff  and  L.  E.  Moses  (223) .  Recently,  many  appli¬ 
cations  of  the  decision  theory  to  structural  engineering 
have  been  proposed  by  several  authors  (108,109,  224-228). 
Also,  the  application  of  decision  analysis  and  probabil¬ 
istic  concepts  to  design  codes  was  explored  by  several  in¬ 
vestigators  (229-232) . 

In  1763,  Reverend  Thomas  Bayes  suggested  the  use  of 
Bayes'  Theorem  in  combining  the  probability  based  on 
hunches  with  that  based  on  relative  frequencies.  Gener¬ 
ally,  the  Bayesian,  or  subjective,  approach  to  probabil¬ 
istic  methods  incorporates  or  introduces  intuitive  judge¬ 
ments  and  feelings  into  the  formal  analysis  of  a  decision 
problem.  The  resulting  choice  of  an  action  is  said  to  be 
consistent  with  the  decision  maker's  preference  for  vari¬ 
ous  consequences  as  well  as  his  judgement  about  the  un¬ 
certainties  involved  in  the  problem. 
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To  begin  the  decision  analysis,  it  is  desirable  to  con¬ 
struct  a  "decision-flow  diagr2un"  or  "tree."  The  engineer 
has  a  choice  at  the  very  beginning:  to  be  involved  or  not  to 
be  involved  (e.g.,  he  can  quit  the  job  in  case  he  chooses 
not  to  be  involved) .  If  he  decides  to  be  involved  in  the 
decision-making  process,  he  must  decide  then  whether  to  ob¬ 
tain  more  information  through  experimentation  (eQ,e^,...). 

For  each  of  these  experiments,  there  are  possible  outcomes 

. . . )  »  which  depend  on  chance.  At  this  point,  the  en¬ 
gineer  Coin  choose  to  take  the  actions  (aj^,a2, . . .)  .  Then,  by 
chance,  the  state  of  nature  may  be  (6^, @2 »•••).  A  typical 

decision  flow  diagram  is  given  in  Figure  18.  Following 
Raiffa  (220)  ,  a  "decision  fork"  is  denoted  by  a  small  square 
and  a  "chance  fork"  by  a  small  circle.  Note  that  there 
exists  a  "consequence"  c  for  each  possible  "path"  in  the 
tree. 


At  each  chance  fork,  it  is  important  to  know  the  prob¬ 
ability  that  any  one  of  the  alternative  branches  will  be 
taken  by  chance.  In  other  words,  we  need  to  know  the  prob¬ 
abilities  of  experimental  results,  i.e.,  P(u)^),  and  the  con¬ 
ditional  probabilities  of  states  of  nature  given  that  a 
given  experimental  outcome  has  occurred,  i.e.,  P(0j|a)^). 

To  compute  these  probabilities  P(0j|a)^),  (from  quantities 
P(9j)/  P(a)^|0j),  the  following  Bayes'  Theorem  can  be  used. 


P(9j 


P(uj^|  9^)P(0^) 

^  P(u^l0j)P(0j) 

j 


(226) 


■a 


I 


i 

I 


Figure  18. 


Typical  Decision-Flow  Diagreun 
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•  •  • 


For  quantitatively  representing  consequences,  it  is 
first  necessary  to  define  a  preference  relation  "  "  over  the 
set  of  consequences  C  =  c  such  that  (a)  if  c^jC^eC,  then 

either  c^J-c^,  c^^c^,  or  both  and  (b)  if  c^fc^,  and  c^^Cj^, 
then  c^fCj^.  The  decision  maker  can  then  express  his  pref¬ 
erence  for  consequences  by  a  real- valued  "utility  function" 
u ( • )  such  that 

Ci^Cj  if  and  only  if  u(c^)  >  u(Cj) 

The  choice  of  experiment  e*  and  action  a*  should  be  made 
such  that 

Epio  tvi(e*»w,a*,0)  ]  =  max  EQoIu{e,(D,a,e)  ] 

aeA  (227) 

eeE 


APPENDIX  C 

SEISMIC  BEHAVIOR  OF  CONCRETE  AND  METAL  STRUCTURES 

General 

The  behavior  of  steel,  reinforced  concrete,  prestressed 
concrete,  and  composite  structures  has  been  tested  by  the 
recent  strong  motion  earthquakes  of  Skopje,  Yugoslavia,  in 
1963;  Niigata,  Japan,  in  1964;  Alaska,  U.S.A. ,  in  1964;  and 
Caracas,  Venezuela,  in  1967.  Especially  the  earthquakes  of 
Alaska  and  Caracas  have  revealed  valuable  information  re¬ 
garding  the  behavior  of  different  materials  as  well  as  the 
overall  behavior  of  multistory  buildings.  Recently  it  has 
been  customary  for  tecuns  of  structural  engineeers  to  go  to 
the  site  of  ein  earthquake  immediately  after  its  occurrence 
and  to  study  the  various  fonns  of  failure  and  behavior  of 
different  materials.  The  reports  of  these  engineers  are 
therefore  very  useful  for  locating  areas  of  further  re¬ 
search  . 

First  of  all,  these  damage  reports  of  the  four  afore¬ 
mentioned  earthquakes  are  examined  to  help  locate  the  prob¬ 
lem  areas.  Then  the  research  which  has  been  done  concern¬ 
ing  these  problem  areas  will  be  cited  and  discussed. 

Damage  Reports  from  Niigata  Earthquake 

1.  Steel  Building  Structures  (Minou  Makina,  Building  Re¬ 
search  Institute)  (233) 

Most  of  the  steel  structures  in  the  city  took  the 
earthquake  quite  well.  Where  there  were  fissures  and  de¬ 
pressions  and  if  the  subsoil  conditions  were  bad  there  was 
considerable  damage;  this  was  the  case  for  the  Niigata 
Timber  Manufacturing  Center.  There  were  occasional  failures 
at  beam-column  connections  (Plywood  Factory  in  Niigata, 
Pipeline  Station  for  Petroleum  Gas  in  Tsukiji  Village) , 
gusset  plates  (Agricultural  Warehouse  of  Niigata  Economic 
Union  in  Niigata) ,  and  fracture  of  tensile  members  (a  fac¬ 
tory  building  in  Nikaho  Town) .  None  of  these  buildings  was 
multistory  and  poorness  of  subsoil  conditions  was  shown  to 
be  the  primary  reason  for  most  failures.  A  boiler  house  in 
Niigata  which  is  a  steel  gabled  frame  house  (35  m.  high) 
and  which  sits  on  a  mat  foundation  with  8  m.  concrete  piles 
was  structurally  undamaged  even  in  a  destructive  area  of 
subsoil. 

2.  R/C  Constructions  (Koichiro  Ogura,  Meiji  University) 

(234) 

Out  of  a  total  number  of  340  normal  reinforced  concrete 
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buildings,  236  had  some  damage  and  104  had  no  damage.  Again 
out  of  77  box-type  reinforced  concrete  buildings,  only  17 
experienced  some  damage  while  60  experienced  none. 


Table  4  summarizes  the  influence  of  settlement  or  slant 
on  the  amount  of  cracking  of  concrete  in  all  the  damaged  re¬ 
inforced  concrete  structures. 


Table  4 

Effect  of  Settlement  on  Structural  Damage  (234) 


Notation 

Settlement 

Cracks  on 

Number  of 

and  SleUit 

the  Body 

Buildings 

Building  had 

Big  cracks  oc¬ 

shown  sink¬ 

curred  on  the 

A 

ing,  unequal 

body  over  the 

30 

settlement, 
or  slant  due 

whole  building. 

to  a  dis¬ 

Small  cracks  oc¬ 

B 

placement  or 

curred,  or  de¬ 

deformation 

flection  at  all 

46 

of  the 

the  expcmsion 

ground. 

joints  occurred. 

Cracks  scarcely 

Damage 

C 

occurred  or  no 
cracks  were  found. 

140  353 

Building 

Big  cracks  oc¬ 

D 

had  shown 

curred  (likely 

12 

no  sinking. 

due  to  vibration) . 

unequal 

settlement 

Small  cracks  oc¬ 

or  slant. 

curred,  cause  of 

25 

Cl 

which  could  not 
be  seen. 

No 

No  or 

No  cracks  or 

164 

damage 

c 

slight 

settlement 

slight  cracks. 

or  slant.* 

*Slight  slant  means  the  slope  is  less  than  one-hundredth 
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Reinforced  Concrete  Structures,  Syuzo  Takada  (Chiba 
University)  (235) 

From  Table  4  the  following  information  is  obtained. 

Nine  such  (reinforced  concrete)  buildings  were  observed  and 
four  different  modes  of  failure  were  considered:  (a)  damage 
to  body,  (b)  damage  to  finish,  (c)  overturning,  and  (d)  un¬ 
equal  sinking.  These  buildings  are  of  moderate  height  (6-11 
stories)  and  none  of  them  experienced  any  "damage  to  body." 
Only  one  was  observed  to  have  "damage  to  finish."  However, 
three  of  them  experienced  both  "overturning"  smd  "unequal 
sinking"  (same  three) . 

4.  Building  Damage  and  Soil  Condition,  Yorihito  Ohsaki 

(Building  Research  Institute)  (236) 

After  a  big  fire  in  1955  in  Niigata  City,  reinforced 
concrete  construction  expanded  greatly  and  at  the  time  of 
the  1964  earthquake  there  were  a  total  of  1,530  reinforced 
concrete  structures.  Of  these,  340  were  more  or  less  dam¬ 
aged  by  the  earthquake.  The  following  is  an  excerpt  from 
the  report  of  Y.  Ohsaki: 

"In  every  case  of  earthquakes  in  the  past,  wooden 
buildings  used  to  suffer  heavy  damage  and  reinforced  con¬ 
crete  structures  remained  safe.  However,  during  the  earth¬ 
quake  of  Niigata  last  year,  the  situation  was  completely 
reversed,  a  fact  which  may  be  regarded  as  the  most  charac- 
terstic  featvre  of  this  earthquake.  Considering  this  situ¬ 
ation,  it  sec.ns  appropriate  to  focus  attention  mainly  upon 
the  deimage  to  reinforced  concrete  buildings." 

The  Venezuela  Earth<^uaike  of  1967 
(Magnitude  6.5  on  Richter  Scale) 

Reinforced  concrete  construction  is  by  far  the  most 
widespread  type  of  construction  in  Venezuela.  In  the  Los 
Palos  Grandes  district  where  four  buildings  collapsed,  all 
the  multistory  construction  is  reinforced  concrete  with  the 
exception  of  the  two  30-story  steel  frame  towers  of  the 
Simon  Bolivar  Center.  Therefore  the  discussion  will  be  fo¬ 
cused  on  the  performance  of  reinforced  concrete  structures 
in  Caracas  during  the  1967  earthquake. 

The  concrete  used  in  construction  was  of  very  good 
quality  and  varies  in  strength  from  2400  psi  to  4300  psi. 

The  ordinary  reinforcing  steel  has  a  minimum  yield  strength 
of  34000  psi  which  corresponds  to  a  working  strength  of 
about  17000  psi.  A  type  of  high  strength  reinforcing  steel 
(Heliacero)  which  is  commonly  used  has  a  minimum  yield 
strength  of  56000  psi  and  a  working  strength  of  about  28000 
psi.  This  higher  strength  is  obtained  by  cold  working  the 
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steel  and  alters  the  stress-strain  diagram  by  eliminating 
the  flat  plateau  between  the  yield  point  and  the  start  of 
the  strain  hardening  region. 

The  typical  construction  of  multistory  buildings  in 
Caracas  consists  of  reinforced  concrete  frames  with  hollow 
tiles  as  fillers.  Most  of  the  design  practices  would  also 
correspond  to  the  provisions  of  the  U.S.  building  codes. 

One  exception  may  be  the  absence  of  top  reinforcing  in  the 
center  of  joists  and  slcib  steel  parallel  to  the  joist. 

This  would  be  considered  as  normal  in  areas  where  the  seis¬ 
mic  risk  is  very  small  but  not  in  relatively  severe  seismic 
areas  such  as  California  in  the  United  States.  One  major 
problem  which  was  completely  overlooked  in  design  was  the 
interaction  of  the  reinforced  concrete  frame  with  the  filler 
wall. 


Of  the  157  multistory  (from  4-22  stories,  most  between 
10-20  stories)  buildings  in  Los  Palos  Grandes  district  of 
Caracas,  Venezuela,  four  of  them  (10-12  stories)  completely 
collasped,  31  suffered  major  and  15  suffered  minor  damages 
(after  subjective  evaluation  of  Hanson  and  Degenkolb) . 

The  absence  of  the  top  reinforcing  steel  in  the  joists 
may  have  been  the  reason  for  the  failure  of  three  of  the 
collapsed  buildings.  This,  coupled  with  deviations  from 
original  plans  and  calculations  led  to  the  collapse  of 
four  buildings. 

Causes  for  major  structural  damage  can  be  classified 
in  three  groups .  The  first  one  is  due  to  the  unfavorable 
interaction  between  the  reinforced  concrete  frame  and  the 
filler  hollow  tile.  This  localized  the  damage  to  the  "soft" 
stories  where  no  filler  was  used  (especially  first  stories 
which  were  used  for  parking) .  However,  the  frames  in  these 
"soft"  stories  were  not  designed  to  carry  the  larger  forces 
and  therefore  failed.  The  second  type  of  failure  was  due 
to  the  vinexpectedly  high  overturning  moments.  These  moments 
caused  compression  failures  in  the  columns  and  in  some  cases 
ruptured  the  less  ductile  Heliacero  reinforcing  steel  in 
tension.  Last,  but  not  least,  is  the  failure  due  to  poor 
detailing  of  connections.  This  is  one  area  which  is  very 
importamt  and  common  to  both  reinforced  concrete  and  steel 
structures  in  earthquake  resistant  design. 

The  above  discussion  was  based  on  a  report  by  Hanson 
and  Degenkolb  for  the  American  Iron  and  Steel  Institute 
(237) .  The  PCA  report  (238)  also  finds  the  Caracas  codes 
and  design  procedures  comparaible  to  U.S.  practices  for  a- 
seismic  design.  It  is  concluded  that  "...most  of  the  dis- 
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tress  resulted  not  from  inadequate  codes,  but  from  a  lack  of 
understanding  of  the  dynamic  behavior  of  structures." 

Some  statistics  are  given  about  the  number  of  "unsafe” 
and  collapsed  buildings.  "Out  of  7000  buildings  higher  them 
four  stories,  approximately  180  buildings  were  damaged  and 
listed  as  temporarily  unsafe.  Of  these  180,  approximately 
40  had  structural  distress.  Five  out  of  the  1000  buildings 
higher  than  10  stories  collapsed." 

As  the  AISI  report  also  points  out,  most  of  the  high- 
rise  buildings  in  Caracas  were  reinforced  concrete  con¬ 
struction.  One  of  the  interesting  and  plausible  reasons  of 
failure  of  these  reinforced  concrete  frame  buildings  is  dis¬ 
cussed  in  the  PCA  report.  "High-rise  buildings  in  Caracas 
are  typically  reinforced  concrete  frautie  construction  with¬ 
out  shear  walls.  Although  the  buildings  were  designed  as 
frames,  the  interior  partitions  and  exterior  walls  are  con¬ 
structed  with  brittle  clay  tile  of  low  strength  but  high 
rigidity.  These  walls  give  the  buildings  the  stiffness  and 
response  of  shear  wall  buildings  and  influenced  the  distri¬ 
bution  of  earthquake  forces  within  the  structure." 

This  adverse  effect  of  frame-partition  interaction  was 
also  pointed  out  in  the  AISI  report  and  shows  how  different 
structures  can  behave  under  dynamic  loading. 

Skopje  Earthquake  of  1963  (239) 

According  to  the  data  gathered  by  the  Statistical  In¬ 
stitute  of  the  Socialist  Republic  of  Macedonia  the  damage 
inflicted  upon  the  housing  due  to  the  earthquake  can  be 
stunmarized  in  Table  5: 


Table  5 

Earthquake  Damages  in  Skopje  (239) 


Degree  of  Damage 

Flats 

Damaged 

(%) 

Living  Space 
Affected 
(%) 

Population 

Affected 

(%) 

Collapsed 

8.5 

7.05 

8.5 

Heavily  deimaged  (most  build¬ 
ings  have  to  be  demolished) 

33.6 

29.9 

36.4 

Moderately  damaged 

36.3 

39.9 

30.6 

Slightly  damaged 

19.0 

19.8 

20.3 

Undamaged 

2.6 

3.4 

4.2 
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Ambraseys  (240)  divided  the  buildings  in  Skopje  into 
four  categories:  (a)  old  adobe  construction  with  or  without 
timber  bracing;  (b)  load-bearing  brick  wall  construction 
supporting  reinforced  concrete,  or  wood  floors  supported 
partly  by  masonry  walls  and  partly  by  reinforced  concrete 
columns  and  beams;  (c)  reinforced  concrete  skeleton  build¬ 
ings  with  and  without  concrete  shear  walls;  (d)  there 
was  also  one  roof  structure  of  prestressed  concrete  and  a 
few  isolated  cases  where  use  of  prefabricated,  prestressed 
or  ordinary  elements  were  made. 

The  new  steel  mill  in  the  outskirts  of  Skopje  was  the 
only  steel  frame  construction. 

The  adobe  and  the  mixed  construction  performed  badly 
as  expected.  The  reinforced  concrete  skeleton  buildings 
were  far  better  and  only  two  collapsed.  The  taller  build¬ 
ings  (up  to  15  stories)  in  this  category  performed  much 
better  since  they  were  designed  for  wind  forces.  The  pre¬ 
stressed  construction  totally  collapsed  when  the  supporting 
columns  failed.  The  steel  mill  which  was  still  under  con¬ 
struction  had  only  minor  damage.  None  of  the  buildings  was 
designed  for  earthquake  forces. 

Alaskan  Earthquake  of  March  27,  1964  (241) 

The  damage  caused  during  this  earthquake  is  of  special 
importance  because  the  many  buildings  affected  by  the  earth¬ 
quake  were  designed  according  to  the  earthquake  provisions 
of  the  Uniform  Building  Code.  The  damage  was  caused  either 
by  the  seismic  sea  wave  (tsunami)  effects,  or  by  earth  vi¬ 
bration  effects.  Most  of  the  buildings  of  structural  inter¬ 
est  were  in  the  Anchorage  area,  which  experienced  primarily 
the  earth  vibration  effects  and  therefore  this  aspect  will 
be  considered  here.  Another  big  cause  of  damage  was  the 
landslides  which  took  place  because  of  the  earthquake;  a  de¬ 
tailed  study  of  the  soil  conditions  leading  to  these  land¬ 
slides  can  be  found  in  the  report  made  by  the  Shannon  and 
Wilson  firm  (242) .  This  earthquake  also  demonstrated  the 
correlation  between  the  epicentral  distance  to  the  struc¬ 
tures  mostly  affected  by  the  earthquake.  The  epicenter  of 
the  Alaskan  earthquake  was  cibout  75  miles  from  Anchorage 
and  therefore  the  dominant  periods  were  quite  long  (more 
than  0.5  sec.).  For  this  reason,  taller  and  more  flexictle 
structures  attracted  more  force,  and  this  fact  should  be 
kept  in  mind  when  the  damage  reports  are  studied.  The  fol¬ 
lowing  is  a  description  of  the  damage  on  the  metal  and  con¬ 
crete  structures  as  quoted  by  Steinbrugge  (241) . 


"All  metal  buildings  performed  excellently  as  may  be 
typified  by  the  usual  absence  of  damage  to  all  steel  gasoline 
service  stations. 

"Some  damage,  however,  was  noted  in  all  metal  Structures 
having  large  elevated  masses.  Equipment  shifted,  one  bin 
fell,  X-bracing  broke,  metal  skin  buckled,  and  steel  columns 
twisted  at  the  Chugach  steam  power  generating  plant  located 
in  the  Ship  Creek  section  of  Anchorage;  differential  settle¬ 
ments  indicated  that  foundation  problems  probably  accentuated 
this  damage.  Nearby,  the  Elmendorf  Air  Force  Base  had  its 
principal  damage  restricted  to  the  steel  connections  from  a 
large  elevated  pit  to  the  main  structure.  Overall  damage  to 
both  power  plants  was  estimated  to  be  slight,  and  overall 
equipment  deunage  probably  could  also  be  classified  as  being 
slight.  A  third  power  plant,  owned  by  the  City  of  Anchorage 
and  also  located  in  the  Ship  Creek  area,  had  apparently  neg¬ 
ligible  damage. 

"Hollow  concrete  block  was  a  common  construction  mater¬ 
ial  for  small  merccuitile  structures  as  well  as  for  industri¬ 
al  structures.  Roof  and  supported  floors  were  usually  wood. 
When  small  in  area,  when  one  story  high  and  when  not  located 
in  the  land  movement  areas ,  such  damage  as  sometimes  occurr¬ 
ed  was  usually  no  more  than  slight  to  moderate.  Some  para¬ 
pets  fell,  unanchored  roofs  punched  out  sections  of  the 
hollow  concrete  block  walls,  but  collapse  was  uncommon. 

Wire  webbing  for  reinforcement  was  usually  laid  in  selected 
horizontal  joints  of  the  hollow  concrete  block.  The  place¬ 
ment  of  vertical  steel  was  inconsistent.  This  reinforcement, 
in  general,  appeared  to  be  of  a  size  and  amount  that  would 
not  be  considered  to  be  fully  adequate  in  other  sections  of 
the  United  States  which  are  considered  seismically  very  ac¬ 
tive.  Workmanship  was  often  poor  where  collapsed  walls 
were  noted;  the  concrete  grout  did  not  fill  all  of  the  cells 
in  many  observed  instances.  i 

"Poured-in-place  reinforced  concrete  wall  construction 
performed  well  for  small  buildings.  Usually  the  roof  and 
supported  floor  materials  were  wood  for  the  small  buildings > 
although  not  always.  Instances  of  metal  deck  and  metal 
open  web  joist  roof  and  floor  systems  were  found.  The  per¬ 
formance  of  these  structures  was  generally  good  and  somewhat 
better  than  similar  size  structures  of  different  masonry 
materials  when  not  located  in  the  land  movement  areas." 

The  five-story  Penney  Building  which  was  entirely  of 
reinforced  concrete  (both  poured-in-place  and  precast) 
collapsed.  Poor  detailing  especially  of  the  poured-in-place 
to  precast  connections  was  a  major  cause  for  the  collapse. 
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Also,  large  torsional  forces  were  developed  in  the  upper 
stories  because  second  floor  and  above  was  structurally  a 
U-shaped  building. 

The  Four  Seasons  Apartment  Building  was  six  stories  and 
was  only  structurally  completed  (it  had  not  yet  been  finish¬ 
ed)  .  It  had  prestressed  concrete  (in  both  directions)  slabs 
for  floors  which  were  supported  by  steel  columns.  It  also 
had  two  cores  of  poured-in-place  reinforced  concrete.  These 
two  cores  failed  and  the  building  collapsed.  Faulty  splic¬ 
ing  and  shear  connections,  and  the  lack  of  grouting  of  the 
prestressing  cables,  were  pointed  out  to  be  the  major  causes 
of  failure.  The  building  was  very  closely  located  to  the 
main  graben  of  the  "L"  Street  landslide. 

Hill  Building  which  was  eight  stories  had  a  central 
core  of  poured-in-place  reinforced  concrete,  a  steel  frame, 
and  at  the  upper  stories  poured-in-place  reinforced  concrete 
one-way  slabs.  The  major  damage  was  concentrated  at  the 
central  core.  Lc±>oratory  tests  for  the  concrete  indicated 
excessive  organic  material  which  had  resulted  in  extremely 
low  strength  concrete. 

Steinbrugge  (241)  completes  his  report  with  the  follow¬ 
ing  conclusion:  "Precast  reinforced  concrete  performed 
rather  poorly  in  too  many  instances.  However,  failures  were 
almost  inevitably  associated  with  the  connections  to  the  pre¬ 
cast  elements.  It  would  appear  obvious  that  the  material 
was  usually  not  at  fault,  rather  a  lack  of  sound  engineering 
judgment  was  the  principal  cause.  There  was  no  evidence  to 
indicate  that  one  construction  material  was  superior  to  an¬ 
other  when  given  comparsdile  design  and  construction  atten¬ 
tion. 

After  presenting  these  damage  reports  of  the  four  earth¬ 
quakes,  we  shall  now  turn  to  the  research  done  in  the  areas 
of  steel,  reinforced  concrete,  prestressed  concrete,  and 
composite  structures  under  strong  motion  earthquakes.  A 
comparison  of  the  performance  of  different  materials  is  un¬ 
warranted.  Steinbrugge ' s  quotation  in  the  last  paragraph 
quite  accurately  sums  up  the  situation.  However,  one  point 
should  be  made  regarding  concrete  structures.  Concrete  is 
a  newcomer  to  the  field  of  earthquedce  engineering.  "Ductile 
concrete"  is  a  very  recent  concept  (50) .  Therefore  the 
majority  of  the  literature  concerns  reinforced  or  prestress¬ 
ed  concrete  behavior  under  earthquake  loading. 

One  of  the  biggest  problems  common  to  both  steel  and 
concrete  structures  concerns  detailing.  Be am- to- column 
joints  are  cne  of  the  major  problem  areas.  The  following 


101 


will  coii^rise  a  detailed  discussion  of  connections,  both  con¬ 
crete  and  steel,  and  also  a  representative  paper  will  be  pre¬ 
sented  for  each  material  emphasizing  both  the  merits  and  the 
shortcomings.  Other  research  will  be  cited  briefly. 

Reinforced  Concrete 

Splicing  is  one  of  the  important  details  of  beam-column 
connections.  At  the  University  of  Illinois  a  series  of  ten 
tests  were  carried  out  on  lapped  splices  in  beam-column 
connections  (243) .  The  main  objective  was  to  study  the  ef¬ 
fects  of  splice  length,  splice  location,  and  stirrup  spacing 
on  the  behavior  of  tension  splices  in  column  specimens.  The 
following  is  an  excerpt  from  the  conclusions  of  the  investi¬ 
gators  (243) : 

"1.  For  equal  splice  lengths,  splicing  away  from  the 
point  of  maximum  moment  increases  both  the  ultimate  deflec¬ 
tion  and  the  steel  stress  at  ultimate. 

"2.  Placing  the  through  bar  on  the  inside  of  the  cross 
section  resulted  in  a  lower  load  and  deflection  at  ultimate. 

"3.  A  closer  spacing  of  stirrups  over  the  splice  length 
resulted  in  greater  ductility  and  also  larger  maximiam  steel 
stresses  for  the  same  concrete  stress. 

"4.  The  ultimate  deflection  of  all  the  spliced  speci¬ 
mens,  except  C-10,  was  less  than  half  that  of  a  similar  un¬ 
spliced  column." 

Another  valuable  study  of  reinforced  concrete  beam- 
column  joints  was  carried  out  by  Hanson  and  Connor  (244)  . 

This  study  specifically  concerned  the  seismic  resistance  of 
the  joints.  The  authors  take  the  recommendations  of  Blume, 
Newmark,  and  Corning  (50)  as  a  basis  for  their  investigations 
and  summarize  them  as  follows: 

"Crucial  recommendations  for  earthquake  resistant  de¬ 
sign  are  the  following; 

"1.  Sufficient  treinsverse  and  shear  reinforcement  to 
provide  a  shear  strength  greater  than  the  flexural  strength. 

"2.  Limitations  on  the  amount  of  tensile  reinforcement 
or  required  use  of  compression  reinforcement,  to  ensure 
ductility  and  energy  absorbing  capacity. 

"3.  Confinement  of  the  concrete  by  hoops  or  spirals  of 
critical  sections  such  as  beam  column  connections,  to  in¬ 
crease  the  ductility  of  columns  under  combined  axial  load 
and  bending. 

"4.  Special  attention  to  details,  such  as  splices  in 
reinforcement  and  exclusion  of  planes  of  weakness  that  would 
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result  from  bending  or  terminating  all  bars  at  the  same  sec¬ 
tions  . " 


With  these  recommendations  in  mind,  the  investigators 
devised  a  test  to  determine  joint  reinforcement  required  to 
ensure  maintaining  ultimate  capacity  for  cast-in-place  beams 
euid  columns  subjected  to  multiple  reversals  of  loading  of 
major  earthquake  magnitude.  An  exterior  beam-column  connec¬ 
tion  was  selected  for  study  since  this  is  the  most  critical 
in  a  multistory  structure. 

The  following  conclusions  were  drawn  by  the  authors: 

"The  series  of  tests  demonstrate  that  properly  designed 
and  detailed  cast-in-place  reinforced  concrete  frames  can 
resist  moderate  earthquakes  without  damage  and  severe  earth¬ 
quakes  without  loss  of  strength.  Adequate  energy  absorption 
is  provided  by  ductility  of  reinforcing  steel.  Joints  con¬ 
necting  beeuns  and  columns  need  special  attention  in  design: 

"1.  Hoops  are  required  for  unconfined  (isolated)  beam- 
column  joints.  A  design  procedure  for  hoops  based  on  sup¬ 
plying  adequate  confinements  emd  shear  resistance  will  pro¬ 
vide  safe  designs. 

"2.  Corner  joints  with  beams  on  only  two  column  faces 
should  be  designed  as  unconfined  joints  requiring  hoops  un¬ 
til  tests  are  made  for  this  case. 

"3.  Hoops  are  not  required  for  exterior  joints  con¬ 
fined  on  at  least  three  sides  by  beams  or  spandrels  of  ap¬ 
proximately  equal  depth  suid  meeting  ACI  318-63  requirements 
for  concrete  strength  needed  to  transfer  the  column  load 
through  the  joint. 

"4.  The  cumulative  ductility  of  a  test  specimen  pro¬ 
vided  a  measure  of  the  ability  of  a  structure  to  withstand 
seismic  deformation.  Well  detailed  joints  sustained  high 
values  of  cumulative  ductility  while  they  maintained  their 
strength.  Omission  of  important  hoop  reinforcement  reduced 
amount  of  cumulative  ductility  which  could  be  sustained. 

"These  details  of  joint  reinforcement  are  not  limited 
to  seismic  design.  Good  design  practice  requires  that 
joints  be  designed  as  strong  as  the  adjoining  members. 

These  tests  indicate  that  for  the  rare  case  of  an  isolated 
joint,  a  design  for  the  usual  forces  of  wind,  dead,  and  live 
load  will  require  hoops." 

More  specialized  literature  includes  the  earthquake 
effects  on  reinforced  concrete  chimneys  (245,246),  walls, 
(247),  shells  and  plates  (248),  freunes  (249),  cUid  buildings 
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(250-252) .  The  dynamic  characteristics  of  plain  concrete 
(253,254),  and  the  effect  of  reinforcement  on  the  ductility 
of  concrete  members  (255) ,  are  also  relevant  topics  to  earth- 
qucdce  engineering. 

Since  the  topics  of  dynamic  torsion  and  dynamic  sta¬ 
bility  have  not  been  explored  extensively,  static  torsion 
(256,257)  and  axial  loads  (258)  must  serve  as  guidelines 
for  the  time  being. 


Prestressed  Concrete 

Before  the  Alaskan  earthquake  of  1964,  prestressed 
concrete  structures  had  not  been  tested  by  a  strong  motion 
earthquake.  Unfortunately,  during  this  earthqueJce  some  pre¬ 
stressed  concrete  structures  performed  very  badly.  It  was 
later  discovered  that  the  failures  were  mainly  due  to  faulty 
aseismic  design  and  foundation  failure,  rather  than  the  in¬ 
adequacy  of  the  material.  Nevertheless,  these  failures 
coupled  with  the  comparative  lack  of  knowledge  about  the 
seismic  behavior  of  prestressed  concrete,  made  many  engineers 
uneasy  about  this  material  for  use  in  earthquake  resistant 
structures.  Most  of  the  aseismic  design  codes  around  the 
world  do  not  make  explicit  reference  to  prestressed  concrete. 
An  exception  for  this  is  the  French  code  entitled  "Rdgles 
Paraseismiques  P.S.  64,"  which  explicitly  admits  prestressed 
concrete  without  providing  any  special  limitation.  There 
are  misconceptions  and  conflicting  ideas  about  the  seismic 
performance  of  prestressed  concrete,  and  extensive  experi¬ 
mental  data  is  urgently  needed  on  this  subject. 

In  a  paper  by  Despeyroux,  some  of  the  questions  which 
most  frequently  arise  regarding  the  seismic  behavior  of  pre¬ 
stressed  concrete,  are  discussed  (259). 

On  the  subject  of  resistcince  to  reversed  loadings, 
Despeyroux  points  out  that  as  long  as  the  reversed  loading 
does  not  exceed  the  yield  stresses,  fatigue  is  not  a  problem. 
For  reinforced  concrete  elements  the  resistance  to  reversed 
loading  remains  uncheinged  up  to  80%  of  the  ultimate  strength 
(260,261).  Therefore,  he  argues,  there  is  no  reason  not  to 
expect  the  same  performance  from  prestressed  members.  How¬ 
ever,  this  theory  has  not  yet  been  confirmed  by  experimental 
results. 

The  second  question  pertains  to  the  excessive  flexibil¬ 
ity  of  prestressed  concrete  members  compared  to  members  of 
the  same  bearing  capacity  made  with  other  materials.  This, 
in  turn,  leads  to  deimage  in  brittle  non-structural  elements. 
Despeyroux  first  points  out  the  advantages  of  flexibility 
in  aseismic  design.  The  major  adveintage  is,  of  course,  the 
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fact  that  flexible  structures  in  general,  attract  less 
earthquake  forces  than  more  rigid  structures.  Then  Despey- 
roux  recommends  using  stricter  code  provisions  against  drift 
and  hcunmering  to  prevent  excessive  unwanted  deformation. 

The  third  question  concerns  the  ductility  and  energy¬ 
absorbing  capacity.  Using  the  experimental  results  of  Caul¬ 
field  and  Patton  (260)  on  prestressed  concrete,  and  his  own 
theoretical  moment-curvature  relationship  for  reinforced  con¬ 
crete,  the  author  concludes  that  both  the  energy-absorbing 
capacities  (areas  xinder  moment-curvature  curves) ,  and  the 
ductilities  (ratios  of  the  ultimate  deflection  to  yield 
deflection)  of  prestressed  and  reinforced  concrete  are  of 
the  same  order  of  magnitude.  Below  is  the  graph  from  which 
he  draws  his  conclusions. 


Figure  19.  Moment  vs  Curvature 

This  conclusion  of  Despeyroux  is  also  experimentally 
confirmed  by  Guyon  (261)  and  Sutherland  (262). 
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As  for  the  joints  of  prestressed  concrete  men±)ers, 
Despeyroux  cautions  the  design  engineer  to  use  his  good 
judgment.  He  says  that  these  problems  are  well  known  and  re¬ 
fers  to  some  literature  on  the  subject  (263) . 

In  a  similarly  written  paper  Lin  (264)  concurs  with  just 
cd>out  every  point  that  Despeyroux  makes.  Other  useful  mater¬ 
ial  on  the  subject  includes  some  experimental  model  studies 
by  Nakano  (265)  and  determination  of  dynamic  response  of  pre¬ 
stressed  concrete  beams  by  Hamilton  (266) . 

Steel 

Most  engineers  view  steel  as  the  most  suitable  and  re¬ 
liable  material  for  multistory  earthquake  resistant  build¬ 
ings.  Both  the  strength  and  the  ductility  of  steel  struc¬ 
tures  tend  to  be  greater  than  similar  buildings  of  different 
material.  It  is  also  true  that  during  strong-motion  earth¬ 
quakes,  steel  buildings  suffered  relatively  little  damage. 
However,  one  very  important  point  in  steel  design  as  well  as 
design  of  buildings  with  any  other  material,  is  the  connec¬ 
tions.  Low-cycle  fatigue  and  brittle  failure  of  beam-column 
joints  of  steel  structures  may  have  occurred.  For  this  rea¬ 
son  an  experimental  work  by  Bouwkamp  (267)  on  girder- to- 
column  steel  connections  will  be  discussed  below. 

This  experiment  was  a  detailed  study  on  four  large  full- 
size  girder-to-column  connections.  Each  connection  was  com¬ 
posed  of  a  riveted,  36-inch  by  36-inch  built-up  column  with 
bolted  T-sections  amd  welded  plates  which  acted  as  moment 
resistant  supports  for  42-inch-deep  welded  girders.  The 
material  used  in  the  design  of  girders  and  columns  was  a 
combination  of  A7  and  A373  steel.  The  two-way  framing  sys¬ 
tem  without  any  interior  columns  made  the  length  of  the 
girders  93  feet.  There  were  two  objectives  of  this  project 
and  they  are  summarized  as  follows  in  the  report; 

"The  first  objective  was  to  determine  the  stiffness  of 
the  T-Sections  which  accomodate  the  flanges  of  the  93-foot 
long  floor  girders.  This  information  was  important  in  order 
to  evaluate  the  end-restraining  effect  (rotation)  of  the 
large  span  girders  and  to  increase  if  necessary  the  rigidity 
of  these  connections  in  order  to  reduce  the  floor  deflec¬ 
tions.  Furthermore,  an  acceptable  rigidity  would  favorably 
contribute  to  the  overall  lateral  stiffness  of  these  build¬ 
ings.  The  second  objective  was  to  determine  the  yield  and 
ultimate  strength  of  the  1-1/2-inch  and  2-inch  thick  welded 
plates  between  the  column  flanges  and  web.  These  plates  are 
30  inches  wide  and  H.  S.  bolted  to  the  flanges  of  the  43- 
inch  deep  girders  spanning  the  short  distance  between  the 
columns  in  the  exterior  frames.  Since  no  previous  informa- 
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tion  was  available  about  the  strength  of  these  welded  plates  • 
as  affected  by  residual  stresses  due  to  welding,  two  geomet¬ 
rically  identical  specimens  of  the  types  shown  in  Figures  2 
and  3  were  to  be  investigated,  one  non-stress  relieved  and 
one  stress  relieved.  Should  the  tests  indicate  the  need  for 
stress  relieving,  all  girder  to  column  connections  were  to 
be  treated  to  relieve  the  residual  stresses . " 


The  following  was  the  test  program: 

"1.  Specimen  1  (Type  I,  third  floor  connection,  non¬ 
stress  relieved) , 

2.  Specimen  2  (Type  I,  third  floor  connection,  stress 
relieved) , 

3.  Specimen  3  (Type  II,  fourteenth  floor  connection, 
non-stress  relieved) , 

4.  Specimen  4  (Type  II,  fourteenth  floor  connection, 
stress  relieved) . " 


The  yield  cind  ultimate  loads  for  the  welded  plates  were 
as  follows; 

"Specimen  1.  non-stress  relieved;  welded  plate 
30x2  in? 

Py=1800  k,  fy*30  ksi;  P^j^^=2805  k,  ksi 

Specimen  2.  stress  relieved;  welded  plate  30x2  in? 

Py=2200  k,  fy=37  ksi;  P^j^^=2750  k, 

Specimen  3.  non-stress  relieved,  welded  plate 
30x1-1/2  in? 

Py=1400  k,  fy=31  ksi;  P^j^^=2200  k,  ksi 

Specimen  4.  stress  relieved;  welded  plate  30x1-1/2  in? 

After  load  of  1600  k  was  reached  end  weld  was  cut 

P  ,  =1344  k,  f  h.^44  ksi 

ult  '  shear  ult 


"The  results  for  Py  indicated  that  the  stress  distri¬ 
bution  in  the  moment  plates  was  improved  by  stress  relieving. 
Similar  conclusions  were  derived  from  the  stress  distribu¬ 
tion  along  critical  sections  close  to  the  welds.  These 
welds  indicated  that  the  distribution  of  the  load  over  the 
several  welds  was  better  for  the  two  stress  relieved  connec¬ 
tions  than  for  the  non-stress  relieved  specimens." 


The  prototype  for  this  experiment  was  the  Health 
Sciences  Instruction  and  Research  Buildings  for  the  Univer¬ 
sity  of  California,  San  Francisco  Medical  Center,  as 


107 


designed  by  Ried  and  Tarics,  Architects  and  Engineers,  San 
Francisco.  This  complex  consists  of  two  buildings:  the  Main 
East  Building  which  is  15  stories  high  and  the  West  Building 
which  is  16  stories.  There  is  also  a  connecting  corridor 
between  the  two  buildings,  two  mechanical  service  towers, 
and  an  elevator  tower.  The  Main  East  Building  has  been  the 
site  for  some  more  dynamic  testing  done  by  the  University 
of  California,  Berkeley  (268) .  These  involve  frequency  re¬ 
sponse,  time  response,  mode  shapes,  energy  transfer,  and 
damping  characteristics. 

Some  of  the  important  conclusions  are  listed  below: 


1.  "The  steel  frame  of  the  service  tower  had  very  little  effect 
cn  the  c^amic  behavior  of  the  East  Building  in  the  sutmer  of  1964,  and 
the  East  Bioilding  was  tested  virtually  in  isolation.  However,  vrfien  the 
steel  frame  had  been  encased  in  concrete,  the  service  tower  did  eiffect 
the  behavior  of  the  East  Building;  in  fact,  both  buildings  then  formed 
a  new  structural  syston.  This  interconnected  system  had  much  larger 
danping  ccpacity  then  any  of  the  isolated  buildings." 


2.  "The  increase  of  damping  after  the  summer  of  1964 
was  caused  mainly  by  the  connections  between  buildings." 

3.  "A  standard  open  frame  model  was  found  to  reproduce 
accurately  the  dynaunical  behavior  of  the  East  Building. 

This  model  was  subjected  to  the  ground  acceleration  record 
of  the  El  Centro  earthquake.  It  was  found  that  the  amount 

of  damping  in  the  model  played  an  important  role  in  determin¬ 
ing  the  extent  to  which  yielding  occurred." 

Damping  of  this  building  is  discussed  in  more  detail  in 
reference  (269). 

Behavior  of  welded  connections  (270)  and  inelastic  re¬ 
sponse  of  steel  structures  (271-275)  also  have  application 
in  the  field  of  earthqueUce  engineering. 

More  experimental  research  is  necessary  to  determine 
the  properties  and  the  behavior  of  materials.  The  parameter 
that  is  used  in  the  equations  of  motion,  and  the  idealiza¬ 
tion  that  is  made  in  the  analysis,  depend  on  the  material 
properties  to  a  great  extent.  Also  as  was  mentioned  earlier, 
the  behavior  of  details  and  connections  of  structural  ele¬ 
ments  made  from  different  materials,  is  very  important. 

The  dynamic  effect  can  be  divided  into  two  effects  as 
applied  to  materials.  The  first  effect  is  that  of  the 


strain  rate  t.  As  the  strain  rate  increases,  the  yield 
stress  and  the  modulus  of  elasticity  (very  slightly)  in¬ 
crease.  However,  this  is  much  more  pronounced  in  blast  type 
loadings  and  the  increases  can  be  ignored  for  earthquake 
considerations.  The  second  one  is  the  effect  of  reversed 
action  or  what  may  be  called  the  Bauschinger  effect.  This 
results  in  deterioration  or  decrease  in  both  the  yield  stress 
and  the  modulus  of  elasticity.  This  is  a  problem  in  earth¬ 
quake-  type  loading  and  should  be  investigated  further.  Some 
literature  on  the  subject  (276-280)  is  presently  available. 

Idealization  of  Mechanical  Behavior 

No  matter  what  material  is  to  be  used,  there  must  be  an 
idealization  of  the  material,  the  cross  section,  the  member, 
and  the  structure.  Nonlinear  and  inelastic  idealizations 
proved  to  be  the  most  suitable  for  strong  motion  earthquake- 
resistcuit  structural  euialysis.  Therefore  these  topics  have 
received  considerable  attention  in  recent  years  (13,14,281- 
288) .  It  is  well  known  that  the  damage  due  to  repeated  de¬ 
formations  is  cumulative  and  could  lead  to  low-cycle  fatigue 
failures  (81,82).  Preliminary  results  indicate  that  it  is 
possible  to  have  low-cycle  fatigue  failures  of  seismic  struc¬ 
tures  (18)  . 

Below  is  a  summary  of  different  idealizations  and  their 
deviations  from  the  "true  behavior"  (289). 

The  following  relationships  must  be  studied: 


a 

vs. 

e 

for 

the 

material 

(stress 

vs.  strain) 

M 

vs. 

<t> 

for 

the 

cross  section 

(moment 

vs.  curvature) 

M 

vs. 

0 

for 

the 

member 

(moment 

vs.  rotation) 

H 

vs. 

6 

for 

the 

structure 

(load  V'' 

.  deflection) 
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Figure  20.  Stress-Strain  Relations 


There  are  two  problems  concerning  the  stress-strain  relation¬ 
ship. 

A.  Does  the  strain  rate  e  have  any  effect  on  a  vs.  e?  This 

effect  is  usually  considerable  in  blast  loading,  but  for 
earthquakes  the  increase  in  E  and  is  negligible. 

B.  What  is  the  effect  of  the  reversed  action  (the  Baushing- 
er  effect) ?  This  may  be  a  problem  in  earthquake  type  load¬ 
ing. 


1.  M  vs,  4>  of  a  cross  section 


Figure  21.  Moment  vs.  Rotations 
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Figure  21  (continued) 
Moment  vs.  Rotations 
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2.  M  vs.  6  for  the  member 


Figure  22.  Free-Body  Diagram  of  a  Beam 

(a)  First  order  theory  (neglect  change  in  geometry) . 

(b)  Second  order  theory  beam-column  problem. 

3.  H  vs.  6  for  the  structure 

(a)  First  order  theory 

(b)  Second  order  theory — P-A  effects 

All  the  combinations  of  assumptions  under  1,  2,  and  3 
are  plotted  on  Figure  23  (289) ,  which  also  includes  the  so- 
called  "true  behavior."  This  diagram  dramatically  illus¬ 
trates  the  differences  in  the  load-deflection  relations  of 
unbraced  frames  under  various  idealizations. 
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Figure  23 .  Behavior  of  Unbraced  Frames 
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APPENDIX  D 

AN  ILLUSTRATION  OF  THE  DIRECT  APPROACH 


Until  an  extensive  research  program  Is  carried  out,  It  Is  difficult 
to  find  a  numerical  example  demonstrating  completely  the  direct  approach 
In  the  seismic  design  of  building  structures.  Nevertheless,  It  Is  pos¬ 
sible  to  Illustrate  the  concepts  Involved  herein  with  a  simple  example, 
which  Is  described  In  this  appendix. 

Consider  the  seismic  design  of  a  slngle-degree-of-freedom  bilinear 
hardening  structural  frame.  The  seismic  behavior  of  such  a  structure  was 
studied  by  Yeh  and  Tao  (87),  who  found  that  the  maximum  displacement 
response  of  the  bilinear  structure  to  earthquake  ex''ltatlon  can  be  less 
tlun  the  maximum  displacement  of  a  linear  system  with  either  a  soft  or 
a  hard  stiffness.  Moreover,  this  bilinear  hardening  system  compares 
favorably  with  nonlinear  systems  undergoing  plastic  deformations  because, 
with  a  proper  design,  the  members  remain  elastic  to  avoid  possible  low- 
cycle  fatigue  damage  In  the  structure (18). 

For  this  example,  the  direct  approach  portion  of  the  design  tree 
can  be  constructed  as  shown  in  Figure  24.  The  alternative  solution  A2im 
to  the  excitation  question  Q2111  Is  chosen  to  be  the  1940  El  Centro 
earthquake  record.  The  digitized  plot  of  this  acceleration  time-history 
la  shown  along  with  a  response  function  In  Figure  25.  In  a  similar 
manner,  the  second  alternative  solution  A21112  Is  chosen  to  be  the  1952 
Taft  earthquake  record. 

The  alternative  solution  A22I  to  the  model  question  Q22  is  chosen 
to  be  a  discrete  model  as  shown  In  Figure  26.  The  diagonal  bracings  are 
purposely  made  slack  Initially.  As  the  deformation  reaches  a  predetermined 
value  A,  one  of  these  cables  will  become  tight  and  begin  to  contribute  to 
the  frame  stiffness.  This  structural  frame  can  be  represented  by  a 
mechanical  systan  with  bilinear  stiffness  as  shown  In  Figure  27.  The 
equation  of  motion  for  this  system  Is  given  as  follows: 

mu  +  cii  +  kj^u  +  sgn(u)(k2-kj^)  (lul-A)H(lul-A)  “  f(t)  (228) 

where  m  Is  the  mass  of  the  system 
c  Is  a  damping  constant 
kj^  and  k2  are  spring  constants  (k2>kj^) 
f(t)  Is  the  forcing  function 
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A-..,,:  1940  El 
Centro 
(Figure  25) 


Safety 


Figure  24.  Application  of  Design  Tree  In  the  Design  of  a 
Bilinear  Structural  Frame 
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Figure  25.  Acceleration  Record  of  1940  El  Centro  Earthquake 
and  a  Typical  Structural  Response 


Figure  26.  One-Degree-of-Freedom  Bilinear  Structural  System 
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u  is  a  variable  denoting  the  displacement  response 

A  Is  the  distance  at  which  the  displacement-force  relationship  Is 
piecewise  continuous. 


8gn(u)< 


!1,  If  u>0 
0,  If  u-0 
-l,lf  u<0 


H(x)  Is  the  Heaviside  step  function 
For  convenience  the  following  expressions  are  also  Introduced, 
.  58. 


u 


St  k. 


g  Is  a  gravitational  constant 
C 


2  mk. 


For  this  particular  case 


,  ^  m  and  ^  m  ^ 


were  arbitrarily  chosen  to  be  3.23 
rad/sec  and  3.76  rad/sec,  respectively,  and  c,  a  dimensionless  damping 
factor,  was  assumed  to  be  0.0  and  0.0058. 


Suppose  that  we  are  interested  in  selecting  a  design  with  the  best 
value  of  A/ug^  on  the  basis  of  maximum  response  A23}^>  Using  the  Continuous 
System  Modeling  Program  with  an  IBM  360/40  digital  computer,  the  relation¬ 
ship  between  the  maximum  response  and  the  nondlmenslonallzed  A  was  found 
as  shown  in  Figure  28.  On  this  basis,  the  design  engineer  could  choose  a 
value  of  A/ust  0.05,  which  results  In  a  maximum  displacement  of  0.2  of 
the  corresponding  static  displacement  Ug^  for  zero  damping  and  the  1940 
El  Centro  earthquake  excitation. 


This  simple  example  Is  given  herein  to  Illustrate  the  concepts  In¬ 
volved  in  the  direct  approach  method,  particularly  the  use  of  the  design- 
tree  In  seismic  design.  As  it  Is  pointed  out  in  the  text  of  this  report, 
the  seismic  design  of  building  structures  Is  a  complex  problem.  For 
example,  the  proper  choice  of  damping  coefficients  for  various  types  of 
buildings  Is  to  be  determined.  Consequently,  more  research  Is  required 
to  make  the  direct  approach  really  practical  in  seismic  design. 
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